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Preamble 

As in earlier reports, we will continue to break our effort into seven distinct units: 

• Atmospheric Correction Algorithm Development 

• Whitecap Correction Algorithm 

• In-water Radiance Distribution 

• Residual Instrument Polarization 

• Pre-launch Atmospheric Correction Validation 

• Detached Coccolith Algorithm and Post-launch Studies 

This separation has been logical thus far; however, as launch of AM-1 approaches, it must be rec- 
ognized that many of these activities will shift emphasis from algorithm development to validation. 
For example, the second, third, and fifth bullets will become almost totally validation-focussed 
activities in the post-launch era, providing the core of our experimental validation effort. Work 
under the first bullet will continue into the post-launch time frame, but will be driven in part by 
algorithm deficiencies revealed as a result of validation activities. We will continue to use this 
format for CY97. 
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Abstract 


Several significant accomplishments were made during the present reporting period. 

• We developed a new method for identifying the presence of absorbing 
aerosols and, simultaneously, performing atmospheric correction. The al- 
gorithm consists of optimizing the match between the top-of-atmosphere 
radiance spectrum and the result of models of both the ocean and aerosol 
optical properties. 

• We developed an algorithm for providing an accurate computation of the 
diffuse transmittance of the atmosphere given an aerosol model. A module 
for inclusion into the MODIS atmospheric-correction algorithm was com- 
pleted. 

• We acquired reflectance data for oceanic whitecaps during a cruise on the 
RV Ka’tmimoana in the Tropical Pacific (Manzanillo, Mexico to Honolulu, 
Hawaii). The reflectance spectrum of whitecaps was found to be similar to 
that for breaking waves in the surf zone measured by Frouin, Schwindling 
and Deschamps [1996]; however, the drop in augmented reflectance from 
670 to 860 nm was not as great, and the magnitude of the augmented 
reflectance was significantly less than expected [Gordon and Wang , 1994]. 

• We developed a method for the approximate correction for the effects of 
the MODIS polarization sensitivity. The correction, however, requires ad- 
equate characterization of the polarization sensitivity of MODIS prior to 
launch. 
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1. Atmospheric Correction Algorithm Development. 

a. Task Objectives: 

During CY 1996 there are four objectives under this task: 

(i) Complete development of an algorithm module for removing the effects of stratospheric 
aerosol and/or cirrus clouds from MODIS imagery over the oceans. 

(ii) Conduct research on the effects of strongly absorbing aerosols, and their vertical structure, 
on the existing atmospheric correction algorithm. Use the results of this research to develop a 
strategy for their removal. 

(Hi) Develop a detailed model of the diffuse transmittance of the atmosphere and the manner 
in which it is influenced by the angular distribution of subsurface upwelling spectral radiance. Add 
a module for this to the atmospheric correction algorithm. 

(iv) Investigate the effects of ignoring the polarization of the atmospheric light field on the 
performance of the proposed atmospheric correction algorithm. 

b. Work Accomplished: 

(i) This task was placed on hold to free time to accelerate Major Task #4 below in response 
to the higher- than-expected MODIS polarization sensitivity in some ocean bands. However, a 
paper submitted to Applied Optics describing our work in this area was accepted for publication 
(Appendix 2 in the Jan. to June 1996 Semi-Annual Report), and talks with MODIS Science Team 
Member B.-C. Gao were initiated to coordinate effort for cirrus cloud removal. 

(ii) As demonstrated our last Semi-Annual Report (Jan. to June 1996), strongly absorb- 
ing aerosols present a serious problem for atmospheric correction. The nature of the problem is 
two fold: (1) in contrast to weakly- absorbing aerosols, when the aerosol is strongly absorbing its 
distribution in altitude becomes very important; and (2) the technique of distinguishing aerosol 
type through examination of the spectral variation of the radiance in the near infrared, used by 
the MODIS atmospheric-correction algorithm, cannot distinguish between weakly- absorbing and 
strongly- absorbing aerosols. During this reporting period , we have tested a “spectral matching 
algorithm” that, although very slow, is capable of distinguishing between weakly- and strongly- 
absorbing aerosols. It is based on combining an atmospheric model with a water-leaving radiance 
model for the ocean, and effecting a variation of the relevant parameters until a satisfactory fit to 
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the MODIS top-of-atmosphere radiance is achieved. We believe that the algorithm is also capa- 
ble of functioning in the same manner when aerosol vertical structure is an additional parameter. 
(Note that vertical structure is only important when the aerosol is strongly absorbing.) A report 
describing this algorithm is attached as Appendix 1. 

(iii) The basic correction algorithm yields the product of the diffuse transmittance t and the 
water-leaving reflectance p w . However, t depends on the angular distribution of p w . If p w were 
uniform, t would be easy to compute, and this approximation has always been employed in the past. 
In a series of papers Morel and Gentili [Morel and Gentili , 1991; Morel and Gentili , 1993] studied 
theoretically the bidirectional effects as a function of the sun-viewing geometry and the pigment 
concentration. Their simulations suggest that, although the bidirectional effects nearly cancel in the 
estimation of the pigment concentration using radiance ratios, p w can depend significantly on the 
solar and viewing angles. (Our major task number 3, a study of the in-water radiance distribution, 
experimentally addresses this problem.) 

In this reporting period, we completed a study to understand the influence of bidirectional 
effects on the diffuse transmittance t. Through the use of the reciprocity principle, we were able to 
develop a simplified method of computing t, given the upward radiance distribution with direction 
just beneath the sea surface. We showed that the difference between t (the correct diffuse trans- 
mittance) and t * (the diffuse transmittance computed by assuming the subsurface upward radiance 
is uniform) is typically < 4%, and is a relatively weak function of the aerosol optical thickness. 
Thus, considering the error likely to result from the removed of the aerosol path radiance, it ap- 
pears that in the blue t can be replaced by £*, except in waters with low pigment concentrations, 
e.g., ;$ 0.5 mg/m 3 . A paper describing this work has been submitted to Applied Optics, and is 
attached here as Appendix 2. 

We have written a software module to accurately compute the diffuse transmittance, and 
produced the required lookup tables. This module has been delivered to R. Evans for integration 
into the atmospheric correction software. 

(iv) This task was completed as described in the last Semi-Annual Report (Jan. to June 1996). 

c. Data/ Analysis/Interpretation: See item b above. 

d. Anticipated Future Actions: 

(i) We must now implement our strategy for adding the cirrus cloud correction into the existing 
atmospheric correction algorithm. Specific issues include (1) the phase function to be used for 
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the cirrus clouds, (2) the details of making two passes through the correction algorithm, and (3) 
preparation of the required tables. These issues will be addressed during CY 1997 with the goal of 
having a complete implementation ready by the end of CY 1997. 

(ii) We will continue the study of the “spectral matching” algorithm with the goal of having 
an algorithm ready for implementation by the end of CY 1997. As our work has shown that 
a knowledge of the vertical distribution of the aerosol is critical, if it is strongly absorbing, we 
have procured a micro pulse lidar (MPL) system for use at sea on validation cruises, and from 
islands (likely Barbados or the Canary Islands) in the Saharan dust zone, to begin to compile the 
climatology of the vertical distribution required to adopt candidate distributions for use in this 
region. 

(iii) This task is nearly complete; however, we still need to develop a method to include the 
effects of the subsurface BRDF for low-pigment waters in the blue. 

(iv) Task completed. 

Additional tasks for CY97: 

(1) We will initiate a study to determine the efficacy of the present atmospheric 
correction algorithm on removal of the aerosol effect from the measurement 
of the fluorescence line height (MOD 20). 

(2) We will examine methods for efficiently including earth-curvature effects 
into the atmospheric correction algorithm [Ding and Gordon , 1994]. This 
will most likely be a modification of the look-up tables for the top-of-the- 
atmosphere contribution from Rayleigh scattering. 

(3) We will examine the necessity of implementing our out-of-band correction 
[Gordon, 1995] to MODIS. 

e. Problems/Corrective Actions: 

(i) None. 

(ii) None. 

(iii) None. 

(iv) None. 
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(v) None, 
f. Publications: 

H.R. Gordon, Atmospheric Correction of Ocean Color Imagery in the Earth Observing System Era, 
Journal of Geophysical Research, Atmospheres (Accepted). 

H.R. Gordon, T. Zhang, F. He, and K. Ding, Effects of stratospheric aerosols and thin cirrus clouds 
on atmospheric correction of ocean color imagery: Simulations, Applied Optics (In press). 

H. Yang and H.R. Gordon, Remote sensing of ocean color: Assessment of the water-leaving radiance 
bidirectional effects on the atmospheric diffuse transmittance, Applied Optics (Submitted). 
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2. Whitecap Correction Algorithm (with K.J. Voss). 

a. Task Objectives: 

As described in earlier reports, a whitecap radiometer system has been built and tested to 
provide a database for developing and validating the whitecap correction algorithm, as well as for 
providing an estimation of the whitecap contribution to the water-leaving radiance during the post- 
launch validation phase. The database includes spectral information as well as variables associated 
with the formation and occurrence of whitecaps such as wind speed and air/sea temperature. 

b. Work Accomplished: 

The basic design, calibration procedure, operation, and data analysis methods for the whitecap 
radiometer are described in a paper prepared for submission to the Jour . Atm . Ocean. Tech . and 
attached here as Appendix 3. 

Prom 29 March to 18 April 1996 the whitecap radiometer system was deployed on the NO A A 
ship RV Malcolm Baldrige on a cruise from Miami to a test location in the Gulf of Mexico, ap- 
proximately 70 miles off shore from Cedar Key (Florida) in the Apalachicola Bay. The location 
provided relatively warm waters (16°-17°C) with a number of cold fronts moving off the mainland. 
These fronts usually lasted a couple of days bringing strong winds (sometimes as high as 18 m/s) 
and lowering the air temperature to about 12°C. The occurrence of an unstable atmosphere and 
good winds provided an interesting spectral whitewater data set. From 1 to 13 November 1996, 
the whitecap radiometer was operated on a cruise from Manzanillo, Mexico to Honolulu, Hawaii. 
This cruise provided whitecap data under conditions of steady winds (the trades) of essentially 
unlimited duration and fetch. The analysis of these two data sets is provided in a paper prepared 
for submission to Jour. Geophys. Res. and attached here as Appendix 4. 

c. Data/ Analysis/Interpretation 

There have been three significant results from our whitecap research: (1) our measurements 
confirm the spectral fall off of the augmented reflectance in the NIR, although the reduction at 
865 compared to 670 nm was not as large as observed in the surf zone [Frouin, Schwindling and 
Deschamps , 1996]; (2) whitecaps show significant nonlambertian effects, particularly at large solar 
zenith angles; and (3) the augmented reflectance of whitecaps appears to be about one-forth that 
predicted by recent models [Gordon and Wang , 1994; Koepke , 1984]. The details of these results 
are provided in Appendix 4. 

d. Anticipated Future Actions: 
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As the basic objectives of this task have been realized, it is being suspended until the val- 
idation phase. Karl Moore, the post doctoral associated who was responsible for the operation 
of the instrument and the data analysis, has accepted a position at the Scripps Institution of 
Oceanography. 

e. Problems/ Corrective Actions: None 

f. Publications: 

K.D. Moore, K.J. Voss, and H.R. Gordon, Whitecaps: Spectral reflectance in the open ocean and 
their contribution to water-leaving radiance, Ocean Optics XIII, Halifax, Nova Scotia, October 
22 - 25 , 1996 . 

K.D. Moore, K.J. Voss, and H.R. Gordon, Spectral reflectance of whitecaps: Instrumentation, 
calibration, and performance in coastal waters, Jour. Atmos. Ocean. Tech . (Submitted). 

K.D. Moore, K.J. Voss, and H.R. Gordon, Spectral reflectance of whitecaps: Fractional coverage 
and the augmented spectral reflectance contribution to the water-leaving radiance, Jour . Geophys. 
Res . (Submitted). 
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3. In-water Radiance Distribution (with K.J. Voss). 

a. Task Objectives: 

The main objective in this task is to obtain upwelling radiance distribution data at sea for 
a variety of solar zenith angles to understand how the water-leaving radiance varies with viewing 
angle and sun single. 

b. Work accomplished: 

We acquired upwelling radiance distribution data with the RADS camera system during a 
cruise with Dennis Clark during November. It was very windy during this cruise, but data was 
acquired in the configuration expected to be used during validation cruises. 

c. Data/ Analysis/Interpretation: none 

We are currently performing the post-cruise calibrations of this instrument, and will be reduc- 
ing the cruise data when this is completed. 

d. Anticipated future actions: 

We will be completing the calibrations and reducing data acquired during the cruise. We are 
also planning on making some minor changes to the instrument to improve its operation at sea. 
Finally we will be acquiring more data during another cruise with Dennis Clark this spring. 

e. Problems/Corrective actions: None. 

f. Publications: None. 
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4. Residual Instrument Polarization. 

a* Task Objectives: 

The basic question here is: if the MODIS responds to the state of polarization state of the 
incident radiance, given the polarization-sensitivity characteristics of the sensor, how much will this 
degrade the performance of the algorithm for atmospheric correction, and how can we correct for 
these effects? 

b. Work Accomplished: 

We have developed a formalism [ Gordon , 1988] which provides the framework for removal of 
instrumented polariz at ion- sensitivity effects, and a method for approximately correcting for them. 
The correction method is presented in a paper submitted for publication in Applied Optics , and 
attached as Appendix 5. The main difficulty that we see now in correcting for the polarization- 
sensitivity effects is the requirement of an adequate pre-launch characterization of the polarization 
sensitivity. 

c. Data/ Analysis/Interpretation: 

See Appendix 5. 

d. Anticipated Future Actions: 

Although this task is now basically complete (a correction algorithm has been developed), 
we still need to prepare a module for including the polarization-sensitivity correction algorithm 
in the MODIS atmospheric correction algorithm. Also, as operation of the polarization-sensitivity 
correction algorithm requires an adequate pre-launch characterization of the polarization sensitivity, 
we will continue to work with MCST and SBRS to insure that proper characterization is realized. 

e. Problems/ Corrective Actions: None 

f. Publications: 

H.R. Gordon, T. Du, and T. Zhang, Atmospheric correction of ocean color sensors: Analysis of the 
effects of residual instrument polarization sensitivity, Applied Optics (Submitted). 
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5. Pre-launch Atmospheric Correction Validation (with K.J. Voss). 

a. Task Objectives: 

The long-term objectives of this task are four-fold: 

(i) First, we need to study aerosol optical properties over the ocean in order to verify the 
applicability of the aerosol models used in the atmospheric correction algorithm. Effecting this 
requires obtaining long-term time series in typical maritime environments. This will be achieved 
using a CIMEL sun/sky radiometer that can be operated in a remote environment and send data 
back to the laboratory via a satellite link. These are similar the radiometers used by B. Holben in 
the AERONET Network. Thus, the first objective is to deploy a CIMEL Automatic Sun Tracking 
Photometer in a suitable location for studying the optical properties of aerosols over the ocean. 

(ii) Second, we must be able to measure the aerosol optical properties from a ship during 
the initialization/calibration/validation cruises. The CIMEL-type instrumentation cannot be used 
(due to the motion of the ship) for this purpose. The required instrumentation consists of an all-sky 
camera (which can measure the entire sky radiance, with the exception of the solar aureole region) 
from a moving ship, am aureole camera (specifically designed for ship use) and a hamd-held sun 
photometer. We have a suitable sky camera and sun photometer and must construct am aureole 
camera. Our objective for this calendar year is (1) to assemble, characterize and calibrate the solar 
aureole camera system, (2) to develop data acquisition software, and (3) to test the system. 

In the case of strongly- absorbing aerosols, we have shown that knowledge of the aerosol vertical 
structure is critical [ Gordon , 1996]. Thus, we need to be able to measure the vertical distribution 
of aerosols during validation exercises. This can be accomplished with ship-borne LIDAR. We have 
to procure a LIDAR system and modify it for ship operation. 

(iii) The third objective is to determine how accurately the radiance at the top of the atmo- 
sphere can be determined based on measurements of sky radiance and aerosol optical thickness at 
the sea surface. This requires a critical examination of the effect of radiative transfer on “vicarious” 
calibration exercises. 

(iv) The forth objective is to utilize data from other sensors that have achieved orbit (OCTS, 
POLDER, MSX), or are expected to achieve orbit (SeaWiFS) prior to the launch of MODIS, to 
validate and fine-tune the correction algorithm. 

b. Work Accomplished: 
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(i) During the first part of this period we were operating the CIMEL in it’s location in the 
Dry Tortugas. In October this instrument was removed for recalibration. At the same time the 
AERONET network, run by B. Hobren, decided to upgrade the CIMEL instruments with more 
stable interference filters and small hardware changes. We will be reinstalling this instrument when 
it is returned from NASA. 

(ii) The sky camera system and aureole system was used on three cruises during this pe- 
riod. The first cruise we participated in occurred during the NASA TARFOX experiment between 
Bermuda and New York. This cruise took place on a cruise liner, and we participated in all three 
weeks of the cruise effort. During this time we made measurements with the sky camera (including 
polarization), aureole camera system (to measure the sky radiance near the sun), and a hand held 
sunphotometer. In general the weather was not as good as expected (a hurricane when through the 
area during the cruise) but there were several clear periods during the cruise. The second and third 
cruises were with Dennis Clark off of Hawaii (during September and November). In the first of these 
cruises we obtained data with all three systems, but in the third Dennis Clarks group provided the 
sun photometer data. In addition to participating on the cruises we performed calibrations of all 
the systems pre and post cruise. We are currently rewriting the data reduction programs for the 
sky radiance distribution system to allow data reduction to take place on the cruise, to speed this 
process. In addition this is giving an opportunity to review this process and to perform tests of this 
data. We have reduced the aureole data from the first two cruises, and are currently evaluating 
this data. 

To address the problem of vertical distribution of aerosols we have acquired a Micro Pulse 
Lidar from SSEI. This system was delivered at the very end of the period (December 16th), so we 
have not had a chance to do anything with the system yet. 

(iii) As described in our last Semi-Annual Report (Jan. to June 1996), we have completed a 
study of the accuracy with which one can compute the radiance at the top of the atmosphere from 
sky radiance measurements made at the sea surface. The results suggest that the bulk of the error 
is governed by the uncertainty in the sky radiance measurements. Furthermore, it was shown that 
the largest error in the radiative transfer process was that due to the use of scalar radiative transfer 
theory, and that improvement would require the use of vector theory, and thus, measurement 
of the the polarization of the sky radiance. We have started to analyze the use of polarization 
measurements at the surface and it appears that when polarization is included in the sky radiance 
inversion algorithm, the radiative transfer error can be made very small. Furthermore, we are now 
examining the extent to which the full linear polarization of the top-of-atmosphere radiance can 
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be deduced from surface measurements. This may be very important for validating the pre-launch 
polarization- sensitivity characterization of MODIS. 

(iv) We have been in contact with personnel involved with SeaWiFS, OCTS, and MSX to 
acquire data formats, and satellite data from these instruments to assess the validity of the at- 
mospheric correction algorithm. We have procured an SGI R10000 Workstation (same chip set 
and operating system as used by MODIS SDST). This will provide the necessary image processing 
capability for the pre- and post-launch era. 

c. Data/ Analysis/Interpretation: 

(i) The CIMEL instrument was at NASA for most of this period for recalibration and refitting 
of the interference filters. Thus there has been little data analysis. 

(ii) We have been working to reduce the data from the sky camera and aureole camera acquired 
this period but have not finished this process. Thus there is little to report in terms of analysis or 
interpretation. We also do not have any data from the LIDAR yet. 

d. Anticipated Future Actions: 

(i) We will be reinstalling the CIMEL in the Dry Tortugas at the first opportunity after its 
return from NASA. We are also working on a better method of acquiring the data through NASA. 
This will enable us to look at the sky radiance data in a more timely manner. 

(ii) We will finish the data reduction work with the sky camera system in the next quarter. 
We are also reworking portions of this system to allow more automation of the data collection, and 
fix minor problems which developed during the last cruise (specifically overheating of the system 
computer and corrosion on the computer backplane). The reduced aureole data will be merged 
with the sky radiance data to provide a complete sky radiance distribution during this next period. 
We will also finish reducing all of the aureole data during this next period. 

Much of our effort in the next severed months will be spent learning the new MPL Lidar 
system, and making the modifications necessary to use the system at sea. We would like to have 
this system ready for operation during the next cruise opportunity, if possible ready for a field tried 
in February. 

(iii) We will continue development of a sky radiance inversion edgorithm that utilizes the 
f ull vector radiative treinsfer equation. This should remove the letrgest radiative tremsfer error in 
predicting the top-of-atmosphere radieince from the bottom-of-atmosphere radieince. 
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(iv) We will continue to try to obtain ocean color data from other sensors to assess the 
correction algorithm. Now, OCTS appears to be the most likely candidate. We will install R. Evans’ 
implementation of our atmospheric correction algorithm on the R10000 computer to facilitate fine- 
tuning the algorithm. 

e. Problems/corrective actions: None. 

f. Publications: 

H.R. Gordon and T. Zhang, How well can radiance reflected from the ocean- atmosphere system be 
predicted from measurements at the sea surface?, Applied Optics 35, 6527-6543 (1996). 

D.K. Clark, H.R. Gordon, K.J. Voss, Y. Ge, W. Broenkow, and C. Trees, Validation of Atmospheric 
Correction over the Oceans, Jour. Geophys. Res. (Accepted). 

K.J. Voss and Y. Liu, Polarized radiance distribution measurements of skylight: Part 1, system 
description and characterization, Applied Optics (Submitted). 
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6. Detached Coccolith Algorithm and Post Launch Studies (W.M. Balch). 

During the second half of 1996, we participated in one cruise to the Gulf of Maine, a well known 
region for mesoscale coccolithophore blooms. We measured coccolithophore abundance, produc- 
tion and optical properties. We also completed a flow cytometer experiment in which we sorted 
field-derived calcite particles to measure their calcite specific scattering coefficients. A thorough 
understanding of the relationship between calcite abundance and light scatter, in coccolith-rich and 
coccolith-poor regions, will provide the basis for a generic suspended calcite algorithm to be used 
with MODIS data. 

a. Task Objectives: 

The algorithm for retrieval of the detached coccolith concentration from the coccolithophorid, 
E. huxleyi is described in detail in our ATBD. The key is quantification of the backscattering coeffi- 
cient of the detached coccoliths. Our earlier studies focussed on laboratory cultures to understand 
factors affecting the calcite-specific backscattering coefficient. As with algorithms for chlorophyll, 
and primary productivity, the natural variance between growth related parameters and optical 
properties needs to be understood before the accuracy of the algorithm can be determined. To this 
end, the objectives of our coccolith studies have been to define the effect of growth rate on: 

(1) the rate that coccoliths detach from cells (this is also a function of turbulence 
and physical shear); 

(2) rates of coccolith production; 

(3) morphology of coccoliths; and 

(4) volume scattering and backscatter of coccoliths. 

For perspective on the directions of our work, we provide an overview of our previous activities. 
During 1995, we focussed on all of the above objectives using chemotstat cultures (in which algal 
growth rate is precisely controlled). During the latter half of 1995, our work focused on shipboard 
measurements of suspended calcite and estimates of optical backscattering as validation of the 
laboratory measurements. We participated on two month-long cruises to the Arabian sea, measuring 
coccolithophore abundance, production, and optical properties. During the first half of 1996, we 
focused again on objectives 2 and 4, during two Gulf of Maine cruises, one in March and one in 
June. During the second half of 1996, we participated on another cruise to the Gulf of Maine, and 
further addressed objectives 2 and 4. 
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b. Work Accomplished: 

We have continued data analysis for the Arabian Sea cruises as well as the first two Gulf of 
Maine cruises of 1996. The current state of the data are as follows: 

(1) We completed our third cruise of 1996 to the Gulf of Maine in November. During 
this cruise, we sampled for toted and calcite-dependent backscattering (continu- 
ously), suspended calcite concentrations, calcification rates, chlorophyll concen- 
trations, and coccolithophore and coccolith counts. 

(2) We completed our last set of flow cytometer experiments in which we sorted 
field-derived calcite particles into vials of pure seawater, measured their volume 
scattering functions, and measured their calcite concentrations. 

(3) Suspended calcite samples from all Arabian Sea cruises and two cruises to the 
Gulf of Maine (March and June) have been run in the graphite furnace atomic 
absorption spectrometer at the University of Maine. We now only have samples 
from our most recent cruise in November. 

(4) The 400 cell and coccolith counts from the southwest Monsoon cruise (summer 
’95) have been completed and the data are being hand-entered into spreadsheets 
at this time. We are currently working on the Arabian Sea Samples from the ’95 
intermonsoon cruise. 

(5) All calcification data from the Arabian Sea cruise have been processed to units of 
gC m“ 3 d“ l and integrated over the water column at each station and they have 
been processed into complete sections. 

(6) The Arabian Sea biooptical data is now processed for temperature, salinity, pH, 
fluorescence and backscatter (with and without calcite) averaged over each kilo- 
meter of all trips. All calibration checks have been done for the underway data. 
Hydrographic plots have been made in which light scattering is plotted as a func- 
tion of temperature and salinity and the optical and pigment data superimposed. 

(7) We have examined the relationship between the calcite-dependent backscattering 
(6£) and the concentration of suspended calcite or concentration of detached coc- 
coliths for our previous work in the Straits of Florida, Arabian Sea, and North 
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Atlantic. We have also analyzed these data relevant to our flow cytometer re- 
sults. This is of major relevance to our MODIS algorithm efforts. Besides actu- 
ally checking our algorithm, it allows us to define the accuracy and precision of 
the algorithm. This is exceedingly important for subsequent interpretation. Our 
results suggest that for the Gordon reflectance model, the algorithm will have a 
precision of ± 25,000 coccoliths/mL (or in terms of carbon equivalents, ± 5 fig C 
as calcite per liter). 

c. Data/Analysis/Interpretation: 

As expected, calcite-dependent backscattering was low in the Gulf of Maine during March, but 
it was still measurable. Typically, calcite scattering accounted for 5-10% of toted backscattering. 
A very different picture was observed in June. Acid-labile scattering increased to 30-40% of total 
backscattering in Wilkinson Basin, a stratified basin in the middle of the Gulf of Medne. Acid-labile 
scattering dropped over Georges Bank as the predominant populations were diatoms, and values 
increased again in the Northeast Channel, similar to previous blooms that we have observed. The 
observations are consistent with the calcite being produced in the more stable Wilkinson Basin 
with subsequent advection around the NE flank of Georges Bank. We completed another Gulf of 
Maine cruise at the end of October, and again took the flow-through system. Much of these data 
are still in the process of being worked-up. 

Calcification measurements from the March 1996 cruise were remarkably high, given that this 
was at the beginning of the Spring bloom. We were finding >10% of the carbon being fixed into 
coccoliths. This also explains the relatively high fraction of calcite-dependent light scattering seen 
during this time. Calcification rates are still being processed for the June cruise but preliminary 
results suggest that they were quite high. 

d. Anticipated Future Actions: 

Work in CY97 will address several areas: 

(1) Processing of the suspended calcite samples from the November ’96 cruise. 

(2) Data will be collated from the Arabian Sea and Gulf of Maine cruises in order 
to calculate turnover of the calcite particles. This cam only occur after step 1 is 
completed. 

(3) We are continuing analysis of our latest flow cytometer results from the fall ’96 
experiment. The suspended calcite samples from that last experiment will be run. 
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(4) The underway data from the Gulf of Maine will be merged with our calibration 
measurements (vicarious calibrations are periodically made at sea and these data 
must be processed to verify whether instrument calibrations changed). 

(5) Hydrographic plots of the Gulf of Maine data will be made in which light scat- 
tering and chlorophyll are plotted in temperature salinity space. 

e. Problems /Corrective Actions: None 

f. Publications: 

W. M. Balch, J. J. Fritz, and E. Fernandez, Decoupling of calcification and photosynthesis in 
the coccolithophore Emiliania huxleyi under steady-state light-limited growth. Marine Ecology 
Progress Series , 142 87-97 (1996). 

E. Fernandez, J. J. Fritz and W. M. Balch, Growth- dependent chemical composition of the coccol- 
ithophorid Emiliania huxleyi in light-limited chemostats, J . Exp. Mar. Biol. Ecol. (In press). 

J. J. Fritz and W. M. Balch, A coccolith detachment rate determined from chemostat cultures of 
the coccolithophore Emiliania huxleyi, J. Exp. Mar. Biol. Ecol. (In press). 

K. J. Voss, W. M. Balch, and K. A. Kilpatrick. Scattering and attenuation properties of Emiliania 
huxleyi cells and their detached coccoliths, Limnol. Oceanogr. (Submitted). 
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8. Other Activities. 

The bulk of the PFs effort during this reporting period was focused on five activities. The 
first was the preparation of revisions of the Water-leaving Radiance and Coccolith concentration 
ATBD’s. These revisions were completed and delivered to the EOS Senior Project Scientist on 
August 15, 1996. The second was a detailed revision of the MOCEAN Validation Plan on behalf 
of the Group. This occurred after attending the MODIS Ocean Group (MOCEAN) Meeting in 
July. The third was participation in the MCST audit of progress toward the Level IB algorithm 
in the VIS/NIR bands on September 5, 1996. The forth was participation in the MODIS ATBD 
review November 20, 1996. The fifth was an intensive effort to assess the effects of the larger-than- 
expected MODIS polarization sensitivity on MODIS products, and to understand the anomalous 
polarization- sensitivity characterization tests. 
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9. Publications, submissions, and abstracts for CY 96. 

H.R. Gordon, Atmospheric Correction of Ocean Color Imagery in the Earth Observing System Era, 
Journal of Geophysical Research , Atmospheres (Accepted). 

H.R. Gordon, T. Zhang, F. He, and K. Ding, Effects of stratospheric aerosols and thin cirrus clouds 
on atmospheric correction of ocean color imagery: Simulations, Applied Optics (In press). 

H.R. Gordon and T. Zhang, How well can radiance reflected from the ocean- atmosphere system be 
predicted from measurements at the sea surface?, Applied Optics , 35, 6527-6543 (1996). 

T. Zhang and H.R. Gordon, Columnar aerosol properties over oceans by combining surface and 
aircraft measurements: sensitivity analysis, Applied Optics (Accepted). 

D. K. Clark, H.R. Gordon, K.J. Voss, Y. Ge, W. Broenkow, and C. Trees, Validation of Atmospheric 
Correction over the Oceans, Jour. Geophys. Res. (Accepted). 

H. Yang and H.R. Gordon, Remote sensing of ocean color: Assessment of the water-leaving radiance 
bidirectional effects on the atmospheric diffuse transmittance, Applied Optics (Submitted). 

H.R. Gordon, T. Du, and T. Zhang, Atmospheric Correction of Ocean Color Sensors: Analysis of 
the Effects of Residual Instrument Polarization Sensitivity, Applied Optics (Submitted). 

K.D, Moore, K.J. Voss, and H.R. Gordon, Whitecaps: Spectral reflectance in the open ocean and 
their contribution to water-leaving radiance. Ocean Optics XHI, Halifax, Nova Scotia, October 
22-25, 1996. 

K.D. Moore, K.J. Voss, and H.R. Gordon, Spectral reflectance of whitecaps: Instrumentation, 
calibration, and performance in coasted waters, Jour. Atmos. Ocean. Tech. (Submitted). 

K.D. Moore, K.J. Voss, and H.R. Gordon, Spectral reflectance of whitecaps: Fractional coverage 
and the augmented spectral reflectance contribution to the water-leaving radiance, Jour. Geophys. 
Res. (Submitted). 

K.J. Voss and Y. Liu, Polarized radiance distribution measurements of skylight: Part 1, system 
description and characterization, Applied Optics (Submitted). 

W. M. Balch and K. A. Kilpatrick, Calcification rates in the equatorial Pacific along 140°W, Deep 
Sea Research, 43 971-993 (1996). 
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Appendix 1 


A method for atmospheric correction in the 
presence of strongly absorbing aerosols 



1. Introduction 


The Coastal Zone Color Scanner (CZCS) demonstrated the feasibility of measuring marine 
phytoplankton concentrations from earth-orbiting sensors. 1 ’ 2 Based the success of the CZCS, sev- 
eral similar instruments with higher radiometric sensitivity and a larger number of spectral bands, 
e.g., the sea-viewing wide- field- of- view sensor (SeaWiFS), 3 the moderate resolution imaging spec- 
troradiometer (MODIS), 4 etc., will be launched in the near future. These ocean color instruments 
will actually measure the chlorophyll a concentration in the water as a surrogate for the phyto- 
plankton concentration. In fact, the CZCS measured the sum of the concentrations of chlorophyll 
a and its degradation product phaeophytin a. This sum was referred to as the pigment concentra- 
tion, C. Phytoplankton pigments have a broad absorption maximum in the blue (~ 435 nm) and a 
broad absorption minim um in the green (~ 565 nm), and the CZCS algorithms derived C from the 
ratio of the radiances backseat tered out of the water (the water- leaving radiance,!^) near these 
two wavelengths. 5 ’ 6 Typically, L w is at most 10% of the total radiance, L t , exiting the top of the 
atmosphere (TOA) in the blue and < 5% in the green. Therefore, it is necessary to extract L w 
from L t to derive C. This process is called atmospheric correction. 

The atmospheric correction algorithm developed for CZCS 6 ~ 10 is not sufficiently accurate for 
the new generation of sensors with higher radiometric sensitivity. Atmospheric correction of these 
sensors requires incorporation of multiple-scattering effects. Gordon and Wang 11 ’ 12 developed such 
an algorithm for SeaWiFS, and found that the multiple- scattering effects depended on the physical 
and chemical properties of the aerosol (size distribution and refractive index). Therefore, incorpo- 
ration of multiple scattering into atmospheric correction required the introduction of aerosol models 
in the algorithm. 

The Gordon and Wang algorithm is very simple to describe. The spectral variation in L t in 
the near infrared (NIR) spectral region where L w % 0 is used to provide information concerning 
the aerosol, as L t there is due principally to Rayleigh scattering (which is known) and to aerosol 
scattering. The Rayleigh scattering component is then removed, and the resulting spectral variation 
is compared to that produced by a set of candidate aerosol models in order to determine which 
two models of the candidate set are most appropriate. These models sire then used to estimate 
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the multiple-scattering effects. Gordon 12 has shown that this algorithm can provide L w with the 
desired accuracy as long as the aerosol is weakly absorbing (more accurately, the aerosol must be 
weakly absorbing and it must follow the relationship between size distribution and refractive index 
that is implicitly implied in the choice of the candidate aerosol models). Unfortunately, strongly- 
absorbing aerosols, e.g., aerosols from anthropogenic urban sources or mineral dust transported from 
desert areas to the ocean, can possess size distributions similar to the weakly-absorbing aerosols 
typically present over the oceans. As the spectral variation of aerosol scattering depends mostly 
on the aerosol size distribution, and only weakly on the index of refraction, the spectral variation 
of scattering in the NIR is not sufficient to distinguish between weakly- and strongly-absorbing 
aerosols. Furthermore, in the case of mineral dust an additional complication arises: the dust is 
colored, i.e., its absorption is a function of wavelength. 13 ’ 14 Even if it were possible to estimate 
the absorption characteristics of mineral dust aerosol in the NIR, one would still not know the 
extent of absorption in the visible. This is a particularly serious problem, as regions contaminated 
by mineral dust are often highly productive and thus important from a biogeochemical point of 
view. In fact, dust deposition may actually be providing nutrients that enable the phytoplankton 
to bloom. 15 

The difficulty in detecting the presence of strongly absorbing aerosols is that the effects of 
absorption become evident only in the multiple scattering regime. In the single scattering regime, 
the reflectance of the aerosol is proportional to the product of the single scattering albedo (a>o) 
and the aerosol optical thickness (r 0 ), i.e., at small r a there is no way to distinguish nonabsorbing 
aerosols (u;o = 1) with a given r a from absorbing aerosols (u>o < 1) and a larger r a . Retrieval 
of information concerning aerosol absorption requires multiple scattering; however, this multiple 
scattering need not be aerosol multiple scattering — when a low concentration of aerosols exists 
in the presence of strong Rayleigh scattering, e.g., in the blue, multiple Rayleigh scattering can 
increase the length of photon paths through the aerosol and enhance the chance of absorption. Also, 
if distributed vertically in the atmospheric column, the aerosol can reduce the Rayleigh- scattering 
component, which is otherwise large in the blue. Thus, the possibility of inferring aerosol absorption 
is increased as one progresses from the NIR into the visible, but unfortunately L w is not known there 
(that is why atmospheric correction is required in the first place). The inescapable conclusion is 
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that the SeaWiFS algorithm 1 1 must fail when the aerosol is strongly absorbing unless the candidate 
aerosol models sire restricted to those with similarly strong absorption properties. 12 Unless one 
could use a global aerosol climatology, i.e., knowledge of the aerosol expected to be present in a 
given region and time, to restrict the candidate aerosol models to the appropriate set, an alternate 
approach is needed. 

In addition to atmospheric correction, there is compelling interest in studying the global dis- 
tribution and transport of aerosols because of their role in climate forcing and biogeochemical 
cycles. 16 ’ 17 Furthermore, not only is the aerosol concentration required, it is also important to 
know their absorption properties to understand their climatic effects. There has been continuing 
interest in measuring aerosol concentration from earth-orbiting sensors. 18-27 Over the oceans these 
sensors generally utilize spectral bands for which the ocean can be assumed to be black (L w = 0) 
or at least to have constant reflectance. By analogy to the atmospheric correction problem above, 
estimation of aerosol absorption properties from space fails for these sensors. The one exception is 
the retrieval of spatial distributions of strongly- absorbing aerosols using the Total Ozone Mapping 
Spectrometer (TOMS) measurements in the ultraviolet. 28 

In this paper, we describe an alternate approach to the problem of estimating oceanic bio- 
physical properties, and aerosol physical- chemical properties, using space-borne ocean color sensors. 
The approach is to utilize all of the spectral bands of the sensor. This insures sufficient multiple- 
scattering (Rayleigh scattering in the blue) to enable identification of aerosol absorption. In order 
to separate the effects of aerosols from radiance backscattered from beneath the sea surface (L t ), an 
ocean color model, in which the reflectance is related to the phytoplankton pigment concentration 
and the scattering properties of the phytoplankton and their associated detrital material, is used. 
As with the SeaWiFS algorithm, severed candidate aerosols models are employed: nonabsorbing, 
weakly absorbing, and strongly absorbing. Through a systematic variation of the candidate model 
aerosols, phytoplankton scattering, C, and r a , a “best” fit to simulated spectral L t data is obtained. 
It is found that the algorithm can successfully discriminate between weakly- and strongly-absorbing 
aerosols, and can provide estimates of C, r G , and u; 0 with an accuracy that is nearly independent 
of For consistency with earlier work, 12 we specifically examine a situation that would be 
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encountered off the U.S. East Coast in summer: polluted continental air transported by the winds 
to the Middle Atlantic Bight. However, this situation is used only as an example to demonstrate 
the approach, we believe that it could be applied to ocean regions subjected to mineral dust as well 
as aerosol resulting from biomass burning. 

We begin with a discussion of the approach and the modeling of the various quantities required 
for implementation. Next, we test the efficacy of the algorithm using simulated SeaWiFS 3 data. 
Recall SeaWiFS has 8 spectral bands centered at 412, 443, 490, 510, 555, 670, 765, and 865 nm. 
Finally, we examine the degradation of the performance of the algorithm in the presence of L t - 
measurement error. 

2. The algorithm approach and implementation 

Neglecting influence of direct sun glitter and whitecaps, the total upwelling reflectance exit- 
ing the top of the atmosphere p f (A) consists of the following components: 11 ’ 12 the pure Rayleigh 
(molecular) scattering contribution p r (A), the pure aerosol scattering contribution p a (A), the con- 
tribution due to the interaction effect between air molecules and aerosols p ra (A), and the desired 
water-leaving contribution t(A)p tu (A), i.e., 

^t(^) = Pr(A) + Po(^) + />ra(-M + (1) 

where t(A) is the diffuse transmittance of the atmosphere. From the satellite image, we have the 
spectrum of the upwelling reflectance p t (A). As p r (A) depends only on the surface atmospheric 
pressure, 9 ’ 29 it is not difficult to remove the pure Rayleigh scattering contribution p r (A) from the 
total reflectance p*(A), 


[p«(^) - Pr(A)] — [p a (A) + Pra(A)] + *(A)pu,(A). (2) 

The known reflectance spectrum of [p t (A) - p r (A)] consists of two parts which are hard to separate, 
the water-leaving reflectance term t(A)pu,(A) and the aerosol contribution [p a (A) + p ra (A)] (which 
includes the interaction term between aerosols and air molecules). The goal of atmospheric correc- 
tion is to retrieve the water-leaving reflectance p„,(A) from the known reflectance [p t (A) - p P (A)]. 
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Because of the high spatial and temporal variability of the physical, chemical, and optical prop- 
erties of aerosols, it is difficult to estimate the aerosol contribution [p a (A) + Pra{ A)] to the total 
up welling reflect ance. 

The basic assumption of the proposed algorithm is that for each aerosol and pigment con- 
centration there is a unique and distinctive spectrum characteristic of its upwelling reflectances 
[/>a(A) + Pra(A)] and [t(A)p^(A)] . In a given sun- viewing geometry, similar (or close) spectra 
to [p t (A) — p r (A)] can only be obtained from the atmosphere-ocean system by a combination of 
aerosols having similar optical properties to the actual aerosol, and a pigment concentration similar 
to that actually present in the ocean. That is, when we estimate [t(A)p 1i; (A)] and [p a (A) + Pra(A)] 
separately and form [p a (A) + Pra^)]* and [t(A)p w (A)] r , where here and henceforth the primes will 
refer to computed or tried estimates, the computed reflectance spectrum 

[MA) - />r(A)]' = [*,(A) + /»,.(A)]' + [*(A)p.( A)] 

will fit the true reflectance spectrum [p*(A) - p r (A)] in the visible and near infrared only if both the 
computed water-leaving reflectance [t(A)/> u; (A)] / and the computed aerosol contribution [p a (A) + 
Pra(^)] # fit their true values individually. In order to implement this idea, we need to be able 
to obtain estimates of [^(A)/3 1i; (A)] / and jp 0 (A) + pra(^)] , « We now describe how the use of this 
assumption is accomplished. 

2.A The water component: tp w 

The prediction of the water-leaving reflectance, /^(A), is effected using the semi-empirical bio- 
optical radiance model, developed by Gordon et a/. 30 for Case 1 waters 5 ’ 31 , i.e., waters for which the 
optical properties sire controlled by the water itself and by the concentration of phytoplankton and 
their decay products. As we use the pigment concentration C as a surrogate for the phytoplankton 
concentration, one would expect that the absorption and scattering properties of the particles 
would depend only on C\ however, it is found for such waters that for a given C the total scattering 
coefficient varies by roughly a factor of two. 5 Thus, a second scattering-related parameter 6° 
which ranges from 0.12 to 0.45 m* 1 with a mean value of 0.30 m -1 (when C has units of mg/m 3 ) 
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Figure 1 Spectrum of normalized water-leaving reflectance [p w (A)]^ for pig- 
ment concentrations of 0.1, 0.5, and 1.0 mg/m 3 . For each pigment concen- 
tration, from the upper to the lower curves, the values of coefficient 6° are 
0.45, 0.30, and 0.12 m -1 , respectively. 

is introduced. Gordon et al. 30 found that by using such a model, the water-leaving radiance 
dependence on C in Case 1 waters could be explained. Similar results were also obtained by 
Bricaud and Morel. 32 This bio-optical ocean color model actually provides the normalized water- 
leaving reflectance, 6 ’ 12 [pu,(A)]jv, defined by 

[/>w(A)]jv = p w ex p[+(r P /2 + r Ol )/cos0 o ], 

where r r and t 0z are the Rayleigh and Ozone optical thicknesses of the atmosphere, respectively, 
and 6 0 is the solar zenith angle. Sample spectra of [p„(A)]iv as a function of C and b° are shown 
in Figure 1. Clearly, the normalized water-leaving reflect ance is very sensitive to the pigment 
concentration C for short wavelengths (412 and 443 nm) and small pigment concentration (C < 0.4 
mg/m 3 ). For longer wavelengths (555 and 670 nm) or large pigment concentrations ( C > 0.8 
mg/m 3 ), [/>u,(A)]jv does not depend significantly on the pigment concentration C. It is taken to be 
zero at 765 and 865 nm. 
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It should be noted that the Gordon et a/. 30 reflectance model above does not take into ac- 
count the bidirectional effects of the sub-surface upwelled spectral radiance, i.e., it assumes that 
the upwelling radiance beneath the sea surface is totally diffuse. Morel and co-workers 33 ~ 36 have 
demonstrated that this is not the case; however, as Morel and Gentili 36 have shown, bidirectional 
effects can be easily introduced into the model and described as a function of C . 

Once p w (A) is determined, it is necessary to propagate it to the top of the atmosphere (TOA). 
As mentioned earlier, this is accomplished using the diffuse transmittance f(A). Tanre et a/. 37 
and Gordon et a/. 6 provided simple expressions for t( A) that include the effects of both aerosol 
and Rayleigh scattering. Later, Yang and Gordon 38 provided a detailed analysis of t( A) based on 
precise computations. They showed that (1) bidirectional effects play a role in t only in the blue 
and only at low C, (2) aerosols play a significant effect on t only if they are strongly absorbing, (3) 
t is independent of the aerosol vertical structure even if the aerosol is strongly absorbing, and (4) 
given an aerosol model it is simple to predict the correct value of t for any aerosol concentration 
and viewing geometry. The value of t can be computed precisely given C (to provide bidirectional 
effects), an aerosol model (to provide the aerosol properties), and the aerosol optical thickness r a (A) 
(to provide the aerosol concentration). However, for the purposes of this paper, we will approximate 
t by assuming it is independent of the aerosol. In this case, t(A) is given by 

t{ A) = exp — (t p (A)/2 + t 0z )/ cos0„] , 

where 6 V is the angle between the zenith and a line from the sensor to the pixel under consideration. 
Thus, the simulated t(A)p w (A) is given by 

*(A)p»(A) = [/»«,( A)] w exp [-(T P (A)/2 + r 0a ) (— + — ^-)j. (3) 

L \ COS U o COS U v ) J 

2.B The aerosol component: p a + p ra 

Gordon and Wang 11 have shown that the multiple scattering effects in p a + p ra depend sig- 
nificantly on the physical and chemical properties of the aerosol, i.e., their size distribution and 
refractive index. Thus aerosol models had to be introduced to incorporate multiple scattering effects 
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in atmospheric correction. Similarly, aerosol models are also required to retrieve aerosol proper- 
ties from space observations. 39 ’ 40 Gordon and Wang 11 used aerosol models that were developed 
by Shettle and Fenn 41 for LOWTRAN-6 42 . These models consist of particles distributed in size 
according to combinations of two log-normal distributions, and are described in detail in Ref. 12. 
Briefly, four models at four different relative humidities are used here. These are the Maritime (M), 
the Coastal (C), the Tropospheric (T), and the Urban (U). The relative humidities used are 50%, 
70%, 90%, and 99%. We denote a particular model by a letter and a number, e.g., M99 refers to 
the Maritime model at 99% relative humidity. There is an increasing amount of absorption as one 
progresses through M,C,T, to U. For example, at 865 nm the aerosol single-scattering albedo, u; 0 , is 
0.9934, 0.9884, and 0.9528, respectively, for the Maritime, Coastal, and Tropospheric models (RH 
= 80%), while in contrast, c^o = 0.7481 for the Urban model. Here, the Urban model is intended 
to represent strongly absorbing aerosols that might be present over the oceans near areas with 
considerable urban pollution, e.g., the Middle Atlantic Bight off the U.S. East Coast in summer. 

We employ these sixteen aerosol models as candidates to test the algorithm. For a two-layer 
atmosphere, with the aerosol confined to the bottom layer, the scalar radiative transfer equation 
(polarization ignored) was solved for each of the 16 candidate aerosol models (M, C, T, U aerosols 
with RH = 50%, 70%, 90%, 99%) with eight values of r a (A) in the range of 0.05 to 0.8 at each 
wavelength A, for solar zenith angle 0 O = 0° to 80° in increments of 2.5°, and for 33 different viewing 
zenith angles with 0 V in the range of 0° to 90°. As it is difficult to have this large computational set 
of values of [p a (A) + Pra(^)] available for image processing for all of the aerosol candidates, aerosol 
optical thicknesses, sun- viewing geometries (0 O , 0 V ), and spectral bands; in a manner similar to 

the Gordon and Wang algorithm, 12 lookup tables (LUTs) are used to provide [p a (A) + p ra { A)] . In 
the Gordon and Wang algorithm the lookup tables related [p a (A) + Pra(^)] to p aj (A), the single 
scattered aerosol reflectance. Equivalently, we relate the term [^ a (^) + Pra(^)] to the aerosol optical 
thickness r a ; our simulated values of [p a (A) + /? ra (A)] are fit to 

[pM + *.(A)] = o( A) r a + 6(A) r\ + c(A) r* + d(A) r* (4) 

using least-squares. To further reduce storage, coefficients a(A), 6(A), c(A) and d{ A) were expanded 
in Fourier series in the azimuth view-angle <p v , and only the Fourier coefficients were stored in 
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the LUTs. Samples of the fit of [/> a (A) + /? ra (A)j to the aerosol optical thickness are presented in 
Figure 2 for a sun-viewing geometry with 6 0 = 60°, 8 V « 46°, and <(> v % 93°. This geometry has 
relatively large errors in the fits compared to the other geometries used here. It can be seen that 
the method of using Eq. (4), and Fourier expanding its coefficients, is appropriate for computing 
[/> a (A) + Pra(A)] for any candidate aerosol model, at any sun-viewing geometry, aerosol optical 
thickness, and wavelength. The largest fitting errors occurred at the largest wavelength (A = 865 
nm) and small aerosol optical thickness (r a % 0.1). They were of the order of 1% to 2%. 



Figure 2(a) Curve fits of [p a (A) + Pra(A)] vs. 
r a for aerosol models of M70, C70, T70 and 
U70 at A = 443 nm with $q = 60.0° t 0 v = 
45.92°, and <fi = 93.49°. 



Figure 2(b) Curve fits of [pa(A) + p Ta (A)] vs. 
r a for aerosol models of M70, C70, T70 and 
U70 at A = 865 nm with 0o = 6O.O°,0* = 
45.92°, and </> = 93.49°. 


2.C The algorithm implementation 


We experimented with severed approaches for implementing the algorithm. The one we found 
most effective is summarized as follows: 


First, for the given sun-viewing geometry (0o, &v, 4>v), we vary the value of aerosol optical thick- 
ness at 865 nm, r a (865), for each candidate aerosol model to provide the aerosol component 
[p a (A) + Pra(A )] ' ■ We then vary the pigment concentration C and the scattering-related coef- 
ficient b° to provide the water-leaving reflectance [t(A)p u ,(A)]\ These yield a trial upwelling 
reflectance [/>t(A) - /^(A)]* at each of N bands of the ocean color sensor. 
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Second, we compute the percent deviation 6* of this simulated spectrum [p*(A) — p r (A)] / from 
the measured true spectrum [p*(A) - p r (A)j for each test set (A,r a , (7, 6°)', where A labels 
the candidate aerosol model. The percent deviation 6(A, r a , (7, 6°)' is defined in the sense of 
root-mean-squares, 


*(A,r a ,C,fc 0 )' = 


1 N 
1 V 

- ^r(Af)] - [p t (A0 - />r(A,)]'l 

N - 1^1 

i=l 

[*(*<)- MA*)] J 


( 5 ) 


Third, we sort the deviations S(A , r a , <7, 6°)' to find 10 best sets of (A, r a , <7, 6°)' which yield 10 
smallest percent deviations £(A, r a , <7, b 0 )'. 


Fourth, as the correct aerosol model is unlikely to be identical to one of the candidates, we 
assume the characteristics of aerosols and pigment concentration can be adequately described 
by these 10 best sets of parameters (A, r a , (7, 6°)'. The retrieved single scattering albedo u> f 0 
and optical thickness of the aerosols, and 6 o/ and the pigment concentration C l in the ocean 
are then computed by averaging u>q, t^, b and C 1 over the 10 best sets of (A, r a , <7, 6°)'. 

3. The algorithm’s performance 

In this section, we examine the performance of the spectrum-matching algorithm by applying it 
to SeaWiFS. 3 The sun- viewing geometries sire taken as those used in Refs. 11 and 12: viewing at the 
center of the scan (viewing zenith angle 0 V % 1°) for solar zenith angle $q = 20°, 40° and 60°, and 
viewing at the edge of the scan ( 0 V % 45°) near the perpendicular plane for &o = 0°, 20°, 40° and 60°. 
These cover much of the the range of sun-viewing geometries available to SeaWiFS. Pseudo data 
are provided by solving the scalar radiative transfer equation for a two-layer atmosphere system 
with a specified aerosol confined in the lower layer. The pseudo water-leaving reflectance /^(A) 
was provided for 6° = 0.30 m -1 (the mean value for Case 1 waters) and pigment concentrations 
C — 0.1, 0.5, and 1.0 mg/m 3 . 

The algorithm attempts to match the pseudo data spectrum of pt( A) by varying the testing 
aerosol model among its 16 candidates ( N A = 16), the aerosol optical thickness r a (865) from 0.01 


11 



to 0.40 in increments of 0.01 ( N r = 40), the pigment concentration from 0.05 to 1.50 mg/m 3 in 
increments of 0.05 mg/m 3 ( N c = 30), and finally, b° from 0.12 to 0.45 m -1 in increments of 0.03 m -1 
( N b = 12). The total number of elements in the test set (A, r tt , C, b 0 )' is N = N A x N r x N c X N b = 
16 x 40 x 30 x 12 = 230400. 

For the first test of the algorithm we examined cases in which the aerosol optical properties 
of the pseudo atmosphere were included in the candidate aerosol models, i.e., the aerosol optical 
properties in the atmosphere system were taken from M70, C70, T70 and U70. The optical thickness 
at 865 nm was taken to be r a (865) = 0.1, 0.2 or 0.3. The main purpose of this was to test the 
code for implementation of the algorithm. In all cases, for the best set (smallest S') the correct 
aerosol model and the correct values of the parameters were chosen. In fact, S' for the correct 
set was a small fraction of a percent and ~ 10 to 30 times smaller them the second best set. The 
residual error was due to small errors in the LUTs. Even the averages over the ten best sets were 
excellent, providing close values of r 0 (865), 6°, and C. As the aerosol single scattering albedo o> 0 
is a weak function of wavelength A, we use its the retrieved value at 865 nm, u> 0 (865), (averaged 
over the ten best sets) as an indication of the algorithm’s ability to distinguish between weakly- 
and strongly-absorbing aerosols. The derived values of u/ 0 (865) showed that weakly and strongly 
absorbing aerosols are easily recognized by the algorithm. 

As it is not likely for the aerosols in the atmosphere will have exactly the same optical properties 
as any one of candidate aerosol models, we tested more realistic cases in which the aerosol models 
were similar to, but not the same as, any of the sixteen candidate aerosol models. Following 
Gordon 12 the aerosol models M80, C80, T80 and U80 (Shettle and Fenn models with relative 
humidity 80%) were chosen for this purpose. We begin by describing the results obtained from 
averaging the parameters from the sets with the ten smallest values of S'. The averaged ^(865), 
7^(865), and C 1 for a given geometry are taken to be the retrieved values of these parameters. 
To estimate the performance on a more global scale, i.e., for all sun-viewing geometries, we then 
average over all seven sun- viewing geometries and compute the mean and the standard deviation 
in the retrieved parameter values. The mean values of retrieved aerosol single scattering albedo 
u>o(865) are provided in Table 1 for aerosol models M80, C80, T80 and U80. It can be seen from 
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Table 1 that the retrieved results for u> 0 are very good for each of four pseudo aerosol models. 
Large percent deviations (the standard deviation over the seven geometries divided by the mean) 
in the range of 3% to about 8%, are encountered for the strongly absorbing U80 aerosol model, 
nevertheless the algorithm can distinguish between the weakly absorbing aerosols (M80, C80, T80) 
and the strongly absorbing aerosol (U80) without difficulty. 

Since the ultimate goal of ocean color remote sensing is to estimate the phytoplankton pig- 
ment concentration, we now examine the retrieval of C using the algorithm. Table 2 presents the 
mean values of retrieved C\ which are averaged over seven sun- viewing geometries and also over 
the four test aerosol models, M80, C80, T80, and U80 (28 cases in all). It can be observed that 
the retrieved results of pigment concentration are reasonable for all three tested aerosol optical 
thicknesses [r a (865) = 0.1, 0.2, and 0.3] and all three pigment concentrations [C = 0.1, 0.5, and 
1.0 mg/m 3 ]. For the small pigment concentration, C = 0.1 mg/m 3 , or for small aerosol optical 
thickness, r a (865) = 0.1, the spectrum-matching algorithm still works very well. With an increase 
in either pigment concentration or aerosol concentration, the percent deviations and percent er- 
rors in the retrieved C l become larger. For comparison, Table 3 provides similar results for the 
weakly-absorbing aerosols only, using the Gordon and Wang correction algorithm. 12 Note that the 
present algorithm behaves as well as the Gordon and Wang algorithm, even when strongly- absorbing 
aerosols are included. Had strongly-absorbing aerosols been included in Table 3, the results would 
have been significantly poorer, e.g., in some cases it would have been impossible to even compute 
C because one or both of the required p w ' s would likely be negative. 

Detailed retrievals of the pigment concentration C are tabulated in Table 4, which gives per- 
centages of cases with relative error in the retrieved, \AC\jCy less than 5%, 10%, 20%, and 30%, 
respectively. For the smallest pigment concentration, in all of the 84 cases ex amin ed (three aerosol 
optical thicknesses, four aerosol models, and seven sun- viewing geometries), |AC|/C < 20%, and 
even < 5% for about 90% of the cases. For a pigment concentration of C = 0.5 mg/m 3 , and for 
small aerosol optical thickness, r a (865) = 0.1, all individual simulations have | AC\/C < 30%, while 
for r a (865) = 0.2 there are about 89% of the cases having |AC|/C < 30%, and for r a (865) = 0.3 
about 75% of the cases having \AC\/C < 30%. At the highest pigment concentration examined 
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(1.0 mg/m 3 ), for small aerosol optical thickness the algorithm still performs very well with only 
one case with |AC|/C > 30%. However, as r n (865) increases, the errors become larger, and about 
71% of the cases have |AC|/C < 30% for r a (865) = 0.2, and only about 60% for r a (865) = 0.3. 

The algorithm clearly works better in estimating C for smaller pigment concentrations. This is 
explained by the relationship between the (normalized) upwelling water-leaving reflectance and the 
pigment concentration (see Figures 1 and 2). As we stated in Section 2, the water-leaving reflectance 
depends strongly on pigment concentration C when C is small (C < 0.4 mg/m 3 ). At small C, a 
small change in C (0.05 mg/m 3 in the algorithm) will result in a significant change in the upwelling 
water-leaving reflectance. But when pigment concentration is as large as about 1.0 mg/m 3 , the 
upwelling water-leaving reflect ance is only a weak function of pigment concentration, and a small 
change of 0.05 mg/m 3 in C is not going to result in any significant change in the upwelling water- 
leaving reflectance. This causes the large percent deviation in C when the algorithm is applied to 
large pigment concentrations. 

Mean values of the retrieved aerosol optical thickness r^(865) over the seven sun-viewing ge- 
ometries and four testing aerosol models (M80, C80, T80 and U80) are presented in Table 5. 
The mean values are close to their corresponding “true” aerosol optical thicknesses, the percent 
deviations are ranged from about 6% to about 11%. 

Figure 3 provides samples of the three best sets (A,r a ,C, b°) determined by the algorithm for 
aerosol models of M80, C80, T80 and U80 with r 0 (865) = 0.2 and C = 0.5 mg/m 3 , for a single 
sun-viewing geometry ( 9 0 = 20°, 8 V = 45.92°, <f> v = 90°). As the pseudo aerosol models (M80, C80, 
T80 and U80) are similar to the candidate models used in the algorithm, (M, C, T, and U with RH 
= 50%, 70%, 90% and 99%), but are not identical to any of 16 candidates, there is no correct aerosol 
model for the algorithm to choose to match the upwelling reflectance [/> a (A) -f p ra (A) -M(A)/j w (A)]. 
The figure shows that the aerosol models which have similar optical properties to those of pseudo 
aerosol models are selected first by the algorithm. Even though there are some errors caused by 
picking the incorrect aerosol models, the pigment concentration chosen by the algorithm is close 
to its true value (0.5 mg/m 3 in Figure 3). For the four pseudo aerosol models tested, the percent 
deviations of the best match for whole spectrum 6(A y r a , C, 6°)' varied from about 0.7% to 1.1%. 
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Unlike the case when the test aerosol was one of the candidates, there was no significant increase 
in 6* from the best set to the second best set. For the best 10 sets, the largest percent deviation 
for C = 0.5 mg/m 3 is about 1.5%. 


M80. t 8 =0.2, C=0.S. b°=0.30, 0^=0.993 C80. ^=0.2, C=0.5, b°=0.30. .988 




Figure 3(a) Reflectance spectrum matching 
for aerosol model M80 and pigment concen- 
tration C = 0.50 mg/m 3 with sun- viewing 
geometry of $ 0 = 20.0°, $ v = 45.92°, and 
<t> = 90.0°. 


Figure 3(b) Reflectance spectrum matching 
for aerosol model C80 and pigment concen- 
tration C — 0.50 mg/m 3 with sun-viewing 
geometry of 0q = 20.0° , 0 V = 45.92° , and 
<f> = 90.0°. 


m tj=0.2, CM). 5, b°=0.30. (1)0=0.953 



Figure 3(c) Reflectance spectrum matching 
for aerosol model T80 and pigment concen- 
tration C — 0.50 mg/m 3 with sun- viewing 
geometry of 0q = 20.0°, 0 V = 45.92° , and 
<t> = 90.0°. 



Figure 3(d) Reflectance spectrum matching 
for aerosol model U80 and pigment concen- 
tration C = 0.50 mg/m 3 with sun- viewing 
geometry of $o = 20.0°, $ v = 45.92°, and 
<p = 90.0°. 


Prom these tests of the algorithm we conclude that it can detect the presence of strongly 
absorbing aerosols successfully. Whenever the optical properties of aerosol in the atmosphere 
are the same (or very close) to that of any of the 16 candidate aerosol models employed in the 
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algorithm, the retrieved pigment concentration C will be excellent, meeting the requirements of 
SeaWiFS and MODIS. If the optical properties of the aerosol in the atmosphere are s imil ar to 
that of any candidate aerosol model, the retrieved results for the pigment concentration C will still 
be good in the presence of small pigment concentration or small aerosol optical thickness. When 
both large aerosol optical thickness [r 0 (865) around 0.2 to 0.3] and large pigment concentration 
[C around 0.5 to 1.0 mg/m 3 ] are present in the atmosphere-ocean system, the performance of the 
algorithm will be degraded. Still, for the worst scenario examined here, r a (865) = 0.3 and C = 1.0 
mg/m 3 , about 60% of the individual simulations have |AC|/<7 < 30%. 

These simulations suggest that the success of the algorithm depends on the appropriateness of 
candidate aerosol models and the bio-optical model employed in the algorit hm Although we will 
not be able to know the percent errors in the retrieved aerosol single scattering albedo Wq( 865) and 
pigment concentration C' in processing a satellite image, we can compute S(A,T a ,C,b°Y over all 
bands which are used for ocean color remote sensing. Generally, small percent deviations lead to 
small percent errors, so that the percent deviations for retrieved u>g(865) and C 1 cam also be used 
to assess the performance of the algorithm. 

4. Further tests 

Addition tests of this algorithm remain to be completed. As mentioned in Section 1, when the 
aerosol is strongly absorbing, p t depends significantly on the aerosol’s vertical distribution. In the 
simulations presented here, the correct vertical distribution was assumed, i.e., the pseudo data were 
created using the same vertical distribution as was assumed for the candidate aerosol models. We 
believe that the vertical distribution can be introduced into the algorit hm simply as new candidate 
aerosol models, e.g., the U70 model with all of the aerosol in the marine boundary layer and the 
U70 aerosol model with the aerosol uniformly mixed in the entire atmosphere, would represent 
distinct candidate aerosol models. This hypothesis will be tested. Also, we need to examine the 
sensitivity of the algorithm to errors in p t — p r . 
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Table 1: Mean values of retrieved u/q( 865) for the seven sun- viewing geometries 
and each of four aerosol models (M80, C80, T80, U80). The standard 
deviation divided by the mean is listed in parenthesis. 



C: mg/m 3 



1.000 

M80 

w 0 = 0.993 

t„(865) = 0.100 

0.992 (0.43%) 

0.996 (0.14%) 

0.997 (0.10%) 

t„(865) = 0.200 

0.995 (0.10%) 

0.995 (0.10%) 

0.996 (0.05%) 

r a (865) = 0.300 

0.996 (0.05%) 

0.996 (0.06%) 

0.996 (0.10%) 



C : mg/m 3 

0.100 

0.500 

1.000 

■ 

r a (865) = 0.100 

0.980 (0.97%) 

0.972 (1.60%) 

0.965 (2.26%) 



0.988 (0.31%) 


r a (865) = 0.300 

0.987 (0.25%) 

0.987 (0.27%) 

0.987 (0.37%) 



C: mg/m 3 

0.100 

0.500 

1.000 

■Kiyi 


0.952 (0.78%) 

0.935 (0.13%) 

0.935 (0.00%) 

r a (865) = 0.200 




r a (865) = 0.300 

0.945 (0.31%) 

0.934 (0.04%) 

0.945 (1.82%) 



C : mg/m 3 

0.100 

0.500 

1.000 

U80 

u>o = 0.748 

r a (865) = 0.100 

0.793 (4.21%) 

0.761 (4.51%) 

0.769 (3.36%) 

r a (865) = 0.200 

0.730 (4.76%) 

0.750 (8.25%) 

0.712 (7.44%) 

r a (865) = 0.300 


0.784 (2.56%) 

0.699 (7.34%) 
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Table 2: Mean values of retrieved C' for seven sun- viewing geometries 
and four aerosol models (M80, C80, T80, U80). The standard 
deviation divided by the mean is listed in parenthesis. 


C: mg/m 3 

0.100 

0.500 

1.000 

r a (865) = 0.100 

0.100 (1.90%) 

0.528 (7.30%) 

1.098 (11.7%) 

T a (865) = 0.200 

0.101 (4.56%) 

0.547 (15.2%) 

0.982 (23.9%) 

r a (865) = 0.300 

0.101 (4.42%) 

0.612 (24.1%) 

0.947 (31.3%) 


Table 3: Mean values of retrieved C' for seven sun- viewing geometries 
and three aerosol models (M80, C80, T80). derived using 
the Gordon and Wang 11 algorithm. The standard deviation 
divided by the mean is listed in parenthesis. 


C : mg/m 3 

0.10 

0.47 

0.91 

t 0 (865) = 0.100 

0.101 (1.6%) 

0.466 ( 3.4%) 

0.912 ( 9.1%) 

r a (865) = 0.200 

0.100 (3.1%) 

0.470 ( 4.7%) 

0.940 (12.8%) 

r a (865) = 0.300 

0.098 (5.5%) 

0.493 (15.3%) 

0.936 (25.3%) 
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Table 4: Percentage of retrieved pigment concentration C' within certain 
error limits for aerosol models of M80, C80, T80, and U80. 



A C/C : 

< 5% 

< 10% 

< 20% 

< 30% 

C = 0.1 mg/m 3 

r„(865) = 0.100 

96% 

100% 

100% 

100% 

r a (865) = 0.200 

89% 

92% 


100% 

r a (865) = 0.300 

89% 

92% 

100% 

100% 



AC/C : 

< 5% 

< 10% 

< 20% 

< 30% 

C = 0.5 mg/m 3 

r a (865) = 0.100 

32% 

75% 

92% 

100% 

r a (865) = 0.200 

32% 

50% 

75% 

89% 

r a (865) = 0.300 

25% 

35% 

64% 

75% 



AC/C : 

< 5% 

< 10% 

< 20% 

< 30% 

C = 1.0 mg/m 3 

r a (865) = 0.100 

21% 

50% 

71% 

96% 

r a (865) = 0.200 

28% 

39% 

57% 

71% 

r a (865) = 0.300 

14% 

25% 

50% 

60% 


Table 5: Mean values of retrieved r'(865) for seven sun-viewing geometries 
and four aerosol models (M80, C80, T80, U80). The standard 
deviation divided by the mean is listed in parenthesis. 


C: mg/m 3 




r a (865) = 0.100 

0.102 (8.94%) 

0.101 (10.2%) 

0.102 (11.6%) 

r a (865) = 0.200 

0.201 (6.31%) 

0.199 (8.53%) 

0.199 (8.20%) 

r a (865) = 0.300 

0.300 (6.22%) 

0.294 (8.79%) 

0.300 (9.93%) 
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Abstract 


Two factors influence the diffuse transmittance (t) of water-leaving radiance ( L w ) to the top 
of the atmosphere: the angular distribution of upwelling radiance beneath the sea surface (Z u ); 
and the concentration and optical properties of aerosols in the atmosphere. We examine these 
factors and show (1) that the error in L w induced by assuming L u is uniform, i.e., in treating 
the subsurface reflectance by the water-body as lambertian, is significant in comparison to the 
other errors expected in L w only at low phytoplankton concentration and then only in the blue 
region of the spectrum, (2) that when radiance ratios are used in biophysical algorithms the affect 
of the uniform-T u approximation is even smaller, and (3) the diffuse transmittance is a strong 
function of the absorption properties of aerosol in the atmosphere, but is nearly independent of 
their vertical distribution. An avenue for introducing accurate computation of the unif orm- L u 
diffuse transmittance into atmospheric correction algorithms is presented. In an appendix the 
reciprocity principle for a medium in which the refractive index is a continuous function of position 
is derived. 
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1. Introduction 


The feasibility of measuring marine phytoplankton concentrations from earth-orbiting sensors 
was demonstrated by the proof-of-concept Coastal Zone Color Scanner (CZCS) 1 ’ 2 mission. Based on 
the CZCS experience, several similar instruments with a larger number of spectral bands and higher 
radiometric sensitivity are being prepared for launch, e.g., the sea-viewing wide-field-of-view sensor 
(SeaWiFS), 3 the moderate resolution imaging spectroradiometer (MODIS), 4 etc. The surrogate 
for measurement of the phytoplankton concentration is the concentration of the photosynthetic 
pigment chlorophyll a within the water (actually, within the plants). As chlorophyll a has a broad 
strong absorption in the blue (~ 435 nm) and a minimum of absorption in the green (~ 565 nm), 
the concentration within the water can be estimated from the solar radiance backscattered out of 
the water near these wavelengths, 5 ’ 8 the water-leaving radiance. Unfortunately, the water-leaving 
radiance typically comprises at most 10% of the total radiance exiting the top of the atmosphere 
(TOA). Briefly, the radiance exiting the top of the atmosphere in a spectral band centered at A », 
can be written 

Lt{ AO = LotH'ri K) + (1) 


where Lother{^i) represents the contribution to the radiance from all sources except the water- 
leaving radiance propogated to the TOA, (A*). Sources of L other include scattering of solar 
radiation in the atmosphere, specular reflection of scattered and unscattered radiation from the di- 
rect solar beam by the sea surface, and diffuse reflection from oceanic whitecaps. The atmospheric 
correction algorithm of Gordon and Wang 7 ’ 8 estimates Lother and removes it from the total ra- 
diance, thereby obtaining the water-leaving radiance transmitted to the top of the atmosphere, 
££,(A,). The water-leaving radiance at the top of the atmosphere is related to the water-leaving 
radiance at the bottom of the atmosphere (usually just called the water-leaving radiance) through 
the diffuse transmittance, f(£„), i.e., 


t{L) = 


Kjjv ) 
L w (i v ) 


where is a unit vector directed from the sea surface to the sensor, and L u ((„) is the water-leaving 
radiance just above the surface. In addition to the attenuation of L w (£„) along the path from the 
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surface to the sensor, the diffuse transmittance also accounts for its augmentation by scattering of 
into the direction i.e., atmospheric scattering from £ to £„. 

The water-leaving radiance L w can assume a range of values. Table 1 provides values of L w 
at 443 nm and 550 nm with the sun near the zenith for a range of pigment concentrations 9 (C, 
the sum of the concentrations of chlorophyll a and its degradation product phaeophytin a) for 
Case 1 waters, 5 i.e., waters for which the optical properties are controlled by the water itself, the 
phytoplankton, and the phytoplankton decay products. The goal of the atmospheric correction 
algorithm is the derivation of L w from L t with an uncertainty of < 5% at 443 nm for very clear 
oceanic waters, e.g., the Sargasso Sea in summer (for which C ~ 0.03 mg/m 3 ). As t is of order 
unity, the residual error in L w (after removal of L other from L t ) should be <0.1 mW/cm 2 /rai 
Sr. The Gordon and Wang 7 algorithm is capable of this performance at 443 nm. Clearly, at higher 
pigment concentrations the relative error in L w at 443 nm will be higher, e.g., ~ 25% for C ~ 1 
mg/m 3 (Table 1). Furthermore, the error in the removal of Lother at 550 nm is about ? - 7 that 
at 443 nm, 10 i.e., ~ 0.03 mW / cm 2 /zm Sr. Thus, even in the Sargasso Sea in Summer, the relative 
error in L w at 550 nm will be ~ 10%. 

Table 1: Water-leaving radiance at 443 and 550 nm 
as a function of the pigment concentration. 


c 

(mg/m 3 ) 

I w (443) L w ( 550) 
(mW/cm 2 /im Sr) 

0.03 

1.95-2.20 

0.28-0.30 

0.10 

1.35-1.60 

0.28-0.30 

0.47 

0.40-0.75 

0.28-0.40 

0.91 

0.30-0.50 

0.24-0.50 


To retrieve the water-leaving radiance from I® , we need to calculate the diffuse transmittance. 
However, as we shall see, the diffuse transmittance itself is a function of the water-leaving radiance, 
and we generally need to know the angular distribution of the water-leaving radiance in order to 
calculate the diffuse transmittance. The angular distribution has been shown to possess significant 
bidirectional structure; 11 13 however, in the computation of t it has always been assumed that the 
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water- leaving radiance is nearly uniform (independent of viewing direction), i.e., that the deviation 
of the actual radiance distribution from uniform causes only negligible error in the computation of 
t. The mean objective of this paper is to evaluate the validity of this assumption. As we focus on 
sensors with a spatial resolution of ^ 1 km designed for the open ocean, where the typical scale of 
variability is a few km, we assume that L w is constant over the scene. Thus, the adjacency effect 
of the atmosphere 14 16 is ignored. We compare the diffuse transmittance calculated with several 
realistic water leaving radiance distributions, including one from actual measurements, to that 
with a uniform radiance distribution. The results reveal that errors in the retrieved water leaving 
radiance caused by making the uniform approximation are significant compared to the error in 
L other only in the blue and only at low pigment concentration. We then examine the behavior of t 
as a function of the concentration and optical properties of the aerosol present in the atmosphere. 


2. Computational procedure 

For the purposes of this paper, we assume the atmosphere is divided into two layers, a 
molecular-scattering layer on the top and an aerosol layer at the bottom. The ocean surface is 
assumed to be flat. 


The computation of t is straigthforward: the radiative transfer equation (RTE) can be solved 
for the radiance exiting the top of the atmosphere ( L ® ) with the correct upwelling radiance distri- 
bution Z u (£) incident just beneath the sea surface. Accounting for the transfer of radiance across 
the air- water interface, 


t(i) = 


Kii) 

Lw( o’ 


where 

iM) = 

Tf(£') is the Fresnel transmittance of the interface for radiance incident from below in the direction 
and m w is the refractive index of water. £ and £' are related by Snell’s law. This is the 
direct approach to finding t(£). However, rather than using this approach, we choose to solve the 
reciprocal problem and use the reciprocity principle (Appendix 1) to derive t. 
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In the reciprocal problem the extraterristrial solar beam is incident on the top of the atmo- 
sphere. Let Fq be the extraterrestrial solar irradiance, £o a unit vector in the direction of propaga- 
tion of the solar beam, and Lr ( £) the resulting radiance propagating downward just beneath the 
sea surface in the direction f . Then, it is shown in Appendix 1 that 


*(-&)) 


i 

Fo^o • uo|T/(£o) 



\i*n\L R tf) 


Lu tii 
M-€S) 




( 2 ) 


where L u ( — £) is the upward radiance distribution incident just beneath the sea surface for which 
we want t, £' 0 and £o are related by Snell’s law, and Qj indicates the integral is to be evaluated over 
all downward If L u (i) is uniform, this becomes 


n-io) = 


fol^o • no\Tf(£o) 


f \Z*h\Ln(i)dSl(i) = 


FrUo) 

-folCo • *o\T f (£o) 


( 3 ) 


where Er(£ o) is the downward irradiance just beneath the surface in the reciprocal problem. 


The relative error in the retrieval of water-leaving radiance caused in making the unif orm 
radiance approximation is 


AZ. w (-^q) _ £^(-£o) ~ -£u>(— &>) _ *(~£o) ~ **(—&)) 
L w {- io) £ w (~io) t m (~io) 

where £^(-£o) is computed from L“ (~£ 0 ) using t*(-£o)- Thus, 


L w (-l o) 


1 

E R (io) 



|£«n|Ijt(£) 


M-i'o) 


dSl(i), 


( 4 ) 


( 5 ) 


For the development of the atmospheric correction algorithm, Gordon and Wang carried out 
extensive simulations (~ 33,000) of the reciprocal problem for a variety of aerosol models, aerosol 
optical thicknesses, and £o's. According to Eq. (5), we can use these existing simulations to compute 
the error for any radiance distribution, L u (£), without having to carry out further radiative transfer 
simulations. This is the reason for using the reciprocial approach developed here. 

3. Computation of the error in L w 


In order to compute the error in L w using Eq. (5) we need the radiance backscattered toward the 
surface in the water. Precise computation of this requires coupling the ocean and the atmosphere 
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in a multiple scattering computation. However, because the radiance reflected out of the ocean 
is small (reflectance <0.04), the probability of photons being reflected out of the ocean twice 
is negligible. Therefore, the coupling is not really required and the ocean and atmosphere can be 
treated separately. 17 To further simplify the computations, we shall not consider multiple scattering 
in the water. Rather, we will employ the quasi-single scattering approximation 10 ’ 19 (QSSA) to 
provide given the scattering phase function for the medium. For given optical properties of 
the ocean, the QSSA will provide an upwelling radiance distribution beneath the surface that is 
identical in angular distribution to that produced in single scattering, and as such, will represent 
the upper limit to the departure of L u from a uniform distribution for a given phase function. The 
reason for this is that multiple scattering will always produce a smoother radiance distribution than 
single scattering. 

We introduce a cartesian coordinate system at the sea surface with the z-axis pointed in the 
downward direction (the same direction as n at the surface). In this system we describe the direction 
of £ by the angle 0 measured from the z-axis, and the azimuth angle <f> measured from the z-axis. If 
£ is directed toward increasing depth, 9 < 90°. Thus, the direction of the solar beam in the water 
is specified by 9q and <f>' 0} and for radiance in the water propagating upward toward the sea surface, 
0 > 90. We take <£o = 0, so photons that are exactly backscattered from the ocean- atmosphere 
system, i.e., scattered in the direction -£ 0) have <f> = 180°. 

Considering only the refracted direct solar beam in the water, in the QSSA the upwelling 
radiance distribution just beneath the sea surface is given by 18 ’ 19 

Lu (M) = C(u>0,P)F 0 T(d 0 ) COS /° P(0) . , (6) 

cos 9 0 - cos 9 

where cos 9 < 0 and C is a constant depending mostly on the value of the single- scattering albedo 
u/o, but also weakly on P. P(0) is the scattering phase function of the medium for a scattering 
angle 0, given by 

cos 0 = cos 9 cos 9g + sin 9 sin 9' 0 cos(<f> — <f >' 0 ). 
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When diffuse skylight refracted into the water in the direction £' is considered as well, an additional 
term 


<>o ,P)f Ltf) *»(('). 

«/{'«n>0 COS U — COS v 


where £<*(£') is the sky radiance transmitted through the air-sea interface, and 


cos 0 = cos 9 cos 9* + sin 9 sin 9 1 cos(0 - 0'), 


must be added to the right-hand-side of Eq. (6). 


In general the scattering phase function for the ocean, P(0), will be a combination of that due 
to molecular scattering by the water itself (Rayleigh scattering), Ph( 0), and that resulting from 
scattering by the suspended particles P p (0), i.e., 

p(0)= MM 0) + 6 P Pp( 0 ) ' 

+ b p 

where b r and b p are the Rayleigh and particle scattering coefficients, respectively. In contrast to 
forward scattering angles for which P p (0) > JP r (0), P p and P r are comparable for © > 80°, so 
Rayleigh scattering plays an important role in the deter min ation of P(0) for 0 > 80°. This is seen 
in Figure 1 which compares P T with P p measured by Petzold 20 in San Diego Harbor. Petzold’s San 
Diego Harbor phase function is often taken in ocean radiative transfer simulations as characteristic 
of particles in ocean water. 21-23 

The limits for P(0) are Pr{&) ( b p < 6 r ) and P p (0) ( b p > 6 r ). We shall estimate the error 
in L w given by Eq. (5) in these two limits. In general, the error will fall between these two limits. 
Figure 2 provides the upward in- water radiance distribution L u (— £) due to the direct solar beam 
alone for these two limits for 9o = 40° and 60°. In these figures, x — 9 = 0 implies photons are 
traveling toward the zenith, and for <f> = 180° (x - 0 ~ 0) and 9 = 0q the photons are exactly 
backscattered. Since photons with x — 9 > the critical angle (~ 48°) cannot escape the ocean 
when the surface is flat, the radiance distributions are truncated at x — 9 = 50°. Also, they are 
normalized to unity at their maximum value. In the Rayleigh scattering limit, an observer looking 
into the water with 9 = 9' 0 would observe maximum radiance for 0 = 180°, i.e., with the sun at the 
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observer’s back. In contrast, in the pure particle limit the maximum would be along 0 = 0. As we 
shall see, these radiance distributions lead to different values for the diffuse transmittance t. 

We have applied Eq. (5) and the QSSA to compute A L w / L w , the error in the recovered water- 
leaving radiance made by assuming that L u (9,<t>) is unif orm (constant). To compute and 

£/i(£o) h* Eq. (5) we assume as in Ref. 7 that the atmosphere can be approximated as two layers 
with pure molecular scattering in the upper layer and pure aerosol scattering in the lower layer. We 
use the Shettle and Fenn 24 Maritime aerosol model at 90% relative humidity (M90) to provide the 
aerosol optical properties. The solution of the reciprocal problem was obtained using the successive- 
order-of-scattering method 25 for solving the radiative transfer equation. The computations are 
provided for aerosol optical thickness r a = 0.1 and 0.2, and wavelengths (A) of 443 and 555 nm, 
the principal spectral regions used to estimate the phytoplankton pigment concentration from L w . s 
These values for r a are typical of those in a pure maritime atmosphere. 26-28 

We begin by examining sun- viewing geometry typical of scanning ocean color sensors: viewing 
in the perpendicular plane (0 = 90°) at the center of the scan (9 V zz 0), and near the scan edge 
(9 V zz 45°), where 9 V = 180° — 9 , i.e., 9 V is the angle between the direction the radiometer is viewing 
and the nadir. The computations are performed for 9 0 = 0, 20°, 40°, and 60°, and thus, cover the 
sensor’s full scan as it progresses along the orbit with ever increasing solar zenith angles. 

Figure 3 provides the error A L w /L w = (t — t m )/t m in assuming that L w is uniform. In the 
case of a pure Rayleigh-scattering ocean (Figures 3a and 3b), the maximum error is < 1% for all 
9 V and Furthermore, the error depends only weakly on the wavelength and r a . In contrast, in 
the case of a pure particle- scattering ocean (Figures 3c and 3d), errors as large as 4% are observed, 
along with considerable dependence on wavelength and on r a . The large errors are explained by 
the much larger range of variation in L u (9,<f>) for the particle-dominated ocean compared to the 
Rayleigh-scattering ocean (Figure 2), which increases the affect of the £ u (-<0/^ti(-fo) m ®q. (5). 
The dependence on A and r a is explained by the fact that as r a increases or as A decreases, the sky 
radiance becomes more diffuse, spreading Lr(£) in Eq. (5) over a larger range of directions (£). 
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Similar results axe obtained in other sun-viewing geometries. For example, Figure 4 shows 
AL W /L W as a function of the viewing azimuth angle for situations in which 9 V = 9 0 . The resulting 
errors are similar to those in Figure 3, although in the case of a Rayleigh-scattering ocean, the 
error near <f> v = 0 or 180° (the principal plane) is much greater them that near <f> v = 90° (the 
perpendicular plane). As most ocean color instruments scan closer to the perpendicular plane them 
the principal pleme, this is not considered import emt. Figure 5 provides the error as a function of 
the viewing angle in the perpendicular pleme. The behavior is similar to that in Figure 4, although 
the error does not exceed 4% over the meaningful range of 9 V (0-60°). 

We have adso carried out computations of AL W /L W using actual measurements of radiance 
distributions made by Voss 13 ’ 29 on the R/V New Horizon West of San Diego at 32°40'N and 
121°18 , W. In this case the water itself contributed < 1-2% to the total backscattering in the 
blue. 13 The radiance distributions at 450 and 500 nm are provided in Figure 6 for 9 0 = 60°. Note 
the similarity of these to those for the particle- dominated ocean (Figures 2c and 2d); however, the 
total variation (minimum to maximum) of the measured L u (d, <f>) is smadler, as would be expected 
in the presence of multiple scattering. Figure 7 provides the resulting A L w /L w as a function of 
<f> v with 9 V = So (Figure 7a) and as a function of S v in the perpendicular plane (Figure 7b) for 443 
and 510 nm. For the 443 nm computations, L u at 450 nm was used, while for 510 nm, L u at 500 
nm was used. The resulting errors are similar to those provided in Figures 4 and 5. The strikingly 
different behavior of the error at the two wavelenghts for 9 V = 0 in Figure 7a and for <f> v = 0 and 
60° in Figure 7b results from differences in the two radiance distributions. 

Considering that the goal of atmospheric correction is recovering L w at 443 nm to within ±5% 
in very clear ocean water, we see that the assumption that L u (9 , <f>) is uniform can lead to significant 
error in L w in such situations, i.e., errors similar in magnitude to Lother, when the phase function 
is similar to Petzold’s, or when the radiance distribution is similar to that measured on the New 
Horizon. However, in more productive waters, or in the green, the error in the recovered L w induced 
by this assumption will usually be small compared to that induced by error in removing Lother, 
i.e., the error in L ® itself will be considerably more than the error in L w induced by replacing t by 

r. 
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Finally, most algorithms for relating L w to water constituents involve the use of radiance 
ratios. 5 For algorithms using the ratio I ti; (443)/L w (555), e.g., the CZCS phytoplankton pigment 
algorithm, Figures 3-5 suggest that the error in this ratio induced by replacing t by t m is < 2% with 
the exception of a very clear ocean (Rayleigh scattering) viewing near <j> v = 0 with B v — 0 o = 60°, 
for which the error is ~ 3% (Figure 4b). 

4. Dependence of r on atmospheric parameters 

As it is clear that usually t can be replaced by £*, it is important to understand how this quan- 
tity depends on the optical properties of the atmosphere. Assuming that the water-leaving radiance 
is totally diffuse and spatially uniform, the work of Tanre et al. 14 provides an approximation to t m : 

£*(0) = eip|-(o.48r r + 0.50(1 - <j)r a )/ costfj , (7) 

where r a is the aerosol optical thickness, and g is the cosine of the scattering angle averaged over 
the aerosol scattering phase function. Similarly, Gordon et al. 6 approximated t m by 

t m (6) = exp[-^0.50r r + (1 - u; a F)r a ^ / costfj , (8) 

where u> 0 is the aerosol single scattering albedo, and F is defined by 

where P a {a) is the aerosol scattering phase function (normalized to 47r) for a scattering angle a, 

COS a = fifjl 1 + y/(l - /i 2 )(l - /i' 2 ) cos 

and /i = cos 6. These formulas yield similar results for u; a = 1. Because of the typically large aerosol 
forward scattering, i.e., g — 0.7 and F ^ 0.9, the affect of aerosols on t m (0) is significantly less than 
on the direct transmittance. In the absense of aerosols, both formulas yield £* % exp[— ry/2 cos#]. 
This approximation can be shown to be exact near the single scattering limit, i.e., t*(8) — ► 1 — 
r r /2cos0. The quantity t* was used as an approximation to t m in processing CZCS imagery, 30 i.e., 
r a was set to zero in Eq. (8). 
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We used Eq. (3) to compute t*(0) in the absence of aerosols, t*. This is provided in Figure 8a 
at the centers of severed of the SeaWiFS spectral bands. Figure 8b provides the error in Eq. (8) for 
the spectral band at 412 nm (r r = 0.32). It is seen that Eq. (8) overestimates £* by — 1.5%, which 
is approximately the loss in radiance upon transmission through the air-sea interface at small solar 
zenith angles. 

To investigate the influence of aerosols on t M , we used several of the Shettle and Fenn 24 aerosol 
models. Specifically, we used their Urban model at 50% relative humidity (U50), their Tropospheric 
model at 50% relative humidity (T50), and their Mauitime model at 90% relative humidity (M90). 
M90, T50, and U50 have w a values at 443 nm of 0.9951, 0.9643, and 0.6534, respectively, i.e., 
increasing amounts of absorption. As expected from Eqs. (7) and (8), the variation of £*(0) with 
6 is similar to that observed in Figure 8a. Figure 9 provides the dependence of t * on r a for each 
aerosol model at 0 = 20° amd 60° for 443 and 555 nm. t M is seen to be a strong function of u> a . 

Gordon 10 has shown that in the case of a strongly- absorbing aerosol, the quantity Lother in 
Eq. (1) is a sensitive function of the vertical distribution of the aerosol. This is particularly true 
in the blue region of the spectrum. Thus, it is of interest to understand the influence of aerosol 
vertical structure on t m . Figure 10 compares t m at 443 nm for the U50 model (strongly absorbing) 
computed for a two-layer atmosphere (aerosols at the bottom) and a one-layer atmosphere in which 
the aerosol are uniformly mixed with Rayleigh scatterers. The results clearly demonstrate that t m 
is nearly independent of aierosol verticad structure. 

Examination of Figure 9 shows that the error in replacing t m by is typically small as long 
as r a is not too large and u; a ~ 1, For example, at r a as high as 0.3 the error in t 9 at 6 = 20° is 
only ~ 2 — 3% for the M90 model, which is expected to be readistic over much of the open ocean. 
Furthermore, the error in t* has the same sign at both wavelengths, so the error induced in radiance 
ratios will be even smadler. As r a = 0.3 would be large for such areais, 20 ~ 28 this approximation should 
be excellent for most oceanic areas. In contrast, for U50 the error under the same conditions will be 
aw high as ~ 15%, but the error in radiance ratios will still be small. These resuts suggest that the 
approximation t m = tl will probably be acceptable (considering the expected error in Lother) if the 
algorithms that use L w 's are in the form of radiance ratios, as in the case of CZCS. For adgorithms 
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using L w 's in other ways, the error can be excessive unless the aerosol is included in the t m estimate. 
Fortunately, the computation of t* is simple given the aerosol model. Recall that the Gordon and 
Wang 7 algorithm utilizes that variation of Lother in the near infrared (NIR) to select an aerosol 
model, from a set of candidate models, for estimating Lother in the visible. Once an aerosol model 
is chosen, Lother in the NIR can be used to estimate r 0 at all wavelengths. An aerosol model and 
r a are all that are needed to accurately compute t m . Given the fact that t m depends exponentially 
on r a (Figure 9), i.e., t m (9) ^ A(9) exp[- B(9)r a ] 1 for a given 9 only two parameters are needed to 
compute t*. To effect the removal of Lother , a set of lookup tables (LUTs) have been prepared that 
relate Lother to r a for each candidate aerosol model, sun-viewing geometry, and wavelength, by 
solving the radiative transfer equation (RTE) for a two-layer atmosphere (aerosols on the bottom). 
These existing solutions to the RTE can be used with Eq. (3) to compute the required £*, since 
t m does not depend on the vertical structure of the aerosol (Figure 10). As the Lother LUTs have 
been prepared for 9 0 = 0(2.5°)80°, i.e., 33 values of 9 0 , LUTs relating t m to r a with the same 
resolution require only 66 constants for each wavelength and candidate aerosol model of interest. 
In contrast, the LUTs for Lother require ~ 35,000 constants per aerosol model per wavelength. 
Thus, computation of a precise value of t m (9) appropriate to a given candidate aerosol model is not 
a challenge. 

5. Concluding remarks 

In this paper we have examined the effects of factors influencing the diffuse transmittance 
of the water-leaving radiance to the top of the atmosphere: the angular distribution of upwelling 
radiance beneath the sea surface; and the concentration and optical properties of the aerosol in the 
atmosphere. Several conclusions are possible based on the analysis. 

First, the error in L w made by assuming that L u is uniform is < 4% for viewing geometries 
typical of ocean color observations, and aerosol optical thicknesses typical of the open-ocean marine 
atmosphere. As the error in retrieving L JJ, from Lt will usually exceed 4% by a considerable amount, 
the departure of L u from uniform usually can be ignored. An exception is for low values of C in 
the blue where L ® is large and its removal from Lt can be expected to be in error by < 5%. 
Thus, for atmospheric correction over clear water in the blue, it will be necessary to esitmate the 
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angular distribution of L u in the computation of t in order to insure that the error in L w remains 
< 5%. This can be accomplished using an estimate of C, obtained with in an iterative manner 
as suggested by Morel and Gentili. 31 In contrast, for sensors that are capable of viewing over a large 
range of azimuth angles, e.g., the polarization and directionality of earth reflectance (POLDER) 32 
instrument, the unif orm- L u assumption can lead to errors of more than 6% in L w (Figure 4b) and 
is highly dependent on the viewing azimuth relative to the sun. 

Next, when biophysical products, e.g., C, are derived from ratios of L w at two wavelengths, 
the effect of replacing t by t* is significantly smaller than the effect on L w . 

Finally, in the absence of aerosols, t* can be computed from simple formulas with an, easily 
correctable (Figure 8b), error of < 2%. In the presence of aerosols, t * is a strong function of 
the aerosol absorption (Figure 9) but is independent of the aerosol vertical distribution (Figure 
10). Given that the atmospheric correction algorithm 7 provides candidate aerosol models, accurate 
computation of t* can be made using the models choosen to remove Lother- 
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6. Appendix: The reciprocity principle 


As we were unable to find the result in the literature, we provide here a derivation of the 
reciprocity principle for a medium in which the refractive index m is a function of position. We 
then apply this to derive Eq. (2) in the text. 


In the absence of internal sources, with such a medium occupying a volume V and illuminated 
from the outside, the radiance L(l>£) along a ray is governed by the radiative transfer equation 


(RTE), 


± L{U) 

dl m 2 (/) 





( 9 ) 


where f is a umit vector tangent to the path of the ray, and l is measured along the ray. In this 
equation, (3(l\ £' — ► £) is the volume scattering function for scattering from to £, c(/) is the 
beam attenuation coefficient of the medium, and cfQ(£') is a differential of solid angle around £' . To 
derive the reciprocity principle, we imagine a given medium with two different radiance distributions 
incident on a volume V from the outside. We index the solution to these two problems by the indices 
1 and 2. Thus, for problem 1, we have 


± w 

dl m 2 {l ) 


-c(l) 


£x(M) 

m *(0 



- -MUO 
V m>(l) 


dm'), 


(10) 


with a similar equation for problem 2. Now in problem 2, we reverse the sign of £ everywhere, i.e., 


± I L 2 {i t -j) 

dl m 2 (l) 
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L ,(/, -Q 

m 2 (l) 
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U m 2 (l) 


dm 1 ), 


( 11 ) 


where the minus sign on the left-hand-side is introduced so that the direction of increasing l is in 
the direction +£ in both problems 1 and 2. The integral term in Eq. (11) can be rearranged in the 
following manner: first, since the direction £' is an integration variable, it can be replaced by 
everywhere in the integral; next, we recall that j3(l\ — ( 9 — ► -£) = /}(£;£ — ► £'); and fin silly, we note 
that £') = dfi(f'). Thus, the integral term in (11) can be written 



^ ™ 2 (*) 


dm 1 ). 
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Now, multiply Eq. (10) by L 2 {(, -i)/m 2 {t), and Eq. (11) by Li(l,0/m 2 (i), and subtract. 
Then multiply the result by m 2 {l) dfi({) dV and integrate over all ft(£) and over V. The left-hand- 
side of the result is 


/ dm [ 

Ja* Jv 


r 


ti (<,f) -0 

m 2 (() m 2 (f) 


dV, 


and the right-hand- side is 


Jr" I I «■ - « W ^ ^ 

~Jr "j «(ft. 

Clearly, the terms on the right-hand-side add to zero, and 



m 2 (f) m 2 (f) 


dV = 0 


Now, we write dV = dA(£) dl, where dA(£) is a differential of area with normal in the direction of 
£. Then using the fact that m 2 (t)dA(£) dfl(f) is constant along the path of the ray, 33 


or 



Lxjlj) L 2 {l,-j) 

m 2 (f) m 2 (f) 


m 2 (/)dn(Od^(o 


m 2 (i) 


I *3 


dQ(t)dA(£) 


= 0, 


<U = 0 


where fi is the beginning of the path in V (the position where the ray enters V) and t 2 is the end 
of the path (the position where the ray exits V). If n is the outward normal to the surface of V, 
then dA(£) at t\ or l 2 is £ • hdS , where dS is the associated differential element of area on the 
surface of V, and this equation becomes 



* m 2 (p) 


dm = o, 


( 12 ) 


with p specifying the position of a point on the surface. This is the reciprocity principle for a 
medium with variable m. It can be rearranged to read 



• n<0 


K • "I 


Liipjwjm 

m 2 {p) 


m 2 {p) 
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(13) 
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where the integration over fi(£) is only over directions for which £ • n < 0. When there are internal 
sources in the medium of intensity density Q(r, £), where r specifies the position of a point in the 
medium, a term Q{r,£)/m 2 (r) must be added to the right-hand-side of Eq. (9). In this case, it is 
easy to show that the right-hand-side of Eq. (13) becomes 

m 2 (f) m 2 (r) 



It is straightforward to use this to derive Eq. (2) in the text. Consider a volume with upper 
surface at the top of the atmosphere, and lower surface just beneath the sea surface. Then, for 
problem 1 choose the incident radiance on the top of the atmosphere (TOA) to be that of the solar 
beam, i.e., Li(pr>0 = Fo6(£ - £o)> where pr is a TOA point, and £o is the direction of the solar 
beam. We assume that there is no upward radiance incident at pa, a point just beneath the sea 
surface, i.e., Li(pb,£) = 0 for £ • h < 0. This is referred to in the text as the reciprocal problem. 
For problem 2 (the direct problem in the text), we let Z 2 (pr>0 = 0 for £ • n < 0 (no incident 
radiance on the TOA), and Li{pb^ 0 be specified for £ • n < 0, i.e., a specified upward radiance 
distribution incident on the bottom surface. Then applying Eq. (13), we have 


^ 2 ( PT j -fo) = 


Fo\to 


i— / 

fo • "I Jn d 


I i • «l 




mi 


dn(t), 


where Qj is the full solid angle in the downward direction, and m w is the refractive index of water. 
Letting be the direction of the refracted solar beam in the water for a flat surface, 


h L *id » 1 ^ 
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In the absence of absorption and scattering in the atmosphere, 


dnd). 


(14) 


Li(phJ) = T f (io)F 0 l ^4s(i - fa, 

\Zo • n l 

where T/(£) is the Fresnel transmittance of the interface. Equation (14) then yields 

L2(pt,-£o) — -^y^ L 3 (pb, -£o), 
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which is the familiar relationship for the propagation of radiance across the air-sea interface. If 
we replace L 2 {p~By -£o) m Ec l* ( 14 ) by Lw(Pb, -£o), the water-leaving radiance just above the sea 
surface, i.e., 


m~ 


Li{Pb, -Z'q) - — ~ j- L w (pb, -Co), 

Tf\Zo) 


then, Eq. (14) becomes 


Q- = : — \ [ \i»n\L l {ph,t) L *( p f' &~ dn{£) = 

L w {phy -£o) i'olfo • n\T f m Jn d L 2 (pa, -£o) 


t(—£o) is the quantity we defined as the diffuse transmittance in the text, Eq. (2). 
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Figure Captions 


Figure 1. Comparison of the scattering phase functions for Rayleigh scattering (dotted line) and 
Petzold’s 20 turbid water measurements. 

Figure 2. Subsurface upwelling radiance predicted using the Rayleigh and Petzold phase functions. 
The radiances are normalized to their maximum value. Panel (a) Rayleigh with Oq = 40°; Panel 

(b) Rayleigh with 0o = 60°; Panel (c) Petzold with 0o = 40°; Panel (d) Petzold with 9 o = 60°. 

Figure 3. The error in t , induced by assuming that L u is uniform, as a function of wavelength and 
aerosol optical thickness at the SeaWiFS scan center and scan edge. The M90 aerosol model is 
used in all the computations. Panel (a) Rayleigh at the scan center; Panel (b) Rayleigh at the scan 
edge; Panel (c) Petzold at the scan center; Panel (d) Petzold at the scan edge. 

Figure 4. The error in t, induced by ass uming that L u is unif orm, as a function of wavelength, 
aerosol optical thickness, and viewing azimuth (^„) for 0o = 9 V The M90 aerosol model is used in 
all the computations. Panel (a) Rayleigh with 9q = 40°; Panel (b) Rayleigh with 9q = 60°; Panel 

(c) Petzold with <?o = 40°; Panel (a) Petzold with 0o = 60°. 

Figure 5. The error in t, induced by assuming that L u is uniform, as a function of wavelength, 
aerosol optical thickness, and viewing angle (0„) in the perpendicular plane of the sun. The M90 
aerosol model is used in all the computations. Panel (a) Rayleigh with 9q = 40°; Panel (b) Rayleigh 
with 0 O = 60°; Panel (c) Petzold with 0o = 40°; Panel (a) Petzold with 0o = 60°. 

Figure 6. Subsurface upwelling radiance distributions measured by Voss with the RADS system. 13 ’ 29 
Panel (a) A = 450 nm, 0 O = 58.2°; Panel (b) A = 500 nm, 0 O = 59.7°. 

Figure 7. The error in t, induced by assuming that L u is uniform, as a function of wavelength 
and aerosol optical thickness, using the Voss subsurface upwelling radiance distribution. The M90 
aerosol model is used in all the computations. Panel (a) 0o = 60° viewing in the perpendicular 
plane to the sun; Panel (b) 0„ = 0 O = 60° with viewing azimuth <(> v . 
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Figure 8. Computed t m for a pure Eayleigh-scattering atmosphere. Panel (a) Values computed 
using radiative transfer theory. From bottom to top the curves correspond to A = 412, 443, 490, 
555, 670, and 865 nm. Panel (b), Error in t m resulting from using Eq. (7) in the text. 

Figure 9. Value of f* as a function of aerosol model and optical thickness r a . Upper curves are for 
A = 555 nm , lower curves for A = 443 nm. Panel (a) 9a = 20°; Panel (b) 9a = 60°. 

Figure 10. Infl uence of vertical structure on t m for the U50 aerosol model at 443 nm. Lines 
correspond to placing all of the aerosol in the bottom layer of a two-layer atmosphere. Dots 
correspond to the mixing the aerosol uniformly with gas molecules throughout the atmosphere. 
Curves from top to bottom are for 9a = 20°, 40°, and 60°. 


24 




0 


Figure 1. Comparison of the scattering phase functions for Rayleigh scattering (dotted line) and 
Petzold’s 20 turbid water measurements. 
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Figure 2. Subsurface upwelling radiance predicted using the Rayleigh and Petzold phase functions. 
The radiances are normalized to their maximum value. Panel (a) Rayleigh with 9 0 = 40°; Panel 
(b) Rayleigh with 9 0 = 60°; Panel (c) Petzold with 9 0 = 40°; Panel (d) Petzold with 9 0 = 60°. 
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Figure 3. The error in t, induced by assuming that L u is uniform, as a function of wavelength and 
aerosol optical thickness at the SeaWiFS scan center and scan edge. The M90 aerosol model is 
used in all the computations. Panel (a) Rayleigh at the scan center; Panel (b) Rayleigh at the scan 
edge; Panel (c) Petzold at the scan center; Panel (d) Petzold at the scan edge. 
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Figure 4. The error in t, induced by assuming that L u is uniform, as a function of wavelength, 
aerosol optical thickness, and viewing azimuth (0„) for do = d v The M90 aerosol model is used in 
all the computations. Panel (a) Rayleigh with do = 40°; Panel (b) Rayleigh with do = 60°; Panel 
(c) Petzold with 0<j = 40°; Panel (a) Petzold with d 0 = 60°. 
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Figure 5. The error in t , induced by assuming that L u is uniform, as a function of wavelength, 
aerosol optical thickness, and viewing angle (0„) in the perpendicular plane of the sun. The M90 
aerosol model is used in all the computations. Panel (a) Rayleigh with 0 O = 40°; Panel (b) Rayleigh 
with 9 0 = 60°; Panel (c) Petzold with 6 0 = 40°; Panel (a) Petzold with & 0 = 60°. 









Figure 7. The error in f, induced by assuming that L u is uniform, as a function of wavelength 
and aerosol optical thickness, using the Voss subsurface upwelling radiance distribution. The M90 
aerosol model is used in all the computations. Panel (a) 0$ = 60° viewing in the perpendicular 
plane to the sun; Panel (b) 0 V = 0q = 60° with viewing azimuth <f> v . 
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Figure 8. Computed t m for a pure Rayleigh-scattering atmosphere. Panel (a) Values computed 
using radiative transfer theory. From bottom to top the curves correspond to A = 412, 443, 490, 
555, 670, and 865 nm. Panel (b), Error in t m resulting from using Eq. (8) in the text. 
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Figure 9. Value of t as a function of aerosol model and optical thickness t„. Upper curves are for 
A = 555 nm, lower curves for A = 443 nm. Panel (a) 0 O = 20°; Panel (b) 0 O = 60°. 




Figure 10. Influence of vertical structure on t m for the U50 aerosol model at 443 nm. Lines 
correspond to placing sill of the aerosol in the bottom layer of a two-layer atmosphere. Dots 
correspond to the mixing the aerosol uniformly with gas molecules throughout the atmosphere. 
Curves from top to bottom are for 9 0 = 20°, 40°, and 60°. 
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ABSTRACT 


A measurement system for determining the spectral reflectance of whitecaps in the open ocean is 
described. The upwelling radiance is obtained from a ship by observing a small region of the water 
surface over time using a six channel radiometer (410, 440, 510, 550, 670 and 860 nm) extended 
from the bow of the ship. Downwelling irradiance is simultaneously measured and used to 
provide surface reflectance measurements. The system includes a TV camera mounted beside the 
radiometer which provides a visual reference of surface events. Air/water temperature and wind 
speed/direction are also measured along with global positioning system (GPS) data. Calibration 
procedures and radiometric characterization of the system for operation under different sky 
conditions and solar zenith angles is emphasized so that full advantage is taken of ship time 
whenever whitecap events occur. Examples of the spectral reflectance of different foam types 
(thick, dense, fluffy to thin residual patches) generated by the ship’s bow in coastal regions are 
presented and found to vary spectrally. The presence of submerged bubbles in the foam 
measurement results in a lower reflectance at the longer wavelengths. For wavebands in the visible 
region, the spectral reflectance values tend to equalize with higher reflecting foam from thicker, 
fluffier and denser foam types. 



1. Introduction 


Remote sensing of the ocean by satellite requires accurate atmospheric correction. This is 
accomplished by measuring the upwelling radiance at a wavelength where the ocean is known to 
contribute very little or no radiance (Gordon and Wang, 1994a). This is typically done at wavelengths 
beyond 760 nm (and sometimes at 670 nm) where the absorption properties of water are particularly 
strong. In the presence of whitecaps, whose quantity and fractional coverage of the ocean surface varies 
primarily with wind speed, the ocean surface may no longer provide a negligible background for 
atmospheric correction at these longer wavelengths (Gordon and Wang, 1994b). An understanding of 
the additional or augmented spectral upwelling radiance from whitecaps must be quantified not only at the 
correction wavelength but at other wavelengths particularly in the visible region where spectral 
information provides a means of determining water content. 

Previous authors have used a wavelength independent foam reflectance of 50% - 100% (Payne 
1972; Gordon and Jacobs 1977; Maul and Gordon 1975) in conjunction with estimates of fractional 
coverage to determine the augmentation effect Others have measured the spectral variation of foam in 
the laboratory using clear water and a variety of more turbid water types with various concentrations of 
detritus, sediment, etc. (Whitlock et al. 1982). Koepke (1984) established an efficiency factor that 
accounts for the diminishing reflectance and increase in area of the whitecap as it ages thus providing a 
smaller augmentation effect than that of other studies. 

Estimates of fractional coverage of whitecaps have been earned out previously by photographic 
methods, either from the air (Austin and Moran 1974; Ross and Cardone 1974) or from a stationary 
platform over the ocean surface (Monahan 1971; Toba and Chaen 1973; Koepke 1984; Bortkovskii 
1993). Spectral determination of whitecap and foam reflectance has been determined in the laboratory 
(Whitlock et al. 1982) and from foam generated in the surf zone (Frouin et al. 1996). Whidock’s 
laboratory experiments showed a decrease in reflectance with increasing wavelength beyond ~0.8pm. 
However, Frouin reported a much lower reflectance than Whidock at the longer wavelengths with a 40% 
decrease at 0.87pm, 50% at 1.02pm and 95% at l.65pm reladve to the reflectance at 440 nm. The 
dilference between the laboratory and field measurements is thought to be due to the stronger absorpuon 



properties of water at longer wavelengths acting on light reflected from submerged bubbles forced into 
the water by large waves. 

The quantity of whitecaps covering the surface area of the ocean is detected as an augmentation in 
the upwelling radiance observed by the relatively low resolution of a satellite pixel. To measure the 
spectral influence and magnitude of whitecaps in terms of a radiometric signal, the foam-free water- 
leaving radiance must be differentiated from the radiance contribution of whitecaps. This entails 
measurement of all foam types; from thick fluffy, highly reflecting foam as the wave breaks to barely 
distinguishable thin, residual and fragmented patches of foam as the whitecap decays back into the ocean. 
To measure the spectral reflectance of these foam types, particularly within the life cycle of an individual 
whitecap, requires an instrument of comparably high spatial and temporal resolution. Camera images of 
the ocean surface can provide a varied foam assortment, but determination of the foam-free background 
can be complicated by variation in sky reflectance towards the extreme regions of the image. This is 
particularly true when the imaging system is mounted on a platform and at an angle to cover as large an 
area of the water surface as possible. Use of a high resolution imaging system for spectral measurements 
also requires a time delay between images of different wavebands, and the content of the spectral images 
may no longer be identical. 

Using a high resolution imaging system from an airborne platform provides greater coverage and 
a decrease in sky reflectance variation as the height above the water surface is increased but at the cost of 
surface resolution. High resolution spectral data of whitecaps is more feasible to acquire from a surface 
platform particularly in the open ocean where flying an aircraft long distances and providing 
simultaneous ancillary measurements such as wind speed/direction and air/water temperature may 
become costly. 

In this paper, a whitecap radiometer system for operation on board a ship or other surface 
platform is described. The spectral reflectance of the water surface is measured by a 6-channel non- 
imaging radiometer with a narrow field of view and wavebands at 410, 440, 510, 550, 670 and 860 nm. 
The radiometer is deployed from the bow of a ship along with a TV camera to provide a visual reference 
of the radiometrically measured surface. In conjunction with the radiometer, a deck cell (cosine 



collector) monitors downwelling irradiance with matching wavebands thus providing surface reflectance 
measurements. By making many high spatial resolution measurements of the water surface at a fixed 
periodic rate, spectral data of individual whitecaps and their associated foam types (defined by their 
reflectance magnitude) can be analyzed and also provide an estimate of their frequency. Measurement of 
the frequency of different foam types provides a value of fractional coverage and the augmented spectral 
reflectance contribution of whitecaps to the water-leaving radiance as observed by ocean color satellites. 
In addition, the whitecap radiometer system acquires wind speed/direction and air/water temperature as 
well as global positioning system (GPS) information to provide universal time and location. Estimates of 
fractional coverage and augmented reflectance contribution can then be correlated to these parameters. 

In this paper we describe the whitecap instrument system and give particular attention to 
radiometric calibration for field use under various sky conditions and solar zenith angles in order to 
correct for slight deviations in solid angle response of the deck cell’s cosine collector and the effects of 
spectral filter mismatch between channels in the deck cell and radiometer. Calibration for different sky 
conditions provides a robust system that takes full advantage of ship time and provides confident 
reflectance measurements of whitecaps whenever they occur. In addition, performance of the system is 
discussed and foam data generated by the ship’s bow in coastal waters is examined. 

2. Whitecap Radiometer Instrument Description 

Our strategy has been to measure whitecap reflectance from a surface platform such as a ship; 
either stationary or underway. The whitecap measurement system consists of a 6-channel radiometer 
aimed at the water surface, a downwelling irradiance cosine collector with matching wavebands, water 
temperature and air temperature sensors, an anemometer to determine wind speed and direction, and a 
global positioning system (GPS) to record location, GMT time and determine ground speed and course 
heading. Taking many reflectance measurements of a small area of the ocean surface over time is 
equivalent to capturing a large field of view with high spatial resolution at one time. 

The 6-channel radiometer has wavebands at 410, 440, 510, 550, 670 and 860 nm and is aimed 
at the water surface with each channel covering ~ 1° field of view. A TV camera (Sony SSC-C350 color 



CCD camera) with a greater field of view (-40° - 50° using a 6 mm focal length video TV lens) is 
mounted beside the radiometer to provide a visual reference of radiometer measurements and is recorded 
on video tape using a video cassette recorder (Sony EV-C100). A triggerable analog switch has been 
implemented to allow individual TV frames to be date and time marked whenever the radiometer acquires 
data, thus providing a frame by frame reference. 

The radiometer and TV camera assembly are extended over the bow of the ship as a unit by 
means of a boom. The boom is comprised of a number of segments of equal length and any number can 
be deployed to provide extension from -1 m to -10 m depending on the sea state and at the same time 
ensure that optical shadowing and wind sheltering ship effects are minimized. 

To provide reflectance measurements of the ocean surface, downwelling irradiance is measured 
by the cosine collector (deck cell) at 6 matching wavebands, positioned on the ship to avoid shadowing. 
Wind speed and direction are measured by an anemometer placed in an unobstructed position towards the 
bow of the ship. Water temperature is measured with an accuracy of 0. 1°C close to the surface by means 
of a small thermistor with signal conditioning circuitry to provide a linear and stable output The 
thermistor’s signal and power cable is fed down the center of a polypropylene rope with the thermistor 
exposed roughly 10-20 cm from the end of the rope. The buoyancy of the polypropylene rope keeps the 
thermistor at or near the surface and the remaining few centimeters of the rope are flayed to ensure the 
thermistor does not skip out of the water as the ship moves. Air temperature is measured with an 
identical thermistor housed in a white plastic louvered cylinder to prevent direct solar heating. It is 
suspended from the bow of the ship close to the water surface but at a sufficient height to prevent contact 
with sea spray. 

The analog signals from the radiometer, deck cell, anemometer, GPS, water and air 
temperature sensors are sampled at a fixed periodic rate. Upwelling and downwelling radiometric data 
along with air/water temperature, wind speed/direction are sampled -7 times a second condnuously for 
about 30 seconds (providing 200 contiguous samples), after which GPS data; GMT time and location, 
are appended to file. The acquisition sequence is repeated until a time determined by the operator. With 
the exception of the GPS data, all parameters including radiometer and deck cell data are measured 



simultaneously (within 20 jis). Making many measurements of the reflectance from small areas of the 

ocean surface as the ship moves through the water allows one to estimate white water coverage as a 
function of wind speed and air/water temperature, as well as providing spectral information. 

All cables from the various measuring components arrive at the electronics console which is 
essentially an interface for individual cable connections as well as providing power and receiving analog 
signals. The analog signals are processed through differential amplifiers, to eliminate ground fluctuations 
that may originate over the long lengths of cable, before being digitized. The pre-conditioned analog 
signals are fed to a National Instruments NB-MIO-16XL acquisition board resident in a Macintosh 
Power PC. This acquisition board is configured to accommodate 16 referenced single-ended inputs and 
provide 16 bits A/D resolution. The acquisition logic, processing sequence and storage to file is 
controlled by Lab VIEW software. The software also provides graphing and numeric display capabilities 
of the data as it is acquired. Data acquisition is allowed to free run and is halted at a time determined by 
the operator. Data files of various sizes are subsequently read into a statistical software package (Matlab) 
to produce detector background corrected and radiometrically calibrated reflectance values of the ocean 
surface along with air/water temperature and wind speed profiles. 

a. The Radiometer 

Each channel of the radiometer consists of an optical tube aligned parallel to one another. Light 
entering each channel passes through a BK-7 window with a broadband AR coating. Behind this 
window, light for each channel is collected by a 25 mm diameter, 100 mm focal length achromat with a 
broadband AR coating. The converging light passes through a 10 nm bandpass spectral interference 
filter. Behind the interference filters, broadband blocking filters are included to enhance out-of-band 
blocking of unwanted light The field of view for all channels (~1° ) is set by a 1.8 mm pinhole aperture 
placed at the focal point of the lens with the active area of the detector positioned 1.5 mm behind the focal 
plane. The detector for each channel is a silicon photodiode (Hamamatsu S 1406-04) packaged with an 
integrated operational amplifier. The voltage output is fed into a low noise precision opamp (Analog 
Devices PM- 1012GP) with a gain appropriate to each channel to offset variation in photodiode spectral 



response before being sent to the electronics console. Power (+/-12 VDC) sent to the radiometer is 
DC/DC converted and further regulated to supply +/-5 VDC to the photodetectors along with +/- 12 VDC 
for the low noise precision opamps. Five of the six channels are arranged in a concentric fashion about 
the long axis of the housing with the 410 nm channel taking up the central position. The complete 
radiometer assembly is housed in a black anodized aluminum cylinder 30 cm in length and 1 1 cm 
diameter (See Fig. 1). These dimensions include an extension hood (5 cm in length) that is provided to 
prevent rain falling on the window and to reduce light that may be reflected and refracted by sea spray 
and salt deposits that inevitably collect on the window. 

b. The Deck Cell 

The deck cell is similar to the radiometer with the exception of the collection optics and its 
housing. The deck cell utilizes a cosine collector with a ~5 cm diameter white diffusive acrylic plastic 
disk (Rohm & Hass, no. 2447) protruding ~ 1 cm beyond a dark rigid PVC base. Moving radially 
outward, the base is terraced, increasing in height with radius to a diameter of 124 mm at which point the 
height is flush with the top of the acrylic disk. The protruding acrylic disk provides a cosine response 

and in conjunction with the terraced base defines a 2 tc field of view as well as ensuring a cosine response 
at the larger zenith angles. 

Five of the photodiodes with interference and blocking filters are arranged in a circle and slightly 
angled to view the center of the acrylic disk from behind, while the sixth photodiode-filter assemblage is 
centered behind the acrylic disk. In this way the cosine response is common to all channels of the deck 
cell. The photodiodes, electronics, and optical filters are identical to the 6-channels in the radiometer but 
with a different gain on the precision opamps and with different optical apertures to accommodate the 
spectral variation of the downwelling light and detector response. 

The 15 cm diameter x 13.5 cm long cylindrical housing is made from rigid PVC plastic mounted 
in a 2-axis gimbal. The bottom base plate of the housing is made from lead and in combination with the 
gimbal mount helps the cosine collector maintain a horizontal attitude despite any low frequency small 
amplitude rocking motion of the test platform. 



3. Radiometric Calibration 

a. Reflectance calibration 

Laboratory radiometric reflectance calibration of the deck cell and radiometer was carried out 
using a calibrated 1000 Watt (FEL) quartz-halogen lamp. The deck cell was positioned 50 cm from the 
lamp with the surface of its cosine collector perpendicular to normal incidence. For the radiometer, a 
Spectralon reflectance plaque was normally illuminated by the lamp 50 cm away and viewed by the 
radiometer at a 45° angle. The reflectance of the plaque is 94.6 +/- 0.2% from 400 to 870 run normally 
illuminated and viewed at an 45° angle. Background dark current readings were taken before and after 
each measurement Absolute calibration was done in this manner before and after deployment in the field 
to verify radiometric stability. 

b. Radiometric performance 

To assess the radiometric performance of the radiometer and deck cell in determining reflectance 
under the variety of natural sky conditions and solar zenith angles occurring during a ship cruise, the 
radiometer and deck cell were tested outside the lab. The radiometer was aimed vertically downward at a 
Spectralon reflectance plaque with the deck cell positioned beside the plaque. The plaque surface was 
level and at the same height as the flat surface of the deck cell cosine collector. The radiometer was 
positioned to avoid direct shadowing of the plaque and the deck cell, and so that the radiometer blocked 
out only a small section of the sky. 

Under overcast conditions with dark and bright patches of cloud passing overhead, variations in 
the measurement of the reflectance value of the plaque were observed. A set of reflectance data were 
taken around noon with sky conditions going from one extreme (very overcast - sky totally covered and 
no distinct shadows observed) to another (direct sunlight breaking through a large clear patch with solar 

elevation angle -60°). The variation (maximum to minimum) in reflectance values observed during these 

conditions for each channel (410 through 860 nm) were: 0.040, 0.045, 0.039, 0.050, 0.058 and 0.063 
with the mean reflectance value for each found to be 0.987, 1.022, 0.969, 0.993, 1.001 and 1. 122 
respectively (See Fig. 5), the reflectance variation of the reflectance plaque and correction will be 



discussed later. Calibration factors determined under laboratory conditions obviously were not sufficient 
to measure real field conditions with a changing angular distribution of the 2k light field. These results 

led to a rigorous calibration and understanding of the radiometric system so that we could be confident ot 
reflectance data acquired under many sky conditions and times of day in field measurements. 

c. Linearity 

The first test was to check the linearity of both the deck cell and radiometer in order to identify 
and eliminate any gross malfunctioning of the system. The photodiodes are reversed biased and as such 
should be highly linear. However, a damaged component, loose wire or bad connection could exhibit 
non-linearity at higher light levels that may not be seen at lower light levels of the laboratory calibration. 
As the gain and throughput of these instruments is specifically set for outdoor conditions, taking them 
into the lab where the proportion of blue light is comparatively low makes it difficult to cover the full 
dynamic range that might be encountered outdoors particularly for the shorter wavelength channels. 

The linearity test was carried out by moving a constant light source specific distances from the 
instrument in question and the output signal recorded. In the case of the deck cell, the cosine collector 
was removed and a diffusive surface of much greater transmissivity put in its place in order to achieve an 
equivalent illumination intensity found outdoors on a clear day with the sun overhead. The diffusive 
surface was required to reduce non-uniform illumination effects observed by the off-axis channels when 
the lamp source (140W) was brought close to the deck cell. 

The radiometer, with a narrow field of view, was tested in much the same way but required a 
1000W lamp source and was tested intact without any alterations. A diffusive surface with low 
transmission was placed a fixed distance in front of the radiometer to provide uniform illumination as the 
lamp source moved position. Both the deck cell and radiometer were found to be highly linear. 

d. Angular response of deck cell and plaque 

From the experiment outdoors, it was noted the that measured reflectance value of the plaque 
would sometimes increase and at other times decrease even though the downwelling intensity remained 



much the same. This suggested a possible sensitivity of the instrumentation to the angular distribution of 
the downwelling light field. A crude test done under a clear sky with the sun high overhead also seemed 
to support this view. An opaque surface sufficiently large to continually block out the direct sun incident 
on both the deck cell and the reflectance plaque was raised and lowered to simulate a changing angular 
light field. By doing this, the reflectance value of the plaque was also observed to change in a consistent 
manner with changing angular light field. 

The angular response of the deck cell and the reflectance plaque was measured in the lab. With 
the lamp source (point) placed ~1.5 meters from the surface of the cosine collector of the deck cell, the 

deck cell was rotated in 10° increments from normal incidence to plus and minus 90°. The reflectance 
plaque could not be viewed from the normal due to the restrictive geometry of placing the detector in the 
same position as the source. Starting with the closest reasonably achievable viewing angle of 2°, 

subsequent measurements were taken at 5° through 85° at 5° increments. Spline curves were then fitted 

to these data points and the results of the angular responses for both the deck cell and the reflectance 
plaque are shown for the 410 nm channel in Fig. 2(a). Also included for reference is the angular 

response of an ideal cosine collector or Lambertian reflector. In Fig. 2(b) the angular response from 0° 

to 85° for both the deck cell and reflectance plaque is shown for all channels in terms of deviation from 
true cosine. 

In most channels, the response of the deck cell and reflectance plaque are similar with slight 
differences becoming apparent at larger angles of incidence. Although both deviate from true cosine, 
they appear to be close in response to each other with the exception of the poorer deck cell performance at 
860 nm. For channels 410, 440, 5 10 and 550 nm the deck cell appears to have an angular response that 
is closer to true cosine than the reflectance plaque particularly at the larger angles. The 670 nm channel 
of the deck cell and reflectance plaque are the most similar yet both deviate from true cosine. In the 860 
nm channel the deck cell response is much lower than the reflectance plaque response. Also, the 
reflectance plaque appears to have a consistent angular response in all channels. 



To better understand the implications of these angular response measurements, performance of 
the deck cell and reflectance plaque operating outdoors under 2 n illumination must be addressed. To 
achieve this, the angular response variation in solid angle with zenith angle must be taken into account. 

In Fig. 3 the angular response has been weighted by 2nsind dd which accounts for the greater solid 

angle viewed with increasing zenith angle (0). Figure 3 shows the solid angle response for the 410 nm 
channel of the deck cell and reflectance plaque and also for a true cosine solid angle response. These 
values are normalized to the ideal cosine sensitivity at 0°. The solid angle response for both the deck cell 

and reflectance plaque in terms of deviation from true cosine is identical to Fig. 2(b). 

For the deck cell and the radiometer to measure correct reflectance of the calibrated reflectance 
plaque outdoors, both the deck cell and reflectance plaque must be corrected for solid angle response to 
the particular downwelling angular light field distribution. To measure true (correct) reflectance, R, in 
each channel, 

r 'T’ 

where L, is true upwelling radiance from a perfect Lambertian reflector, E t is true downwelling irradiance 
measured by an ideal cosine collector, and 

At = L m d r rad ’ 

where L rad is the radiance measured by the radiometer viewing the reflectance plaque and r nd is the 
correction factor, due to deviation of plaque response from true Lambertian, applied to L rad to yield L e 
Similarly, 
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where E * is the downwelling irradiance measured by the deck cell and r* is its correction factor. Also, 
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where TC is the true cosine response, RP the reflectance plaque response, DC the deck cell response, L 
is the downwelling radiance or sky condition and d is the zenith angle. TCR, RPR and DCR are the true 
cosine, reflectance plaque and deck cell responses respectively for a particular sky condition. Therefore, 

R = — ‘sL.Q , 


where C T = — — , and the response correction factor, C r , is dependent on sky condition L(0). For 
RFR 

example, to simulate the overcast sky condition the coroidial expression, L(6) = 1 + 2cos6 (where 9 is 

zenith angle), is used. C r in this case is the product of the angular response of the deck cell and the 
overcast sky distribution, integrated over all zenith angles and divided by the integrated product of the 

angular response of the reflectance plaque and the overcast sky profile over all zenith angles ( 9 = 0° to 
90°). This process is repeated for each waveband channel. See Table 1. 



Other sky conditions have been simulated such as different solar zenith angles on a clear day. 
Data from principal plane measurements taken during clear sky conditions using a sun-tracking 
photometer (Welton, 1996), developed by CIMEL, were used to generate the radiance distribution for 

different solar zenith angles at 0°, 10°, 20°, 30°, 40°, 50°, 60°, 70° and 80°. Radiance measurements 
from the sun photometer-sky radiometer at two different wavebands (440 and 860 nm) passing through 
the sun at a 55° solar zenith angle were used. There was no discernible difference in the correction 

factors, C r , when the two wavelengths were applied, and changing the peak solar radiance by a factor of 
10 only effects the resulting correction factors shown in Table 1 by -0.1%. 

As can be seen from Table 1, applying the laboratory calibration relationship between the deck 
cell and the reflectance plaque to determine reflectance under a clear sky with the sun directly overhead 

(solar zenith angle = 0°, i.e. ‘sun O’) is not valid. 
e. Filter response 

In addition to characterizing the effect of the angular response mismatch of deck cell and 
reflectance plaque under different sky conditions, the spectral filter response must also be examined. A 
monochromator using a 1000W lamp source was used to determine the bandpass characteristics of each 
channel in both the deck cell and the radiometer. The monochromator output bandwidth was set to 1 nm 
to profile the nominal 10 nm bandpass of each filter assemblage. Measurement of the bandpass filters in 
each channel of the radiometer was a straightforward process. However, for the deck cell the amount of 
light passing through the diffusive acrylic material of the cosine collector was too small to afford any 
appreciable signal and was removed as before. 

In order to normally illuminate the filter assemblages in the deck cell, the deck cell housing had to 
be angled with respect to the monochromatic output beam since the filter assemblages and their 
photodiodes (with the exception of the center photodiode-filter assembly) are purposely angled to view 
the center of the cosine collector. This action subsequently raised the question of the viewing angle of 

the interference filters inside the deck cell. The back of the cosine collector subtends a half angle of -28° 



to the photodiode detectors and the field of view for each channel was established by the different sized 
optical apertures placed in front of the filters. They were intended originally to increase light throughput 
at the shorter wavelengths and reduce the throughput at the longer wavelengths. The half angle field of 

view for each deck cell channel was 24°, 23.5°, 18.5°, 18.4°, 16.6° and 16.6° for channels 1 through 6 

(410, 440, 510, 550, 670, and 860 nm). 

From these measurements it became obvious that the interference filters would possess different 
bandpass and transmission characteristics than measured under normal illumination conditions. As the 
actual bandpass profile of the deck cell channel could not be measured directly, due to insufficient light 
output from the monochromator with 1 nm resolution, the bandpasses were calculated. This involved 
accounting for the center wavelength shift to shorter wavelengths, the associated transmission drop off 
and bandpass broadening (Handbook of Optics, 1995). The center wavelength shift, transmission drop 
and bandpass broadening were weighted with the effective increase in light intensity contribution due to 
the increase in field of view of the back side of the uniformly illuminated cosine collector at larger angles. 
The photodiode spectral response was also considered. 

The calibration procedure for the radiometer (viewing the plaque) and the deck cell would account 
for this difference in spectral response only if the spectral profile of the lamp source in the lab was the 
same as the solar spectrum. Obviously it is not, and differences in the spectral profile of the lamp and 
solar spectrum must be accounted for with the broader bandpass of the deck cell channels. In Fig. 4, the 
radiometer and deck cell bandpass profiles (calculated) for each channel are shown along with the 
normalized lamp source and solar spectral profiles. The difference in the estimate of reflectance will be 
affected by going from a laboratory calibration under the lamp source to operation outside under the 
different spectral weighting of the solar spectrum. In Table 2, the bandpass multiplication correction 
factors which have been calculated to offset the error in making reflectance measurements out of doors 

are provided. The bandpass correction factors, C bp (AX), have been arrived at by calculating the ratio of 

the integrated product of the radiometer bandpass and the solar spectrum, to the integrated product of the 
radiometer bandpass and lamp spectrum, divided by the ratio of the integrated product of the deck cell 



bandpass and the solar spectrum, to the integrated product of the deck cell bandpass and lamp spectrum 


or: 


k+Ak 

\R(X)L^(X)d . A 
C = — 

^ BP A + M 




\dc(X)l m „(x) dX 

k 

k + Ak ’ 

jDC(X) L,^(X) dX 

k 


where R is the bandpass response for a particular channel of the radiometer, DC the bandpass response 
for the same channel in the deck cell and L f| _ r and L„,„ are the illumination profiles over Ak for the 
lamp source and sun respectively. 

/ Application of the correction calibration 

In Fig. 5 the overcast sky and 30° solar zenith angle (sun30) correction factors are applied to the 
reflectance data measurements taken of the plaque outdoors (as conditions varied from being very 
overcast to direct sunlight breaking through with 30° solar zenith angle and then eventually to broken 

conditions). The filter bandwidth correction has also been included. Due to the higher gain of the 
radiometer which is specifically set to measure the lower upwelling water reflectance signal and foam 
contributions, the 860 nm channel observing the -95% reflectance plaque saturates when the sun appears 
from behind the clouds. This saturation has not been included in characterizing the 860 nm reflectance 
variation with sky condition. 

Variation in the measured reflectance persists as the real sky condition departs from the modeled 
sky definition. The most noticeable benefit of these correction factors is in the 860 nm channel. 
Regardless of the sky correction, application of the solid angle response and filter bandwidth correction 
results in a mean measured reflectance (over the time series shown) of 0.9944 for overcast or 1.017 for 



direct sun with 30° zenith angle rather than the uncorrected mean reflectance of 1. 122. [It should be 

noted that a reflectance plaque value of 1 has been assumed for simplicity in this analysis but its real 
values of 95.7%, 94.2%, 94.4%, 94.3%, 94.4% and 94.9% for channels 410 through 860 nm are 
applied to the field data]. The best description of an overcast sky and direct sun within the time series are 
shown on the figure. The application of the appropriate sky correction to the data taken under a valid sky 
description depicts a closer approximation to a reflectance value of 1 in all channels. The standard 
deviation from a reflectance of 1 for all channels is 0.0025 for the overcast condition and 0.0028 for the 
direct sun condition. Other regions of the time series are the result of sky conditions that are more 
difficult to describe and therefore to further correct. The standard deviation of variation in reflectance 
over the extreme sky conditions that occurred during the complete time series in Fig. 5 was 0.0081, 
0.0089, 0.0086, 0.0101, 0.0127 and 0.0313 for channels 410 through 860 nm respectively. 

In Table 3, the amount of correction that has been applied to the measured reflectance for each 
channel and for both sky conditions (‘sun 30’ and ‘overcast’) is shown. The resulting accuracy or 
percentage deviation from a reflectance value of 1 is given in the last two columns. For example, the 860 
nm channel data is corrected by 12.77% in the overcast data to read a corrected reflectance of 0.995 (i.e. 
0.005% accuracy). The direct sun data is saturated for this channel. 

g. Correction of field data 

To make reflectance measurements of the water surface, only the deck cell requires solid angle 
response correction for different sky conditions. We have: 



where R is the reflectance we wish to determine, L ws is the upwelling radiance from the water surface 
(and including contributions from below the surface) received by the radiometer and E is the true 
downwelling irradiance. As before. 



where £* is the downwelling irradiance measured by the deck cell and is the sky correction factor that 
must be applied to the deck cell to have true cosine response, i.e.. 
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Table 4 shows the field correction factors , C f , the measured reflectance must be multiplied by in order to 
read correct reflectance. 


4. Experiment Description, Results and Discussion 

The initial measurements with the whitecap radiometer system were made to evaluate its 
performance and ability to provide a dependable data base of air/water temperature, wind speed/direction, 
GPS information as well as radiometric and visual events recorded and date/time stamped onto video 
tape. Initially the system was deployed on the RV McGaw for a 14-day cruise off the coast of Southern 
California. The ship visited a number of stations located in a small area on a daily basis to collect data for 
other oceanographic studies. Weather conditions were far from ideal for whitecap formation, being 
primarily calm with mainly diffuse overcast days. However, data was taken of foam generated by the 
ship’s bow as it moved from station to station. Also, application of the overcast correction to foam data 
acquired during overcast conditions is less sensitive to error than applying a clear sky correction that only 
approximates the true solar zenith angle and also removes the possibility of specular sun glitter in the 
measurements. The radiometer was mounted at different positions over the side of the ship providing a 



data base for radiometric performance of the system and analysis of different stages of foam 
development. 

The control software (Lab VIEW) had originally been set to acquire all data at the same time. In 
this way, radiometric, air/water temperature, wind speed/direction could easily be addressed by UTC 
time or location during later analysis. However, the price for waiting for updated GPS data at every 
sample resulted in a system acquisition rate of about 0.5-1 Hz despite setting the National Instruments 
acquisition board to the optimal acquisition rate. In the rare moments when a whitecap did happen to 
form, the slow update of the system failed to provide enough samples of the event and in many cases 
missed the event altogether. Had the wind speeds being much stronger, yielding a good whitecap 
coverage, the system would have provided sufficient data. 

Another factor to consider in capturing transient white water events was the integration time. The 
0.5-1 Hz sampling rate of the system was not radically affected by varying the integration time due the 
wait period for GPS data. Although the data acquisition rate of the board was optimized (55 kHz 
sequential sampling frequency over 16 input channels), a set number of samples were binned at high 
speed (100 samples) to increase the signal-to-noise ratio, and the average value for each channel written 
to file. This yielded an effective integration time of the radiometer of just under 30 ms, i.e. just less than 
a TV video frame period so visually recorded events could be correlated to radiometric measurements 
with optimal signal-to-noise binning. 

Figure 6 is an example of the bow foam data obtained. In this example, the radiometer samples 
as the ship moves to its next station and half way through this time series the ship slows to complete 
stop. When the ship has slowed to a speed where bow foam is no longer generated in the field of view 
of the radiometer, the water-leaving radiance and sky reflectance can be determined. The sky condition 
was completely overcast and slightly foggy, with a very calm water surface. This measurement of the 
foam-free water could then be subtracted from the foam reflectance to provide the spectral augmented 
foam reflectance contribution although the real meaning of augmented reflectance in this case is lost due 
to the origin of the foam. 



In Fig. 7 the spectral reflectance of foam is shown is terms of different reflectance intervals 
relating to different foam types, i.e. thick fluffy fresh dense foam to residual thin patches, streaky foam 
and submerged bubbles. For this differentiation, any reflectance spectrum with a particular value at 860 
nm that falls within a reflectance interval is averaged. In this example, the set reflectance intervals are 0- 
0.04, 0.04-0.08, 0.08-0.12, 0.12-0.16, 0.16-0.20, 0.20-24, 0.24-0.28, 0.28-0.32, 0.32-0.36, 0.36- 
0.40, 0.40-0.44, 0.44-0.48, 0.48-0.52, 0.52-0.56, 0.56-0.60, and 0.60 - 1.00. The number of spectra 
averaged in any one interval is also shown to give an indication of the measurement accuracy in 
determining the spectral profile for that reflectance level. 

Foam data was also acquired under clear and sunny conditions where the effects of specular sun 
glitter were included. As the ship was there to service test station areas in a small grid, there was 
frequent changes in ship heading, making a fixed radiometric viewing angle with respect to the sun 
difficult to maintain. These changes in ship heading at random times of the day made the quantity of 
specular sun glitter more difficult to estimate. Some specular sun glitter was tolerable but depending on 
the illumination and viewing geometry it could at times radically affect the reflectance signal. In addition, 
it was difficult to avoid ship shadow particularly at lower solar zenith angles. The effects of specular sun 
glitter in radiometrically discriminating white water from nonwhite water events is addressed elsewhere 
(Moore et al. 1997 ). However, much of the data from this cruise was acquired under overcast sky 
conditions. 

In Fig. 8 a collection of foam data taken at different times throughout a very overcast morning, 
starting with local times of 7:45 am, 10:30 am, 1 1:00 am and 12:00 noon (19 October 1995), was 
complied into a single time series. The foam reflectance at the beginning of the time series is highest due 
to thicker fluffier foam layers generated by a higher ship speed. This first part of the time series is the 
same data used in Figs. 6 and 7, followed by starts and stops and slower ship speeds producing a varied 
set of foam reflectances. As before, reflectance intervals relating to different foam types (or foam 
thicknesses) and the foam-free water-leaving radiance and sky reflectance (referred to as the foam-free 
water reflectance) are shown in Fig. 9 along with the number of spectral reflectance values used in the 
estimation of each interval. 



The pronounced drop off at the 670 and 860 nm channels is thought to be due to the strong 
absorbing properties of water at those wavelengths particularly on light reflected from submerged 
bubbles. Also of interest is the drop off in the 440 and more so in the 410 nm channels. These 
suppressed reflectance values at either end of the spectrum become more acute with higher reflectance 
values corresponding, in general, to thicker and more developed foam. 

In Fig. 10 the spectral reflectance of intervals 2, 4, 7, 10 and 12 from Fig. 9 are shown. The 
spectral profiles are normalized at 440 nm along with the background foam-free water reflectance. The 
spectral variation of different foam types can be seen by normalizing their reflectances. The foam-free 
water spectral reflectance profile is different showing a relatively small difference between 410 and 440 
nm. The foam-free water has reflectance values of 0.0431, 0.0466, 0.0509, 0.0535, 0.0370 and 0.0310 
at 410, 440, 510, 550, 670 and 860 nm. The influence of the background foam-free water reflectance 
can be seen particularly in the lower reflectance intervals. The 550 nm value in interval 2 is slightly 
higher than the 510 nm value and the spectral shape appears to be somewhere between the foam-free 
water reflectance and that of a foam reflectance from a higher reflectance interval. While the slope 
between the 670 and 860 nm values for reflectance interval 2 is smaller than for the higher reflectance 
intervals, it is not quite as small as the foam-free water reflectance. The high 510 and 550 nm values in 
both the foam-free reflectance interval and interval 2 seem to eventually disappear in the higher 
reflectance intervals. Since the lower reflectance interval is comprised of the thinnest foam types, with 
patches of single-layered bubbles, measurements of the foam-free water reflectance and water-leaving 
radiance will inevitably be included. With thicker and denser foam types comprising the higher 
reflectance intervals, the measurement of foam-free contributions decreases with higher reflectance 
interval. The probability of including the reflectance contributions from submerged bubbles also 
decreases with higher reflectance intervals. On the other hand, the 410 nm values appear to remain the 
same regardless of reflectance interval. 

In Fig. 1 1 the augmented reflectances for the same foam types used in Fig. 10 normalized at 440 
nm are shown. The measure of augmented reflectance in this case simply means that the foam-free water 
reflectance has been subtracted from the foam reflectance in order to investigate spectral variation with 



foam type only. Many features are similar to the reflectance profiles for these intervals shown in Fig. 10. 
As expected, with the removal of the foam-free water reflectance, the smallest reflectance interval no 
longer shows the similarity to the background foam-free water reflectance as it did in Fig. 10. All foam 
types are essentially identical with the exception of the 670 and 860 nm values which increase with 
higher reflectance interval. Once again the all foam spectra have a lower 410 nm value than the foam-free 
reflectance, and the 440, 5 10 and 550 nm values are almost equivalent, explaining the white appearance 
of foam. However, the values at 410 nm now appear to increase, if only very slightly, with increasing 
reflectance interval. Also the drop off from 440 nm to 860 nm, in these normalized spectra, can now 
easily be seen to decrease in a progressive manner with increasing reflectance interval and therefore of 
foam type; for reflectance interval 2 (with spectral augmented reflectance at 440 nm of 0.0597 taken 
from 357 samples with background foam-free water spectra subtracted) there is a -62% drop, for interval 
4 (spectral augmented reflectance at 440 nm of 0.252 from 94 samples) it is 57%, for interval 7 (spectral 
augmented reflectance at 440 nm of 0.488 from 1 13 samples) it is -53%, interval 10 (spectral augmented 
reflectance at 440 nm of 0.620 from 48 samples) -45% and -41% at interval 12 (spectral augmented 
reflectance at 440 nm of 0.701 from a single sample). 

From this data there appears to be a similar drop off trend in the blue between the 440 nm and 
410 nm values. Once again the greatest drop off occurs for the lowest reflectance interval although not as 
large as between 440 and 860 nm. For interval 2 there is a -20% drop off compared to the foam-free 
reflectance drop off of -7%. At interval 4 a -16% drop, interval 7 a -13.5% drop and at intervals 10 and 
12 a -15 and 17% drop respectively. This increase in drop off in intervals 10 and 12 is possibly due to 
statistical lack of data. 

While the very strong absorption properties of water may explain the lower reflectance values at 
860 and 670 nm, the drop off at 410 nm, on the other hand, raises the question of whether the 
absorption properties of the water in the blue wavelengths are responsible for this effect. The absorption 
spectrum of pure sea water is low from about 550 nm down past 400 nm and becomes less absorbing 
with shorter wavelength. However, the presence of dissolved organic matter, yellow substance, detritus 
and phytoplankton, one or other of which is typically abundant in coastal waters, increases the 



absorption qualities of the water body in the shorter wavelengths. (The spectral absorption coefficient of 
this water was 0.1 1, 0. 10, 0.065, 0.08, 0.475 and 4.3 m' 1 at the corresponding waveleng ths (Roesler 
1995)). In the case of absorption in the red wavelengths, it is the very strong absorption q uali ties of the 
water itself that cause the observed attenuation, and can easily outweigh the effects of any additional 
absorption due to particulates present in the water with the exception of high sediment amounts in turbid 
waters (Whitlock et al., 1982). Even so, the acute drop off at 410 nm cannot be fully explained by water 
absorption when considering the small absorption increase between 410 nm and 440 nm compared to the 
much greater absorption difference between 440 nm and 670 nm for similar values of augmented 
reflectance. 

5. Conclusion 

The whitecap system for determining the spectral reflectance of whitecaps in the open ocean has 
been described and calibrated for optimum use on board a ship or other over water surface platform. The 
calibration for different sky conditions and solar zenith angles, if not exactly represented by the real sky 
conditions at the time of acquisition, shows the degree of error that can be incurred without correction 
and indicates the degree of spectral variation that must be due to real effects in the field. In this paper we 
have also shown how the system works with an understanding of the possible errors that could have 
been incurred in going from laboratory calibration with a lamp source acting as a point source to outdoor 
2jc illumination under the solar spectrum. 

The system has been designed to acquire spectral data of real whitecaps and their frequency rather 
than foam generated from the ship’s bow. By taking many reflectance measurements of the water 
surface over time, a database of the spectral reflectance of different foam types and their frequency of 
occurrence can be assessed, providing estimates of fractional coverage and the augmented or extra 
spectral reflectance contribution of whitecaps and foam. Measurement of wind speed/direction and 
air/water temperature will provide additional information in correlating the augmented signal to these 
formation parameters. 



Nevertheless, field testing the whitecap measurement system has provided interesting information 
on the spectral features of foam. Examples of foam reflectance, which can generally be related to 
different foam types and thicknesses, have been shown with interesting features that differ from the 
foam-tree water spectral profile. The magnitude of these features are significantly greater than the error 
due to illumination condition and solid angle response deviation of the cosine collector, filter bandwidth 
mismatch or reflectance calibration of the system. Of particular interest is the drop off between the 440 
nm and 860 nm values. These results are in good agreement with work reported by Frouin et al. (1996) 
taken in similar waters, although with lower 670 nm values (see Figs. 10 and 11) due to increased 
quantities of submerged bubbles generated by the ship. The greatest drop off is found for the thinnest 
foam layers (Fig. 1 1) where the majority of measurements in this reflectance interval are likely to be 
made of submerged bubbles. In the data acquired here, there is a greater proportion of submerged 
bubbles in the lower reflectance intervals due to the nature of the foam type generated by the relatively 
immense and hydrodynamical shape of the ship bow forcing large quantities of bubbles deep below the 
surface. As noted before, this drop off is considered to be due to the stronger absorption properties of 
water at longer wavelengths, effecting light that is reflected from submerged bubbles through a small 
layer of water between the bubbles and the surface. In the case of surf waves, the majority of 
measurements that fall within a similar reflectance interval may be from a larger proportion of thin surface 
foam, which produces an equivalent reflectance in the visible region to that of dense packs of submerged 
bubbles, but lacks the severity of the 670 and 860 nm attenuation found in ship bow foam. 

The augmented reflectance for different foam types (Fig. 1 1) shows a distinct change in the 
spectra in going from the foam-free water reflectance to the thinnest foam measurement. It also shows an 
increasing trend in the 670 and 860 nm values with respect to the other wavebands as the foam becomes 
thicker and more highly reflecting. On the other hand, the reflectance of foam (Fig. 10), without the 
background foam-free water reflectance subtracted, shows that the water color becomes more prominent 
in measurement of the thinner and lower reflecting foam types. The foam-free water reflectance 
influences the thinner foam spectra due to the inevitable inclusion of foam-free reflectance contributions. 
Since the color of the loam is inextricably related to the color of the water, the use of a normalized foam 



or whitecap spectrum at 440 nm to describe foam types in terms of percentage drop off may not be rdbust 
for every application. Even the use of normalized augmented reflectance, which appears to effectively 
remove the water color influence in these examples, can be limited when using smaller reflectance 
intervals to determine the thinnest foam from the foam-free water. 
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Table 1. Correction factors for each channel under different sky conditions due to the difference in solid 
angle response between deck cell and reflectance plaque. 
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Table 2. Correction factors for different filter bandwidths of the radiometer and deck cell in going from 
laboratory to outdoor illumination conditions. 
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Table 3. Effect of application of sky correction factor for each channel. 
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Table 4. Correction factors for each channel under different sky conditions due to deviation of solid angle 
response of the deck cell. 




































































Figure Captions: 


Fig. i. A schematic of the whitecap radiometer and its components. 

Fig. 2(a). Angular response of the deck cell and reflectance plaque for the 410 run channel. Ideal cosine 
collector response is shown for reference. 

Fig. 2(b). Angular response of the deck cell and reflectance plaque in terms of deviation from true cosine 
response for 410, 440, 5 10, 550, 670 and 860 nm channels. Each interval is 50%, and each 
waveband is offset by 50% for separation. 

Fig. 3(a). Solid angle response of deck cell and reflectance plaque for the 410 nm channel. Ideal cosine 
collector solid angle response is shown for reference. 

Fig. 3(b). Solid angle response of the deck cell and reflectance plaque in terms of deviation from true 
cosine solid angle response for 410, 440, 510, 550, 670 and 860 nm channels. 

Fig. 4. Normalized radiometer and deck cell bandpass profiles with lamp and solar spectra. 

Fig. 5. Field test of system ability to measure a plaque reflectance, as described in text. Note that the 860 
nm channel of the whitecap radiometer saturates when the plaque is illuminated by direct sun. 
Each interval represents a reflectance difference of 20%, and each waveband is offset by 20% to 
enable display. 

Fig. 6. Time series of foam data generated by bow of ship. The lower trace is for 860 nm. The rest 
(upper trace) are for other channels. 



Fig. 7. Spectral reflectance of different foam types. 


Fig. 8. Longer time series of foam data generated by ship showing periods when ship slowed and came 
to rest. 

Fig. 9. Spectral reflectance of different foam types from larger time series. 12, 14, 17, 110 and 112 refer 
to reflectance intervals normalized at 440 nm in Fig. 10. 

Fig. 10. Normalized spectral reflectance at 440 nm for different foam types taken from data used in Figs. 
8 and 9. Note the spectral influence of the background water in measurement of the thinnest 
foam types in interval 2. 

Fig. 1 1. Normalized augmented spectral reflectance at 440 nm for different foam types (used in Fig. 10). 
Note the spectral similarity for all foam types, and the increasing 670 nm and 860 nm values with 
respect to 440 nm for increasing reflectance interval from thicker fluffier foam types. 
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Abstract. A radiometric system, deployed from a ship, is used to directly measure the optical 
influence of whitecaps. Estimates of the augmented reflectance contribution and fractional 
coverage are determined from measurements in the Gulf of Mexico and open ocean waters of the 
Pacific. The accuracy of these estimates are dependent on the ability to radiometrically determine 
the difference between white water and nonwhite water events. Subtle spectral variations in the 
visible are found in the augmented reflectance but are typically smaller than the error in determining 
the white water threshold. Sky reflectance, sky condition, specular sun glitter, wave height and 
wind direction can influence the spectral characteristics of the augmented reflectance particularly at 
low sun angles. Evidence is presented suggesting that whitecaps cannot be modeled as lambertian 
reflectors, particularly at low solar elevations. The augmented spectral reflectance of whitecaps 
measured in the open ocean is found to be 0.00031 at a wind speed of 8 m/s, 0.00058 at 10 m/s, 
and 0.00096 at 12 m/s for wavebands in the visible and 15% - 25% less at 860 nm. 


Introduction 

Due to the relatively small water-leaving radiance signal from the ocean in comparison to the 
total light arriving at an ocean color sensor, the quality of data depends primarily on the accuracy of 
atmospheric correction algorithms. Atmospheric correction over the ocean is determined by 
measurement of the upwelling signal arriving at the satellite at a waveband where the radiance back 
scattered out of the ocean is known to contribute very little or no signal [Gordon and Wang, 

1994a]. Utilization of near infrared wavebands beyond 760 nm (and even 670 nm) can be justified 
in the absence of whitecaps (in case I waters) due to the strong absorption properties of liquid 
water at those wavelengths. In the case of an agitated sea state the presence of whitecaps occurring 
within a satellite image pixel augments the radiometric signal by an amount dependent on the 
whitecap coverage and the reflectance of whitecap foam [Gordon and Wang, 1994b]. This is 
referred to as the augmented reflectance of whitecaps. Previous researchers have used a whitecap 
foam reflectance of 50% - 100% [Payne, 1972; Maul and Gordon , 1975; Gordon and Jacobs , 



1977] and have assumed it to be constant at all wavelengths. Laboratory measurements carried out 
by Whitlock et al. [ 1982] yield a value of -55% in the visible part of the spectrum for dense foam 
in clear water and have shown that it diminishes with increasing wavelength beyond -0.8 pm, with 
5% decrease at 0.85 pm, 10% at 1.02 pm , and 50% at 1.65 pm. Koepke [1984] determined the 
effective reflectance of whitecaps, which takes into consideration the decrease in reflectance and 
increase in area of an individual whitecap as it ages, to be 22% in the visible. This wavelength 
independent estimate of the effective reflectance is then combined with the relative area covered 
with whitecaps as a function of wind speed [Monahan and O’ Muircheartaigh, 1980] in order to 
obtain the total optical influence of whitecaps for different wind speeds. Koepke found that it was 
necessary to introduce this value of effective reflectance since the simple combination of fractional 
coverage with the foam reflectance values of Whitlock, valid for dense fresh foam (-55%), gave 
too high a value for the optical influence of the total whitecap. 

Field measurements taken of the reflectance of foam generated by breaking waves in the 
surf zone [Frouin et al., 1996] indicate a larger spectral variation in reflectance than reported by 
Whitlock. Frouin found the reflectance decreased by 40% at 0.87 pm, 50% at 1.02 pm, and 95% 
at 1.65 pm relative to the reflectance at 440 nm. The greater spectral variation in reflectance found 
in field data, unlike the laboratory foam of Whitlock, is thought to be due to the strong absorption 
properties of water at longer wavelength acting on light reflected from submerged bubbles forced 
into the water column by the breaking waves. 

In addition to the spectral differences between foam measured in the laboratory and in the 
field, estimates of whitecap fractional coverage for different wind speeds are found to be quite 
noisy [Blanchard, 1971; Monahan, 1971; Ross and Cardone, 1974; Wu , 1979; Toba and Chaen, 
1973; Monahan and O’ Muircheartaigh, 1980,1986], Fractional coverage (W) has been typically 

related to wind speed (U) with the form W= ocU^. The values of oc and (3 vary with geographic 

location, air and sea temperature as well as wind speed range. These relationships remain noisy, 
despite attempts to fine tune particular data sets, due to the dynamic nature and interdependence of 
the many parameters involved in whitecap formation. Whitecap coverage is primarily dependent 



on wind speed, but also on other factors such as fetch and duration [Cardone, 1969], water 
temperature [ Miyake and Abe, 1948], air temperature, stability of the lower atmosphere defined by 
the air/water temperature differential [Monahan and O’Muircheartaigh, 1986], salinity [Monahan 
and Zietlow, 1969] and even surface tension variations due to the presence of organic films 
[Garrett, 1967]. A slightly different approach by Wu [1988] relates whitecap coverage to wind- 
friction velocity, being proportional to the square root of the wind stress which depends strongly 
on atmospheric stability conditions. Bortkovskii and Novak [1993] assessed whitecap and foam 
coverage with particular dependence on the sea surface temperature which effects sea water 
viscosity and wind friction. 

In addition, defining a threshold between white water and nonwhite water events is a difficult 
process and usually depends on a subjective interpretation. Further complexities can arise, such as 
the increase in foam streak-to-whitecap ratio with wind speed above 9 m/s [/toss and Cardone, 
1974], As the foam streak approximates a single bubble layer with a low reflectance contribution 
that is close the threshold, a greater error can be induced in the estimate of the augmented 
reflectance contribution at higher wind speeds. 

To date, the optical influence of whitecaps on the upwelling water-leaving radiance has been 
estimated by utilizing and combining the work of different researchers as described above. In this 
paper, a measurement system deployed from a ship is used to directly determine the augmented 
spectral reflectance and fractional coverage of whitecaps in the field. Our measurement system 
provides spectral reflectance data from individual whitecaps - as they grow and decay - and, by 
integrating over longer periods of time, estimates of the fractional coverage and augmented 
reflectance contribution to the water-leaving radiance can be determined. The principal reason we 
undertook the present investigation was to see if the Frouin etal. [1996] spectral measurements in 
the surf zone were applicable to real oceanic whitecaps. A secondary motivation was to validate 
the whitecap reflectance versus wind speed relationship proposed by Koepke [1984] [Gordon and 
Wang 1994b], 



Instrument System 

The whitecap reflectance measurement system is described in detail elsewhere [Moore et al., 
1997], The system consists of a 6-channel radiometer with narrow field of view and nominal 10 
nm bandpasses at 410, 440, 5 10, 550, 670 and 860 nm. It is held over the water surface by 
means of a boom extended from the bow of a ship providing an unobstructed view of the water 
surface. At the same time a cosine collector with matching wavebands measures downwelling 
irradiance so that the reflectance of the water surface can be calculated. A TV camera is mounted 
by the radiometer to provide a visual reference, and whenever the radiometric data is acquired the 
associated video frame is time and date marked and recorded onto video tape to assist in later 
analysis. Air and water temperature, wind speed and direction are measured simultaneously with 
the radiometric data at a rate of ~7 times a second continuously for about 30 seconds (providing 
200 contiguous samples), after which global positioning (GPS) data, universal time, and location 
are recorded. This acquisition sequence is repeated until a time determined by the operator. This 
radiometric rate allows many sample points of an individual whitecap to be captured as well as 
providing an adequate data set over a time period of reasonably consistent sea state and sky 
conditions. 

Once the system has been installed on a ship, only periodic measurements of the dark current 
are required. This requires bringing the radiometer in from the boom and covering both the 
radiometer and the deck cell with light tight caps. This is necessary only when significant 
temperature changes are encountered such as when conditions change from early morning to noon, 
or from dark, overcast, rainy conditions to bright sunshine. The radiometer can also be angled up 
to 20° from nadir in order to minimize the effects of specular sun glitter and bright sky reflections. 

Data Description 

Data from two different cruises are presented in this paper. The data acquired from the first 
cruise is used to introduce the data reduction process, consider the difficulties in radiometrically 
defining white water events and to investigate the effects that low solar elevation, wind direction. 



specular sun reflection and sky conditions have on the estimate of fractional coverage and on the 
augmented spectral reflectance contribution. In the second cruise, data was acquired in the open 
ocean and is used to provide practical measurements of the augmented reflectance contribution 
needed for ocean color and atmospheric correction algorithms. 

For the First cruise, the whitecap radiometer system was deployed in the spring of 1996 on the 
NOAA ship, R/V Malcolm Baldrige for a 2 1-day ship cruise from Miami to a test location in the 
Gulf of Mexico, approximately 70 miles offshore from Cedar Key (Florida) in the Apalachicola 
Bay. The location provided relatively warm waters (-16 - 17° C) with a number of cold fronts 
moving down off the continental mainland, resulting in variable conditions. These fronts usually 
lasted a couple of days bringing strong winds (sometimes as high as 18 m/s) and lowering the air 
temperature to about 12°C. The occurrence of an unstable atmosphere and good winds provided 
an interesting whitecap data set. 

The second cruise involved a -6,000 km transit through equatorial waters of the Pacific Ocean 
from Manzanillo, Mexico to Honolulu, Hawaii. This 13-day cruise took place in November 1996 
abroad the NOAA ship R/V Ka’imimoana. Conditions were more consistent with air/water 
temperature and wind speed/direction changing only by a relatively small amount over time. Fetch 
was typically greater than 1000 km and duration could be considered unlimited. Data was acquired 
far from land and for the most part, water temperature was ~22-23°C with air temperature -20- 
23°C. Wind speeds ranged from -8 m/s to -13 m/s and wind direction changed from an initial 
northerly direction at the beginning of the transect to an easterly and southerly direction as the 
Trade winds were crossed. The sky conditions were mainly broken and scattered with some 
overcast periods. 

The data used in the following section was acquired during one particular cold front passage 
during the First cruise in the Gulf of Mexico. The front approached in the morning with a 
completely overcast sky. Wind speed was initially -6 m/s from the south east but at about 10:40 
am the wind changed direction to the north northwest with wind speed increasing to about 15 m/s 
(gusting to 18 m/s) within 15 minutes. The water surface which had been relatively ‘flat’ with a 



wave height of about 1-2 ft (30 - 60 cm) became ruffled with many small whitecaps appearing. 

As the day progressed the whitecaps became larger but fewer in number as the wave height 
increased and a swell began to develop. Although there were lull periods, the winds picked up 
again during the night blowing between 10 to 15 m/s and changing direction to the east as the front 
passed. By the next day, the wave height had increased to 2 m with a swell frequency of 7-8 
seconds. The sky was clear and the ship adopted a rocking motion with the sea swell. Examples 
of data taken from both days with similar wind speeds but differing sea states, sky conditions, air 
temperatures, duration, etc. are presented in the next section. 

Data Reduction 

Phase I Reduction Process 

One half hour of data collected at the rate described above amounts to about a 1.5 megabytes 
of data storage. To reduce this data, analysis software using Madab was developed to partition and 
separate the data matrix into radiometric data, environmental parameters (wind speed/direction, 
air/water temperature) and GPS data. The radiometric data is then reduced to produce dark 
background-subtracted and radiometrically-corrected reflectance values from the six up- and 
downwelling channels. This radiometric correction process consists of three steps [Moore et ai, 
1997] briefly reviewed here. 

First, ratios of upwelling radiance to downwelling irradiance are formed from the dark 
subtracted signals for each channel. These six ratios are then multiplied by a calibration factor 
determined in the laboratory relating the channels in the radiometer to those in the deck cell in order 
to yield reflectance measurements. This was established using a calibrated 1000W (FEL) quartz 
halogen lamp source and a calibrated Spectralon reflectance plaque. The deck cell was positioned 
50 cm in front of the source and measurements taken. For the radiometer, the reflectance plaque 
replaced the deck cell, and the radiometer viewed the plaque from a 45° angle. 



For the second step, a multiplication factor is required to correct bandpass differences between 
the radiometer and deck cell and for differences in the light spectrum between laboratory and 
outdoor illumination conditions. 

Thirdly, each channel of the deck cell must be corrected for solid angle response for 
operation outdoors under 2n illumination due to the slight deviation in angular response of the deck 
cell cosine collector from true cosine. Correction factors for different angular distributions of the 
downwelling light Field are required for various sky conditions that may be encountered. A set of 
correction factors were established for different sky conditions such as overcast, or clear sky with 
different solar zenith angles, to correct solid angle response. 

After applying these corrections, the corrected data are presented as a time series in blocks of 
200 water surface reflectance measurements along with wind speed/direction, air/water 
temperature, upwelling radiance and downwelling irradiance data. The reflectance data are then 
grouped into reflectance intervals; any reflectance spectra whose average value falls within a 
predetermined reflectance interval is binned. In analysis of foam data obtained with the system in 
earlier work [Moore et al.., 1997], the value of the reflectance at 860 nm was used as the binning 
parameter rather than the mean spectral value used here. However, in estimating the fractional 
coverage and augmented spectral reflectance, it was found that it makes no difference which 
parameter is used. 

For analysis of the data from the first cruise, the reflectances are divided into 48 intervals, 
with high resolution reflectance intervals of 0.0025 reserved for reflectance values between 0.005 
and 0.1. (First interval is 0.0 to 0.005). Above this, the reflectance intervals are incremented by 
0. 1 and continue up to 1.0. As will be seen later, the high resolution intervals used for the lower 
reflectance values below 0. 1 have been implemented in an attempt to accurately determine and 
discriminate between white water and nonwhite water events. Reflectance intervals above 0. 1 
were primarily used to analyze the spectral characteristics of different foam types which may be 
defined by their reflectance magnitude. 



In the work by Frouin et al. [1996] in making spectral measurements of foam in the surf 
zone, reflectance intervals were also used in an attempt to determine spectral variations with foam 
type and thickness. The total reflectance was partitioned into 4 intervals (0 - 0.25, 0.25 - 0.4, 0.4 
- 0.55, and any reflectance greater than 0.55). Koepke [1984] on the other hand required a 
threshold between nonwhite water and white water in analysis of photographic images to determine 
fractional coverage and the effective reflectance, which considers the drop in reflectance and 
increase in area as a whitecap decays. To accomplish this, the white water threshold was set at 0. 1 
with the maximum reflectance in a particular photograph normalized to 1. These images were 
acquired using a 6 cm x 6 cm Hasselbad camera angled between 45° and 60° degrees in viewing 
the surface from a height of 30 m. This threshold was further complicated by the increased white 
threshold towards the horizon due to increased Fresnel sky reflection. With the whitecap 
radiometer used in this paper, the effect of sky reflection variation with viewing angle is drastically 
reduced due to the small field of view of the radiometer optics (1° ), and with 0.0025 reflectance 
resolution more accurate threshold information can be obtained. 

After the reflectance data from each block has been grouped into reflectance intervals, these 
reflectance intervals are added to previous blocks of similarly reduced reflectance data, resulting in 
a set of reflectance intervals for the complete data matrix. UTC time is extracted and displayed at 
the beginning of each block so that reflectance measurements can be correlated with water surface 
events replayed on the VCR. The VCR (Sony EV-C100 (video Hi8» can be advanced frame by 
frame in either direction with the time/date mark appearing on the bottom of the monitor 
approximately every third frame. In this way, individual samples from whitecaps and their 
different foam types can be radiometrically tracked. 

Once all the data from a particular acquisition period has been reduced to this level, the spectral 
reflectance values and the number of samples binned in each interval are stored along with the total 
number of samples taken, mean air/water temperature, wind speed/direction readings and GPS data 
from each block (200 samples) of the original data matrix. 



Phase I Examples 

In Fig. 1(a) an example of a time series with two whitecaps passing under the radiometer is 
shown. These whitecaps were acquired during the onset of the cold front described above. The 
sky was overcast and the data were taken ~1 1:00 am. Wind speed was -12 m/s, air temperature 
16.2°C and water temperature 17°C. The 96 samples (from a block of 200) shown in Fig. 1(a) 
were acquired over a period of - 15 seconds. In this example, a large whitecap suddenly breaks 
into view of the radiometer with thick white foam (sample point 1 1) reaching a peak reflectance of 
-55%. Six traces are plotted representing the six radiometer channels (410, 440, 5 10, 550, 670 
and 860 nm). The lower trace (symbols) corresponds to the 860 nm reflectance. Thick foam is 
temporarily replaced by a region of submerged bubbles and thinner surface foam (-sample points 
13, 14 ,15). Thick foam comes into view again at sample point 17. At sample point 20 and 21 a 
thin layer of foam passes followed by some dense (not thick) older foam and stays in field of view 
to about sample point 35. Sample points from about 35 - 56 show the reflectance of thinning 
residual and fragmented surface foam. From 60 to about 75 the reflectance of the foam free water 
surface is shown and is suddenly followed by another whitecap of smaller magnitude (sample 
point 76) which continues to decay out to about sample point 96. 

In Fig. 1(b) the individual reflectance samples from this sample period have been binned into 
reflectance intervals. The lower reflectance intervals correspond to the nonwhite water reflectance; 
some are higher due to Fresnel reflectance from brighter regions of the sky, others are lower due to 
shadowing by adjacent waves. The variation in the spectral profiles of each interval is partly due to 
the small number of samples used. However, a common feature that stands out is the lower 860 
nm reflectance. This is attributed to reflectance from submerged bubbles modified by the strong 
absorption of water at 860 nm [Frouin et al., 1996; Moore et al„ 1997]. (Other examples of phase 
I data reduction into reflectance levels are shown in Figs. 2(b), 3(b) and 4(b) where larger data sets 
of about 10,000 samples are used). 



Phase II Reduction Process 


In phase II, the reduced data from phase I are processed to determine fractional coverage and 
the augmented spectral reflectance contribution. Correlations between these estimates and the 
environmental conditions responsible for whitecap development can then be investigated. 

Each reflectance interval, resulting from the reduction of a complete data matrix, is multiplied 
by the fraction of reflectance samples that have contributed to that reflectance interval. This 
weighted reflectance spectrum is the effective spectral reflectance contribution of that reflectance 
interval (i.e. the contribution of that interval to the total reflectance). Repeating this for all 
reflectance intervals and summing their effective reflectances together yields the total spectral 
reflectance. This is equivalent to the reflectance observed in a satellite image pixel mapped onto the 
ocean surface. 

The total reflectance contains both white water and nonwhite water contributions. In order to 
determine only the contribution of whitecaps and foam, the background water-leaving radiance, 
sky reflectance, sun glitter and other nonwhite water contributions must be removed. The 
background water color (and sky reflection etc.) has the greatest number of samples being the most 
frequently measured and since the presence of foam on the surface of the water increases the 
overall reflectance signal, the optimal threshold must be at some reflectance interval above it 
However, reflectance intervals above the reflectance interval with the most number of samples may 
contain both thin foam reflectance as well as sky reflectance contributions from brighter regions of 
the sky. Also, reflectance intervals below this may contain thin foam reflectance contributions that 
are in the shadow of an adjacent wave and may register as having a lower reflectance than a sky 
reflectance contribution from an unshadowed foam free water facet. This is true for threshold 
intervals near the reflectance interval with the most number of samples. This inability to 
distinguish between white water and nonwhite water events places a limit on the accuracy with 
which fractional coverage and augmented spectral reflectance contribution can be determined. 

The estimate of fractional coverage and the augmented reflectance contribution is dependent 
on choosing a reasonable reflectance threshold. Once an appropriate threshold has been chosen, it 



is subtracted from each interval above it. This difference in reflectance is assumed to be due to 
white water contributions only. The reflectance difference is then multiplied by the fraction of 
samples that have fallen into that reflectance interval to yield the effective augmented spectral 
reflectance of that particular reflectance interval. All the weighted reflectance intervals above the 
threshold are summed to yield the (total) augmented spectral reflectance contribution. Fractional 
coverage is estimated from the total number of white water samples divided by the total number of 
samples and, notably, is not weighted by the magnitude of reflectance. 

In the following examples, a range of plausible thresholds are implemented to demonstrate the 
effects on the estimate of both fractional coverage and the augmented spectral reflectance 
contribution from data acquired with similar wind speeds but under different sky conditions, solar 
zenith angles, and sea states. 

Phase II examples 
Overcast Sky 

Data presented in Figs. 2(a) - 2(d) were the result of 9,000 surface reflectance measurements 
acquired over a period of -23 minutes as the ship covered a 4 kilometer path during the onset of the 
cold front described above. Mean location was 29.05°N and 83.5 1°W. During acquisition the 
mean wind speed was 1 1.86 m/s with a standard deviation of 3.6 m/s (gusts up to 18 m/s). The 
wind was from the north between 350° and 0°. Ship heading was originally -170° during the first 
half of the acquisition period and then passed through a southerly direction and settled on a north 
westerly direction of 330°. The radiometer was nadir viewing and since conditions were overcast 
there was no specular sun glitter. Air temperature was 16.2°C and water temperature 17.0°C, 
resulting in a slightly unstable atmosphere. Other parameters taken from the ship’s log were 
salinity (35.027 ppt), barometric pressure (1015.4 mB), and relative humidity (93%). In this 
example, during the onset of the cold front, many small whitecaps suddenly appeared on the 
relatively flat and dark water surface. Duration was limited to -15 minutes with whitecap 
development and sea state measured during an abrupt transition period in weather conditions. The 



resulting measurements produced high values of the augmented reflectance and fractional coverage 
estimates. 

In Figs. 2(a) - 2(d) the more interesting portions of the output of data processed in phase II are 
presented. Figure 2(a) shows all the reflectance samples reduced to reflectance intervals. In Fig. 
2(b) a histogram of the number of samples used to determine each reflectance interval is shown. 
Remember that the First 39 levels cover the mean reflectance spectral values between 0 and 0.1. 
Above that the reflectance intervals are incremented in 0. 1 steps causing the bump at interval 40. 

In the histogram there seems to be no obvious distinguishing feature to identify the background 
from the white water, but as mentioned above, a suitable threshold must be above the reflectance 
interval with the most number of samples. 

Figure 2(c) shows the total reflectance. The total reflectance spectrum peaks at 5 10 nm 
suggesting a blue-green color and a maximum value of -0.06180. The nonzero value at 860 nm 
(0.02412) is mainly due to the overcast sky being reflected at the water surface. The graphs in Fig. 
2(d) show the augmented spectral reflectance contribution and an estimate of the fractional 
coverage for different thresholds choices. 

In Fig. 2(d) a range of reflectance intervals (intervals 16 through 23 and where interval 16 
contains the greatest number of samples) are used as threshold reflectance intervals between white 
water and nonwhite water events. In the resulting estimates of the augmented spectral reflectance, 
the spectral shape has changed from that of the total reflectance. The spectral shape of the 
augmented reflectance seems to vary somewhat from one choice of threshold interval to another 
with exception of a common lower 860 nm value and a relatively flat region between 440 and 670 
nm. The 860 nm value, as seen in previous data, is much lower than the other wavelengths. 

For threshold intervals 16, 17 and 18 the fractional coverage estimates are particularly high 
(37.2%, 22.5% and 14.3%) for this wind speed of 1 1.86 m/s. From threshold interval 19 through 
23, the spectral shape has become more settled and appears to be relatively constant with the 
exception of the 410 nm channel which appears to decrease with increasing threshold with respect 
to the other spectral values. In this range of thresholds, the estimates of fractional coverage are 



also more plausible (9.5%, 6.9%, 5.3%, 4.1% and 3.7%). This spectral shape continues for 
higher threshold intervals (not shown here) and eventually breaks down with the highest threshold 
intervals due to a small number of samples. 

It is interesting to note that the augmented spectral reflectance for these threshold intervals ( 19 
through 23) maintain a small peak at 5 10 nm although it decreases with increasing threshold. 
Increasing the threshold moves the analysis into higher reflectance levels and therefore into thicker 
whitecap foam types. Investigation of the spectral variation of different foam types with 
increasing reflectance [Moore et al., 1997] revealed a similar trend with the 5 10 nm value 
decreasing while the 670 nm and 860 nm values increased. The 410 nm value was also observed 
to decrease in the foam spectra. 

The 5 10 nm reflectance maintains the highest value and is reminiscent of the total reflectance 
spectrum which is predominantly comprised of non-foam water reflectance contributions. The 
reflectance values of the wavebands in the visible region tend to equalize with higher threshold 
choice. The augmented reflectance spectral profiles seem reasonable and explain the apparent 
whiteness of whitecaps, where the channels covering the visible region are essentially equivalent 
with an increasing hint of blue-green for lower reflectance intervals. This is similar to what is 
visually observed; white thick fresh foam and blue- green thin foam where the background water 
color penetrates through the thin bubble layers. 

With a consistent spectral profile for the higher threshold intervals, the optimal choice 
becomes difficult. However, it seems that threshold interval 23 is perhaps too far into the white 
water realm and the can be ruled out as the optimum threshold choice. Threshold interval 19 might 
be a reasonable choice but could contain some non-foam bright sky reflectance contributions. It 
also has a particularly high 410 nm value that is not present when higher threshold reflectance 
intervals are applied suggesting the presence of some non foam contributions. Being conservative, 
the optimal choice seems to be around threshold interval 20 or even 21. The augmented 
reflectances are 0.003 15 and 0.00290 at 5 10 nm and 0.00169 and 0.00161 at 860 nm for 
thresholds 20 and 2 1 respectively. The fractional coverage estimates using these thresholds. 



although a little high due to the inclusion of the thinnest foam types, arc 6.9% and 5.3% 
respectively. 

In the augmented spectra for thresholds 20 and 21, changing the threshold by a small 
reflectance of 0.0025 changes the fractional coverage estimate by 1.6% but only by -0.000186 for 
the augmented reflectance in the visible (410 nm -550 nm), and by -0.00008 for the augmented 
reflectance at 860 nm. 

The upwelling signal, particularly in the lower threshold regions, is mixed with contributions 
from brighter sky reflectances, thin foam layers and fragmented foam patches. As stated before, 
there is no decisive definition of threshold; only a comprise between either obtaining a reasonably 
pure white water spectrum and underestimating the augmented reflectance contribution, or 
obtaining perhaps a more complete measure of the augmented reflectance contribution at the cost of 
including some nonwhite water contributions. It is important to include the lowest reflecting and 
thinnest foam types since they persevere longer than the thicker and more highly reflective foam 
and as result have great statistical weight in contributing to the total augmented reflectance. This 
threshold problem makes accurate determination difficult and somewhat subjective. However, it is 
comforting that the augmented reflectance is a weak function of the actual choice of the threshold. 

Clear Sky 

In the next example, whitecap data was taken the next day (-9:05 am) under clear sky 
conditions, and as described above, with a more developed sea. Wave height was about 2 m with 
a periodicity of 7-8 seconds. Mean wind speed was 1 1.7 m/s with a standard deviation of 1 m/s 
coming from the north west at about 340°. The ship heading was 355° and moving through the 
water at 7 knots (2.6 m/s) covering a straight track of 4.3 km. Air temperature was 12°C and the 
water surface temperature was 16.5°C, resulting in an unstable atmosphere. Other parameters 
taken from the ship’s log were salinity (35.0895 ppt), barometric pressure (1019.9 mB) and 
relative humidity (61%). This location (-29. 15° N, 84.02° W) was approximately 20 km north and 
56 km west of the location where the data during overcast conditions in Fig. 2 was acquired. With 



wind speed continually blowing in the 10-15 m/s range for roughly a 24 hour period, the sea had 
become more fully developed, producing a different whitecap morphology and foam distribution 
from the previous day. Reflectance measurements included a larger proportion of foam streaks and 
residual foam patches resulting from larger breaking waves although less frequent than in the 
previous ‘overcast’ case. 

In Fig. 3(a) the reflectance measurements reduced to reflectance intervals are shown and the 
histogram of reflectance interval vs. number of samples is shown in Fig. 3(b). Noticeable 
differences from the overcast sky condition are (i) the tighter distribution of samples about the 
peak, and (ii) the much lower sky reflectance (-0.00276 at 860 nm). This is the result of a low 
sun angle (30° above the horizon) and the viewing angle of the radiometer. The radiometer was 
orientated to view the water surface at 20° from nadir in a direction opposite the sun in order to 
minimize specular sun glitter in the measurements. As a result, the sky reflectance signal arriving 
at the radiometer was from a region of the sky centered roughly 80° from the sun where the sky is 
dark relative to overcast conditions. 

With regard to the histogram of reflectance interval vs. number of samples, it may be tempting 
to fit a symmetrical distribution about the mode and use the difference between the fit and the 
histogram to identify any augmentation. But this would assume an equal probability of darker and 
brighter sky reflectances contributing to the distribution from a symmetrical distribution of wave 
orientations of the water surface reflecting a proportionately gradiented sky. It is not that simple 
since the wind typically blows from one direction causing the waves to progress across the surface 
in rows creating a preference for surface angle orientation and skewing any symmetry in the 
distribution. In addition, in the presence of low sun angles, the distribution of reflectance intervals 
is further complicated by the specific illumination direction of the sun reflected from waves and 
whitecaps. Even if the distribution of reflectance intervals vs. number of samples were perfectly 
symmetric, the finite bandwidth of this distribution indicates an accuracy limitation. Within this 
envelop there is an indiscriminate region were sky reflectance signal overlaps with thin foam 
contribution thus concealing any obvious threshold. 



In Fig. 3(c) the total reflectance has a peak at 5 10 nm but is not as prominent as in the overcast 
data. This is most likely due to the ship being in a slightly different body of water. In Fig. 3(d) 
the augmented spectral reflectance profiles are shown for thresholds 4 through 9 (threshold 4 is the 
mode) . Following similar reasoning as in the overcast case, it seems that the profiles resulting 
from the application of thresholds 6 and 7 are more realistic. The augmented reflectances are 
0.00097 and 0.00085 at 5 10 nm and 0.00071 and 0.00059 at 860 nm for thresholds 6 and 7 
respectively with fractional coverage estimates of 4.8% and 3.4% respectively. 

It is interesting to note the relatively higher 440 nm value appearing in the augmented 
reflectance spectra as it does in the total reflectance profile, whereas in the Fig. 2 data it is relatively 
low in both the total and augmented reflectance spectra with the application of higher thresholds. 
This difference is due to measurements taken in slightly different water bodies and to the inclusion 
of the foam free water color (sky reflection and water-leaving radiance) in measurements of the 
thinnest foam types and submerged bubbles. 

In this example, the augmented reflectance at 670 nm using threshold 6 is marginally highest 
and the 860 nm has substantially increased in comparison with the overcast example. Also, the 
augmented reflectance contribution at 670 appears to be consistently higher than at 550 nm, 
whereas in the overcast example both the 550 and 670 nm values are essentially equal. The 
augmented reflectance profile suggests a stronger influence from the longer wavelengths. This 
effect is even more pronounced in the following example. 

Clear Sky Low Sun 

In Figs. 4(a) - 4(d), data was acquired an hour earlier (at -7:50 am) on the same day as the 
data presented in Fig 3. The sun was lower in the sky at approximately 10° above the horizon and 
at -80° azimuth from true north. Mean wind speed was 12. 1 m/s with a standard deviation of 1.3 
m/s during the acquisition period. The wind was out of the north west (325°) and the ship 
heading was 355°. Air temperature was 12.6°C and water surface temperature 16.55°C (unstable 
atmosphere). Other parameters from the ship’s log were salinity (35.0445 ppt ), barometric 



pressure (1019. 1 mB) and relative humidity (62%). As in the preceding example, the radiometer 
was angled 20° away from the direction of the sun to reduce the possibility of sun glitter and 
specular reflection. The location (29. 1 1° N, 84.02° W) was approximately 2.7 miles due south of 
the location where the data in Fig. 4 was taken. During this period of ~18 minutes, 7,000 surface 
reflectance samples were taken. The whitecap development and foam distribution during this 
acquisition period was similar to the ‘clear sky’ case, only the sun elevation was smaller. 

Figure 4(a) shows the reflectance intervals. In the histogram of the number of samples vs. 
reflectance levels in Fig. 4(b), the shape of the distribution is similar to that in Fig. 3(b). The 
distribution rises sharply and falls with a slight broadening towards higher reflectance. The total 
reflectance shown in Fig. 4(c) is somewhat higher than in Fig. 3(c) data suggesting a possible 
higher white water contribution. The spectral shape is similar with the exception of relatively 
higher 410 and 440 nm values. This can be attributed to a spectrally different angular distribution 
of the downwelling light field. The sky was comprised of a prominently yellow reddish maximum 
at high zenith angles due to the low sun, and a very rich and dark blue sky toward the zenith - 
typical of a clear early morning as the sun rises. Horizontal waves and waves inclined away from 
the low sun reflect the predominantly very blue sky and are most frequently measured by the 
radiometer. This bias towards more blue measurements accounts for the different shape of the total 
reflectance profile in this example. However, when plausible threshold reflectance intervals (Fig. 
4(d)) are applied in order to estimate the augmented reflectance contribution, the opposite is 
observed with a bias towards the longer wavelengths. 

At threshold interval 5 (mode of the histogram) in Fig. 4(d), the augmented reflectance values 
at the shorter wavelengths are lower with respect to the longer wavelengths than they were in the 
total reflectance profile. Using the same criteria from the two previous data sets, white water 
contribution should begin to appear above the background reflectance in the region of reflectance 
intervals 6, 7, 8, 9 and maybe 10. Moving up to threshold interval 6, the resulting augmented 
reflectance spectrum is more accentuated with even lower values towards the blue end and higher 
values in the red. This is the result of the low sun angle ( 10° above the horizon). The augmented 



reflectance spectra for threshold 6 yields a 410 nm value (0.00203) that is less than that at 860 nm 
(0.00265), and with a maximum at 670 nm (0.00313). 

Although the spectral shape of the augmented reflectance profiles from thresholds 6, 7 and 8 
are somewhat similar, their fractional coverage estimates appear to be too high - 18.5%, 10.9% 
and 7.9%. Threshold 9 has a more reasonable fractional coverage of 6.0% with an augmented 
reflectance contribution -0.0017 in the visible. Threshold 10 has a fractional coverage of 4.8% 
and an augmented reflectance of -0.0015 in the visible. Comparing the fractional coverage and 
augmented reflectance contribution estimates with the previous example taken under similar 
conditions (wind speed, air/water temperature, sea state) but with higher sun angle, shows the 
estimates in this example are much higher; nearly 3 times the augmented reflectance and more than 
twice the fractional coverage using similar thresholds. These augmented reflectance profiles and 
their particularly high estimates of fractional coverage suggest that the higher augmented reflectance 
values at the longer wavelengths may include specular reflections from foam free water facets, 
although an additional explanation is provided below. 

Comparing Figs. 3(d) and 4(d) shows that the augmented reflectance becomes relatively more 
enhanced in the red and near infrared for the lower solar elevation. Although whitecaps are usually 
modeled as lambertian reflectors, we believe that this reddening of their reflectance with reduced 
solar elevation provides indirect evidence of their nonlambertian nature. As the solar elevation is 
decreased, the contribution to the incident irradiance from the diffuse sky light is enhanced at the 
expense of the direct sun light For low solar elevations, the incident irradiance in the blue is 
almost totally from diffuse sky light. In contrast, the reverse is true in the red, most of the 
incident irradiance is direct sun light. Let us assume that the whitecaps are lambertian, i.e., they 
have reflecting properties similar to sheets of white paper floating on the water, and that their 
reflectance is spectrally neutral throughout the visible. In the red, whitecaps on the side of the 
waves away from the sun or in the troughs of the waves receive little or no illumination, i.e., they 
are shadowed from the direct sun. In contrast, in the blue, where the illumination is diffuse, all 
whitecaps will be illuminated. Thus, in this example, the augmented reflectance would be expected 



to increase from red to blue because fewer whitecaps are illuminated with red light. In fact, the 
experimental result [Fig. 4(d)] is exactly the opposite! 

Further evidence of the nonlambertian character of whitecaps can be seen by comparing Figs. 
3(a) and 4(a). The highest reflectance levels in these figures will be due to the thickest-youngest 
foam patches. In this case, removal of the ocean background is not relevant to measuring their 
properties. The reddening of the spectrum in the visible for these foam patches from Figs. 3(a) to 
4(a) is very evident Figure 4(a) shows a strong variation in reflectance with wavelength, a 
surrogate here for the variation in the angular distribution in of the incident irradiance. In contrast 
Fig. 3(a) shows a much more neutral or bluer reflectance for the most-reflective whitecaps. Thus, 
even for the brightest whitecaps, the dependence of reflectance on the angular distribution of the 
incident irradiance is clear: they are nonlambertian. However, note that the measurements we have 
presented here cannot shed light on the actual bi-directional reflection distribution function (BRDF) 
of the augmented reflectance. The only conclusion we can make regarding the BRDF is that the 
nadir-viewing augmented reflectance spectrum for diffuse illumination [Fig. 2(d)] and for solar 

elevations > about 30° are similar. 

The nonlambertian nature of whitecaps and its relationship to the solar elevation is easy to 
understand. Whitecaps are isolated volume scatterers, and at low solar elevations the solar beam 
can enter the sides as well as the top of the whitecap and be scattered out toward the zenith. In 
contrast, for high elevations or for diffuse illumination relatively fewer photons will enter the 
sides, and this results in a lower overall reflectance. 

Open Ocean Data 

The data used here were acquired from a transect of the Pacific Ocean from Manzanillo, 
Mexico (19.03° N, 104.20° W) to Honolulu, Hawaii (21.20° N, 157.55° W). Only data which 
were acquired under relatively constant and overcast conditions were used to model augmented 
reflectance with other parameters such as wind speed etc. These data sets have a well behaved 



background in terms of radiometrically determining a threshold between white water and nonwhite 
water events. Overcast conditions also provide a better downwelling light description for 
correction of the deck cell cosine collector [Moore et al., 1997]. Other whitecap data were acquired 
during broken sky and scattered cloud conditions, but have not been included in the modeling. 

The threshold is more difficult to define for these data sets, due to the large variability and random 
nature of the downwelling light field and its effect on surface reflectance measurements. 

The data used are from 4 different days and locations during the transect Wind speed was 
relatively constant and remained within a range of 1-2 m/s for days at a time. Reflectance data was 
acquired for wind speeds between -8-13 m/s. Stability of the lower atmosphere, defined by the 
air temperature / water temperature differential, varied from neutral to slightly unstable. Other 
parameters such as salinity, humidity and barometric pressure were also recorded for possible 
correlation. 

The reflectance interval resolution was increased from the previous value of 0.0025 to 0.001. 
While more threshold choices were generated as result of the higher resolution, it was primarily 
implemented to monitor background stability by observing the distribution of reflectance intervals 
verses number of samples. In cases where the sun came out from behind the clouds or a cloud 
blocked the sun, the foam free water surface reflectance either decreased or increased. This was 
observed as a broadening (sometimes with bumps) in the distribution of reflectance interval vs. 
number of samples for given data set. 

In Fig. 5 the augmented reflectance contribution from 17 data sets is presented as a function 
of wind speed. Each of the 17 augmented reflectance data points were acquired from an average 
acquisition time of 45 minutes (or -800 m 2 surface area). For clarity, the augmented reflectances 
from wavebands 410 nm through 550 nm (ARC 410 . 5J0 ) have been averaged since their values vary 
only slightly with respect to one another. The augmented reflectance at 670 nm (ARC^) which 
also varies only slighdy with the other visible wavebands, and the augmented spectral reflectance at 
860 nm (ARCg^) which is always distinctly lower than the visible wavebands are shown. To 
relate augmented reflectance to wind speed a power law relationship has been applied in the same 



vein as previous authors have related fractional coverage to wind speed [ Blanchard , 1971; Ross 
and Cardone , 1974; Wu, 1979; Toba and Chaen , 1973, Monahan and O’ Muircheartaigh , 
1980,1986], For a given wind speed the augmented reflectance contribution in the 410-550 nm 
region can be related to the wind speed with a correlation of 0.702 by: 

ARC 4[0 . 550 = 9.648 x 10' 7 U 1777 

where ARC 410 . 550 is the augmented reflectance contribution as described above and U is the wind 
speed in meters per second. Thus, the wind speed accounts for about half of the variance in 
ARC 410 _ 55 , For the augmented reflectance contribution at 670 nm and 860 nm we have: 

ARC 670 = 1.131 x 10' 6 U 2 698 

and 

ARC S60 = 1.302 X 10‘ 6 U 2 545 

with correlation factors of 0.699 and 0.647, respectively. Data taken during scattered and partly 
cloudy conditions, although not used in determining the above relationship, are also shown in Fig. 

5. These data were acquired with sun elevation greater than 50°. Reducing this data was more 

complex due to a variable threshold but the results are similar to the data acquired under steady 
overcast conditions with a well defined threshold. 

In Fig. 6 the respective fractional coverage estimates are plotted as a function of wind speed. 
The distribution of fractional coverage estimates for a given wind speed appear to be more scattered 
than for augmented reflectance. 

In Fig. 7 the relationship between augmented reflectance contribution (ARC 410 . 550 ) and 
fractional coverage can be expressed as a linear relationship with a correlation of 0.694 given by: 



ARC 41 , M5(j = 9.252 x 1() 3 W + 0.(XX)195 
where W is the fractional coverage. 

A scatter plot of the percentage drop in augmented reflectance from ARC 410 . 550 to ARC S60 is 
shown in Fig. 8. A linear fit is applied between the data points to indicate an increasing trend in 
the drop off between ARC 410 . 550 and ARC S60 with wind speed rather than to establish a model, as 
the correlation is very low, 0.359. However, the data indicate that for wind speeds in the range 8- 
12 m/s the percentage drop off between ARC 410 . 550 and ARC g60 is approximately 15 - 25%. 

Possible correlations between the augmented reflectance and other parameters such as water 
temperature, air temperature, stability of the lower atmosphere, humidity, salinity, barometric 
pressure, wave height, sea swell were investigated. No correlations could be established mainly 
due to the small variation in the parameters during the transect. 

Discussion 

In this paper we have shown examples of whitecap spectra of different foam types and 
reflectances, both from individual whitecaps and as an integrated sum over a period of time 
expressed as the augmented reflectance contribution and fractional coverage. From the first cruise, 
estimates of the augmented spectral reflectance contribution and fractional coverage have been 
obtained from measurements taken under different sky conditions and sea states but with similar 
wind speeds in order to highlight the variability of whitecap formation and the difficulties in 
providing accurate information on the optical influence of whitecaps. From the second cruise, 
estimates of augmented spectral reflectance and fractional coverage acquired during overcast skies 
and more consistent weather conditions with essentially unlimited fetch and duration have been 
derived for ocean color satellite waveband correction and atmospheric correction algorithms. 

Of particular importance is the determination of the threshold required to accurately determine 
augmented spectral reflectance contribution and fractional coverage. Despite the application of high 



resolution white water threshold reflectance intervals, the threshold choice remains difficult and 
somewhat subjective. Regardless of the threshold resolution that may be applied, the very wave 
nature of the water surface and the distribution of the downwelling light field creates an inseparable 
and indiscriminate region between foam and foam free contributions. As mentioned earlier, the 
reflectance of foam in the shadow of an adjacent wave may result in a lower reflectance than a foam 
free water facet reflecting a bright region of the sky. Choice of threshold is further compounded 
since the thinnest foam types, which are close to the threshold reflectance, are very important due 
to their long endurance. Without adequate discrimination of the thinner foam type, the estimation 
of fractional coverage and augmented spectral reflectance contribution will always incur some 
error. Utilization of a higher threshold resolution and the subsequent inclusion of the less 
reflecting and thinnest foam types results in fractional coverage estimates that are roughly 3 times 
higher than that reported by other researchers [ Blanchard , 1971; Ross and Cardone , 1974; Wu, 
1979; Toba and Chaen , 1973; Monahan and O’ Muircheartaigh, 1980,1986]. The ratio of the 
augmented reflectance contributions in the visible to the fractional coverage estimates determined 
with this system results in an effective reflectance of - 1.25%. This much lower effective 
reflectance than the 22% determined by Koepke [1984] gives an indication of the greater quantity 
of thin foam and foam streaks that have been included in these measurements. 

From the data used in this paper, the fractional coverage has been observed to vary by a larger 
amount than the augmented reflectance contribution in going from one threshold interval choice to 
another. This is due to the definition of fractional coverage which is not weighted by reflectance, 
and as result, the thinnest foam patches are assumed equal in reflectance to the more highly 
reflecting thick, dense foam types. The augmented spectral reflectance contribution is the more 
meaningful measurement for satellite observations whereas fractional coverage must be combined 
with reflectance data to achieve the same result. 

Examples from the first cruise in the Gulf of Mexico reveal the dependence of whitecap 
formation on parameters other than the wind speed measured at the time of observation. The 
examples described from the first cruise all had similar wind speeds. Changing weather conditions 



such as during the onset of a cold front produced higher values of augmented reflectance and 
fractional coverage (‘overcast’) than when conditions had been consistent for a longer period of 
time (‘clear sky’, ‘clear sky low sun’). For the overcast’ example, the lower atmosphere was 
essentially neutral (AT = I M - T water = -0.8°C) and produced high augmented reflectance and 
fractional coverage estimates. By the next day, the sea state under the influence of conditions of 
greater duration produced augmented reflectance and fractional coverage estimates that were 
roughly 1/3 those obtained during the onset of the cold front, despite the atmosphere being 
unstable (i.e. AT = T ar - T watCT = ~-4°C). The whitecaps visually appeared different to those in the 
‘overcast’ case with large breaking waves, although less frequent, and with more residual foam 
patches and foam streaks that were essentially absent in the ‘overcast’ case. Nevertheless, one 
would expect the opposite result with greater whitecap influence during unstable atmospheric 
conditions and less during neutral conditions. On the other hand, the augmented reflectance and 
fractional coverage for the ‘clear sky low sun’, which had similar whitecap and foam conditions to 
the ‘clear sky’ case, produced values roughly twice the ‘clear sky’ values due to the effects of the 
low sun. The dependence on sea state, wind variations with rising or falling trend, change in wind 
direction, minimal fetch and duration along with specular sun reflections, whitecap BRDF effects, 
and low sun angles are responsible for these differences in fractional coverage and augmented 
spectral reflectance for similar wind speeds. The data from the Gulf of Mexico was taken roughly 
1 10 km offshore which is an adequate fetch for 12 m/s wind speeds [/tars and Cardone , 1974], 
The estimate of the augmented reflectance from the ‘clear sky’ case, which was acquired under sea 
state conditions with adequate fetch and duration, and under illumination conditions that were not 
particularly extreme (sun elevation of 30°), produced similar results to the open ocean data taken in 
the Pacific. 

Another interesting aspect resulting from these examples is the spectral change in the 
upwelling signal during low sun angle. The spectral variation in estimates of the augmented 
reflectance contribution seems to be dependent on sea state and the whitecap morphology. In Fig. 
2(d) the decrease trom the 440 nm to the 860 nm value is -41 % from observation of many small 



whitecaps and very few residual and fragmented foam patches. From Fig. 3(d) the drop in 
augmented reflectance from 440 nm to 860 nm is -28% from observation of a few large whitecaps 
but mostly residual and patchy foam. We attribute the significant change in the spectral variation of 
whitecap reflectance from the ‘clear sky’ to the clear sky low sun’ examples to the nonlambertian 
nature of the whitecaps (the BRDF effect). This effect could also influence the magnitude the drop 
from 440 to 860 nm. 

Since the low sun has this influence on the augmented reflectance (and fractional coverage), 
the wind direction with respect to the sun becomes an important consideration. The longitudinal 
axis of waves (and whitecaps) are more or less perpendicular to the direction of the wind. In the 
‘clear sky low sun’ case, the wind is from the north-west at 325° with the sun at 80° azimuth with 
respect to true north. This results in whitecaps and the crests of waves being illuminated at -25° 
incident angle (azimuth). One might expect a maximum effect if the wind had been from the west 
at 260° or perhaps from the east at 80° azimuth and a minimal effect with the wind from the north at 
350° or south at 170° azimuth. 

The goal of this study is to determine the augmented reflectance contribution from whitecaps 
as viewed from satellites and, in particular, its spectral dependence. Useful satellite data of the 
ocean is taken around solar noon and from cloud free regions. The low sun angle examples (10° 
and 30° above the horizon) used here can be considered extreme conditions and are not frequently 
encountered. However, it brings attention the effect low sun angle has on the augmented 
reflectance contribution. This effect may have particular importance for higher latitude waters both 
north and south of the equator where, during the winter months or even in higher latitudes in the 
summer months, the sun maintains a low angle above the horizon at around midday. For higher 
sun elevation, the BRDF effect (although substantially reduced) may still exist. These high latitude 
waters are particularly prolific as regards phytoplankton blooms, triggered by seasonal nutrient rich 
cold upwellings. They are often windy areas with an abundance of whitecaps and a bias of the 
reflectance towards the yellow-red end of the visible spectrum, with low sun angle, brings 
attention to the 670 nm waveband near the phytoplankton fluorescence peak at 685 nm. The 



Moderate Resolution Imaging Spectroradiometer (MODIS) sensor [Salmonson, 1989] will be 
capable of estimating this phytoplankton fluorescence signal [Letelier and Abbott, 1996], This fact 
adds to the importance whitecaps may play in these higher latitudes. And of course, the high 
values at 670 and 860 nm at lower sun angles, from a combination of whitecaps and foam free 
waves, could complicate atmospheric correction algorithms unless adequately characterized. 

The spectral characteristics of the augmented reflectance found in both the Gulf of Mexico and 
from the Pacific Ocean seem to possess some spectral characteristics that are related to the water 
color properties in which they were formed. In Fig. 9 examples of the nonwhite water spectral 
profiles and their augmented reflectances from both water bodies are shown. Both examples were 
taken under overcast sky conditions. For Gulf of Mexico waters, the augmented reflectance from 
threshold interval 20 in the ‘overcast’ case is used and the augmented reflectance and the nonwhite 
water spectrum have distinctly lower 410 nm values than their respective 440 nm values. For the 
very blue and clear Pacific Ocean waters, the 410 nm value is higher than the 440 nm value in both 
the augmented reflectance and nonwhite water spectral profiles. 

The augmented reflectance profile in Fig. 9 for the Pacific Ocean has been complied by 
averaging the normalized augmented reflectance of each of the 17 data samples used in earlier 
modeling. The difference between each channel of the mean augmented reflectance profile shown 
and the mean ARC 410 . 550 is +2.9%, -6.0%, +6.1%, -3.1%, -2.4% and -21.6%. The standard 
deviation for each waveband over the 17 data points is typically less than, or about equal to, the 
size of the symbols in Fig. 9. For the augmented reflectance from the Gulf of Mexico the 
difference between each channel and its ARC 410 _ 550 is -4.5%, -1.8%, +7.0%, -0.8%, -2.1% and - 
42.2%. 

In both cases, the augmented reflectances are similar with the exception of the 410 and 860 
nm values. The 860 nm value is lower in the Gulf of Mexico example due to a different type of 
wave breaking (many small whitecaps suddenly appeared at the onset of a cold front yielding very 
high augmented reflectance and fractional coverage estimates). Although the 410 nm value of the 
augmented reflectance seems to be related to the water in which it was formed, the 5 10 nm value 



does not and would appear to be independent of water color type. A similarly low 410 nm value 
was observed in foam data generated by a ship’s wake in coastal waters off the coast of Southern 
California [Moore et al.. 1997]. The low 410 nm value may be influenced by the higher 
absorption properties of coastal waters at the shorter wavelengths, but the enhanced absorption 
alone cannot explain this effect. Initially it was thought that this coupling between the augmented 
reflectance and the water type in which the white water is formed might be explained by the 
inclusion of very thin foam types where a good proportion of the background water color is also 
measured. But this does not explain the apparent indifference of wavebands above 410 nm (and 
perhaps 440 nm) to the foam free background water color. 

Although most of our open ocean analysis utilized only data acquired during overcast 
conditions (to facilitate more accurate estimation of the background water reflectance threshold, and 
to eliminate sun glint, which is easy to mistake for thin foam), some data for partly cloudy and 
clear skies were examined. This partly cloudy/clear sky data possessed an augmented reflectance 
spectrum similar to the overcast cases, as well as similar values of ARC 410 . 550 (Fig. 5). Thus, we 
believe that our ARC results should be reasonably representative of clear sky situations, and thus 
useful for estimating the whitecap contribution to the reflectance measured by ocean color sensors. 

Conclusions 

To summarize, spectral characteristics found in the augmented reflectance can be influenced by 
a number of factors such as sky condition (e.g. bright horizon and overcast sky), solar elevation, 
sea state: type of wave breaking and direction of wave breaking (particularly at low solar 
elevation). The spectral profile of the augmented reflectance contribution may even be dependent 
on the color properties of the water body in which the whitecap forms, since the detection of very 
thin foam types include upwelling water-leaving radiance contributions. Despite the inaccuracies in 
determining the white water threshold, the spectral profile of the augmented reflectance remains 
consistent over a large number of reasonable threshold choices. Threshold determination is 
inherently complicated and subject to some error due to the overlap of white water and nonwhite 



water reflectance contributions. For most applications, any subtle spectral characteristics found in 
the estimates of the augmented reflectance, besides the distinctive lower 860 nm value, are 
overshadowed by the inaccuracy in determining the correct white water to nonwhite water 
threshold. 

We have determined the augmented reflectance contribution in the open ocean for the 410 - 
550 nm range to be 0.00031 at a wind speed of 8 m/s, 0.00058 at 10 m/s and 0.00096 at 12 m/s 
with a standard deviation of 0.000108. These results are somewhat lower than previously derived 
values of whitecap optical influence [Koepke, 1984; Gordon and Wang , 1994b]. Using Koepke’s 
effective reflectance of 22% and the fractional coverage determined by wind speed for water 
temperature above 14°C [Monahan and O’Muircheartaigh, 1980] yields values that are greater than 
4 times our values at 12 m/s and roughly 3 times at 8 m/s. Applying the fractional coverage 
relationship which takes into account atmospheric stability [Monahan and O’Muircheartaigh , 

1986], our results are still lower even when assuming a very stable atmosphere. Also the 860 nm 
value, which may be an indicator of the type of wave breaking responsible for the amount and 
depth of bubbles forced into the water column, is found to be 15%-20% lower than the optical 
influence of whitecaps in the visible for wind speeds 8-12 m/s in the open ocean, i.e., half of that 
observed by Frouin et al. [ 1996] in the surf zone. We have ascertained the augmented reflectance 
contribution at 860 to be 0.00026 at wind speed of 8 m/s, 0.00046 at 10 m/s and 0.00073 at 12 
m/s. This implies a lower contribution from whitecaps at 860 nm (-5.6 times lower for wind 
speeds at 12 m/s and -3.8 times lower for wind speeds of 8 m/s using the Koepke [ 1984] effective 
reflectance of 22% and the Monahan and O’Muircheartaigh [1980] wind speed to fractional 
coverage relationship). 

Koepke [1984] concluded that the optical influence of oceanic whitecaps was in fact less 
important than had been previously assumed since his effective reflectance reduced the optical 
influence of whitecaps by a factor of 2. Certainly from our direct measurements of the optical 
influence of whitecaps and foam in the 8 - 12 m/s wind speed range suggest that whitecap optical 
influence can be considered to be even less important than Koepke’s optimistic conclusion. 



Notation 


ARC 41l) _ 55l) 

ARC 670 

ARC^ 


AT 


U 

W 


mean augmented reflectance contribution from values at 410, 440, 510 and 550 nm. 

augmented reflectance contribution at 670 nm. 

augmented reflectance contribution at 860 nm. 

air temperature 

water surface temperature 

temperature differential (AT = l w - T wawr ); 

wind speed in meters per second 
fractional coverage of whitecaps and foam 
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Figure Captions 


Fig. 1(a). Time series of the reflectance of two whitecaps and their foam types passing in view of 
the radiometer. 

Fig. 1(b). Reflectance spectra of the two whitecaps of Fig. 1(a) reduced to reflectance intervals. 

Fig. 2(a). Reflectance spectra of -9,000 water surface reflectance samples from ‘overcast’ data 
reduced to reflectance intervals. 

Fig. 2(b). Histogram of the number of reflectance samples per reflectance interval for ‘overcast’ 
data. 

Fig. 2(c). Total reflectance of ‘overcast’ data (the sum of all weighted reflectance intervals). This 
is equivalent to the upwelling signal detected by a satellite assuming no atmospheric 
contribution. 

Fig. 2(d). Augmented spectral reflectance contribution for a range of reflectance interval threshold 
choices. The fractional coverage estimate for each augmented reflectance profile is shown. 
For higher threshold choices the contributions in the visible region equalize with the exception 
of the 410 nm and 860 nm contributions. 

Fig. 3(a). Reflectance spectra of - 10,000 water surface reflectance samples from ‘clear sky’ data 
(sun elevation of -30° above horizon) reduced to reflectance intervals. 

Fig. 3(b). Histogram of the number of reflectance samples per reflectance interval for ‘clear sky’ 
data. 



Fig. 3(c). Total retlectance of ‘clear sky’ data (the sum of all weighted reflectance intervals). This is 
equivalent to the upwelling signal detected by a satellite assuming no atmospheric contribution. 

Fig. 3(d). Augmented spectral reflectance contribution for a range of reflectance interval threshold 
choices. The fractional coverage estimate for each augmented reflectance profile is shown. 
Note the higher 670 nm and 860 nm values due to the influence of the longer wavelengths 
from a low sun. 

Fig. 4(a). Reflectance spectra of -7,000 water surface reflectance samples from ‘clear sky low 
sun’ data (sun elevation of -10° above horizon) reduced to reflectance intervals. 

Fig. 4(b). Histogram of the number of reflectance samples per reflectance interval for ‘clear sky low 
sun’ data. 

Fig. 4(c). Total reflectance of ‘clear sky low sun’ data (the sum of all weighted reflectance 
intervals). This is equivalent to the upwelling signal detected by a satellite assuming no 
atmospheric contribution. 

Fig. 4(d). Augmented spectral reflectance contribution for a range of reflectance interval threshold 
choices. The fractional coverage estimate for each augmented reflectance profile is shown. 
Note the much higher 670 nm and 860 nm and lower 410 nm values due to the influence of 
the longer wavelengths from a low sun. The 410 nm value is lower than the 860 nm value for 
all threshold choices. 

Fig. 5. Augmented reflectance contribution for different wind speeds from open ocean data taken 
in the equatorial Pacific. Data points of the augmented reflectance from different waveband 
regions and their power law relationships with wind speed are shown. ARC 4|() . 550 


is the mean 



augmented reflectance of wavebands 410, 440, 5 10 and 550 nm. ARC 670 and ARC, 60 are the 
augmented reflectance for the 670 and 860 nm wavebands respectively. ARC 410 . 550 PC / 

scattered is for data taken under partly cloudy and scattered conditions (solar elevations -50°) 

and not used in establishing power law relationship. 

Fig. 6. Fractional coverage vs. wind speed for same data used in Fig. 5. A power law 

relationship is applied although the fractional coverage values are more varied for a given wind 
speed. 

Fig. 7. Fractional coverage vs. mean augmented reflectance of 410, 440, 5 10 and 550 nm values. 
A linear fit is applied with a reasonable correlation of 0.694. Wind speed is shown for each 
data point. 

Fig. 8. Percentage drop in the augmented reflectance between ARC 410 . 55o and ARC S60 for wind 
speeds in the 8-12 m/s range. The linear fit indicates an increasing trend with higher wind 
speed. 

Fig. 9. Augmented and nonwhite water reflectance profiles normalized at 510 nm (not to scale) 
are shown for data taken in the Pacific Ocean and the Gulf of Mexico. The augmented 
reflectance for the Pacific Ocean is complied from all data presented in this paper used to 
estimate the optical influence of whitecaps. Note the similarities in augmented reflectance 
profiles regardless of the water color in which the whitecaps form. The 410 nm value seems to 
be the exception and mimics the nonwhite water reflectance in both cases. The 860 nm value, 
which may be used to indicate the type of wave breaking, is lower in the Gulf of Mexico data 
due to type of whitecaps formed by the sudden arrival of a cold front. 
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Abstract 


We provide an analysis of the influence of instrument polarization sensitivity on the radiance 
measured by space-borne ocean-color sensors. Simulated examples demonstrate the influence of 
polarization sensitivity on the retrieval of the water- leaving reflectance (p w ). A simple method for 
partially correcting for polarization sensitivity — replacing the linear polarization properties of the 
top-of-atmosphere reflectance by that from a Rayleigh-scattering atmosphere — is provided and its 
efficacy evaluated. It is shown that it improves p w retrievals as long as the polarization sensitivity 
of the instrument does not vary strongly from band to band. Of course, a complete polarization- 
sensitivity characterization of the ocean- color sensor is required to implement the correction. 
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Introduction 


The spectral reflectance of the ocean atmosphere system is modified by the concentration of 
marine phytoplankton, microscopic plants that constitute the first link in the marine food ch ain 1 
The component of the spectral reflectance that is due to photons backscattered out of the ocean (the 
water-leaving reflectance) is usually termed the ocean color. The flight of the Coastal Zone Color 
Scanner (CZCS) 2 ’ 3 on the satellite NTMBUS-7 proved the feasibility of measuring the concentration 
of the photosynthetic pigment chlorophyll a (a surrogate for the concentration of phytoplankton in 
the water) on a global scale. Based on the CZCS success, several follow-on ocean color missions 
have been planned, e.g., the sea-viewing wide-field-of-view sensor (SeaWiFS), 4 and the moderate 
resolution imaging spectroradiometer (MODIS). 5 

The contribution of the water- leaving reflectance, p w , to the reflectance of the ocean- atmosphere 
system, p t , is < 10%, the remainder being due to photons scattered by the atmosphere and re- 
flected from the sea surface. Extraction of p w from p t is referred to as atmospheric correction. 
As p w is a small component of p t , adequate atmospheric correction can only be effected if p t is 
accurately measured, e.g., a 5% error in p t translates into ~ a 50% error in p w in the blue in waters 
with low phytoplankton concentrations, and a larger error in waters with higher concentrations. 
The goal set for SeaWiFS and MODIS is the retrieval of p w in such waters with an uncert ain ty of 
±5%. 

In ocean color remote sensing, it has been assumed implicitly that the sensor is able to ac- 
curately measure (subject to calibration limitations 6 ) the radiance exiting the top of the atmo- 
sphere (TO A). However, the radiance reflected from the ocean- atmosphere system can be strongly 
polarized. 7 Since all radiometers display some sensitivity to the polarization state of the radiance 
they intend to measure, generally a biased measurement will be obtained. Although the ocean- 
viewing radiometers are generally designed to have low polarization sensitivity, e.g., for SeaWiFS 
and MODIS it was specified that the response vary by < 2% for all linear polarization states of 
the incident radiance, some may not meet the design requirements. Others, e.g., the Spectroscopic 
Imagers (SPIM) on the Midcourse Space Experiment (MSX), 8 were not specifically designed to 
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examine strongly polarized sources, and may be quite sensitive to polarization. Because of this, 
it is of interest to understand the influence of residual instrument polarization sensitivity on the 
retrieved products, and to devise a method to minimize it. This is the focus of the present paper. 

We begin by reviewing the measurements required to specify the polarization sensitivity of 
a radiometer, show how the “measured” TOA radiance depends on the polarization state of the 
light, and provide examples of the polarization properties of the TOA radiance. Next, we review 
the atmospheric correction algorithm and use it to derive the desired water-leaving radiance by 
operating it with the radiance measured by the sensor if the instrument’s polarization sensitivity 
is ignored. Finally, we propose a simple algorithm for making an approximate correction for the 
effects of residual instrument polarization sensitivity, and provide some examples of its efficacy. 

2. Polarization Sensitivity 


2A. Sensor polarization sensitivity 


We assume that the radiometer can be modeled as an optical system, which responds in a 
manner dependent on the state of polarization of the radiance to be measured, and a detector with 
a response that is independent of the state of polarization of the radiance. The radiance is specified 
by the column vector I, the Stokes vector , 9-12 



where 

/ = (EtEj + ErE;), 


Q = {E t El - E r E ;) 9 

u> 

U = (E t E’ r + E r Et), 

V = t (E t E‘ r - E r E;), 

Ei and E r are the components of the electric field in any two orthogonal directions normal to the 
direction of propagation, the superscript * indicates the complex conjugate, and the angle bracket 
denotes the average over time. The first element of the Stokes vector, /, is the radiance that would 
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be measured with a detector that is insensitive to the polarization state of the field. We refer to it 
here as the radiance. The polarization state of the radiation is determined by the other components 
of I, for example, the degree of polarization of the radiation is 


p = v^ 2 + ± v 2 


( 2 ) 


where 0 < P < 1. The limit P = 0 corresponds to a completely unpolarized radiation field, while 
P = 1 corresponds to a completely polarized radiation field. Radiation fields with intermediate 
values of P are partially polarized. Any field with a degree of polarization P can be represented as 
a linear combination of an unpolarized field of radiance (1 - P)I and a completely polarized field 
of radiance PI , i.e., 


I = 





Fields for which V = 0 are said to be linearly polarized. 


The action of the optical system on I is to produce a new Stokes vector I m , given by 


In* = MI, 


where M is a 4 x 4 matrix (the Mueller matrix). The measured radiance, / m , is the top element 
of the column vector I m . The trains form matrix M describes the action of the instrument on I. It 
follows from the action of the instrument on the fields, 



where E r aind Ei aire the electric field components of the beam in the directions perpendicular and 
parallel to a reference plaine, respectively. This reference plane is arbitrary. It is defined using the 
bawis vectors / and r, respectively, parallel and perpendicular to the reference plaine. The plane 
itself is formed by l and the direction of propagation of the radiance, f x /. Starting from the 
Eq. (3) and the definition of /, Q, U , and V in terms of the fields, the derivation of M from A is 
straightforward. 9 ’ 13 In what follows we will omit the explicit dependence of I amd M on wavelength 
A for simplicity. 
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Viewing the earth, the sensor responds to the Stokes vector I f exiting the top of the atmosphere 
(TOA). It is defined with respect to a reference plane determined by the propagation direction of 
the light (specified by the polar angles 6 , <j>) and the vertical, with ( t and f t , respectively parallel 
and perpendicular to this plane. However, the transformation matrix M is defined relative to a 
reference plane (l and r) fixed with respect to the instrument. If these two reference planes are not 
coincident, then a transformation (rotation) of the Stokes vector from one reference plane to the 
other has to be made. Letting ( t • i — cos a, we have 13 

I m = ME(a)I t , (4) 

where 

/I 0 0 0\ 

o / x _ 0 cos 2a sin 2a 0 | 

a — 0 —sin 2a cos 2a 0 I ’ 

VO 0 0 1/ 

with a measured clockwise from l t to t looking toward the source. Assuming that the detector itself 
is not sensitive to the polarization state of the light, Eq. (4) shows that the “measured” radiance 
I m is related to the components of the “true” stokes vector I e by 

7 m = M u It + Afi 2 (cos2a Q t + sin 2a U t ) + M 13 (-sin2a Q t + cos 2a U t ) + M 14 V t , (5) 

where Mu, Mu, Mi 3 and Mu are the elements of the first row of the matrix M. For the radiance 
backscattered to the top of the atmosphere, V t % 0 (refer to Subsection 2B), so I m in Eq. (5) can 
be rewritten as 

I m = Mul t + M \ 2 (cos 2 a Q t + sin 2a U t ) + Mij(— sin 2a Q t + cos 2 a U t ) (6) 

Clearly, the measured radiance I m is not It, I m depends on the characteristics of the instr um ent: 
(Mu, M 12 and M 13 ) and the characteristics of the radiation ( I t , Q t and U t ). Of course, it is 
desirable to have an ideal instrument with M being a unit matrix (M 12 = M l3 = 0 and M u = 
1), i.e., an instrument having no polarization sensitivity. However, this is not practical. All 
instruments, unavoidably, have some polarization sensitivity. On the other hand, if the degree of 
polarization of the radiance received by the instrument is zero ( Q t = U t = V t = 0), the measured 
radiance I m will be the true radiance I t . But, as the upwelling radiance we intend to measure is 
polarized, I m will never be I t . 
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In order to study the effect of the polarization sensitivity of the instrument, we need to know 
only three elements, A/ u , Mu and Mu, in the transform matrix M. We now describe a series 
of hypothetical experiments to determine Af u , Af 12 and Af 13 of the instrument. First, choose a 
reference plane for the instrument, with basis vectors / and r parallel and perpendicular to this 
reference plane, respectively, and t x / in the propagation direction of the beam. Then illumin ate 
the instrument with linearly polarized radiance (P — 1) of radiance Iq in the following ways [here, 
to conserve space, we write I as a row vector, I = (I, Q, U, V )]: 

(1) with the oscillation direction (plane) of the electric vector along /, so the incident Stokes 
vector is I in = 7 0 (1, 1,0,0), for which the instrument records a radiance and from Eq. (6) 
with a = 0 (since Mu, Mu and Mu are defined on this reference plane) we have 

I'i = (Mu + Mu)Iq\ 


(2) with the oscillation direction along r, i.e., Ii n = /o(l, — 1,0, 0), for which the instr ume nt 
records a radiance I+, 


h = (Mu ~ M 12 )Jo; 


(3) with its oscillation direction along the line which bisecting the angle between / and f, i.e., 
I»n = Jo(l,0, 1,0), for which the instrument records a radiance If r , 

If r = (Mu + Miz)Iq. 


From these experimental results, Mu, Mu and M l3 can be found easily: 


Mn = 


_ h + li 

2 /„ ’ 

h ~ h 

Ml ’ = ITT' 


Mn = -f- - Mn. 
-<0 


Considering that M\\ can be determined during instrument calibration using an unpolarized source, 
we use Mn = 1 in this paper. Upon defining q t and u t as 
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Eq. (6) can be rewritten 



Im =/*[! + m 12 (cos 2 a q t + sin 2a u t ) ■+■ m l3 (- sin 2a q t + cos 2a u t )] 


(7) 


with mi 2 and m l3 defined by 


Mi 


m 12 — 
mi 3 = 


12 


Mu’ 

M 13 


Mi 


u 


Even if Iq is unknown in the above experiment, we can still obtain the polarization sensitivity of 
instrument by finding m 12 and mi 3 from the following equations 

h - U 


m 12 — 


m 13 = 


/f + h ’ 

*Jir 

I i + /f 


- 1. 


2B. Polarization properties of I t 

To examine the influence of the instrument polarization sensitivity, we need to know the 
polarization characteristics of the radiance exiting the top of the atmosphere. For remote sensing 
of ocean color, in general I t can be written as 


It — Ir + I a 4* Ira + tl w , (8) 

where I r is the contribution from the Rayleigh scattering, I a is the contribution from the aerosol 
scattering, I ra is the contribution from the interaction between molecular and aerosol scattering, 
tl w is the water leaving radiance diffused transmitted (t) to the top of the atmosphere. This yields 

It = fr + fa + fra + tl w . 

The water-leaving component tl w is at most (in the blue) ^ 10% of the total I t , and considering 
the incident light field on the top of the atmosphere from the sun is completely unpolarized, 
the polarization state of I t is mainly governed by the atmospheric scattering contribution. For 
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Rayleigh scattering, its polarization state is well known and easy to determine. But for aerosols, 
their scattering properties, and the polarization state of the scattered light, are generally unknown, 
as they strongly depend on their particle size distributions and refractive indexes. 


Consider a spherical coordinate system with origin at the TOA, z axis directed downward, and 
x axis directed away from the sun (sun’s rays are in the x-z plane). Let the direction of propagation 
of a photon be specified by the polar and azimuth angles, 6 and #, in this system. The solar beam 
has ff = 0 0 and <p = 0. Photons exiting the TOA have 8 > 90°. Using r to represent the optical 
depth measured from the TOA, the propagation of I at a wavelength A is governed by the vector 
radiative transfer equation (VRTE) 

C os* g^M = _I(A,t,M) + u/ 0 jf R(a)P(A,T,^,^- + O, < A)R , (a')I(A,r,0',^)dn', (9) 

where P is the phase matrix, and R is the rotation matrix. As mentioned in Subsection 2A, the 
reference plane for I is defined by the direction of propagation and the z axis, i.e., the reference 
plane is perpendicular to the ocean surface. It should be noted that there are two rotations among 
three reference planes whenever a scattering occurs, 10 since the phase matrix is defined on the 
scattering plane (determined by the incident light and the scattered light). One must first rotate 
the incident reference plane to the scattering plane, then apply phase matrix P, and finally, rotate 
the scattering plane to the scattered reference plane. The boundary condition at the top of the 
atmosphere is 

I(O;0,<A) 6(cos0 - cos0 o )6(4>- <fo), 8 0 < t/2. 

The boundary condition at the (assumed flat) sea surface is given by 


I(n;8r,tr) = F (tfi)i(n ;*,&), 


( 10 ) 


where the Fresnel reflection Mueller matrix F is 

/p+(&i) p-(9i) 0 0 

-p la \ _ P-(9\) P+{&i) 0 0 

1 x) 0 0 p 3 s(0i) o 

VO 0 0 p 33 (0i). 

with 6 r = t - $i and <p r = 4>i- The factors p±(8i) and p 33 (0i) are given by 


ol 


/ cos 9{ - \/ n 2 - sin 2 9{ \ 2 

( sin 2 9{ - cos 9 x \/n 2 - sin 2 9 X 

Y ±l i 

V cos 9i + y/ n 2 - sin 2 9 X ) 

. V sin 2 9y 4* cos 9{ \/ n 2 - sin 2 0* 



( 11 ) 
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P 3 a( 0 ») = 


cos &i - \/tP~- sin 2 9, 


» \ / sin 6j 

i / V sin 2 9i 


— cos $i y n 2 - sin 2 


cos + \frp- sin 2 7 V sin 2 9 { + cos 9 { y/n 2 — sin 2 $i 
The aerosol scattering phase matrix has the following form for spherical particles: 


Pa(©) 


(P 11 

Pi 2 

0 

0 \ 

P 12 

Pi 1 

0 

0 1 

0 

0 

P 33 

P34 

V 0 

0 

-Pi 34 

P 33/ 


( 12 ) 


where P llt P l2 , P33 and P34 are functions of the scattering single 0 . For the aerosol models used in 
this study, the phase matrices are computed by using the Mie theory. The Rayleigh phase matrix, 
P r , is 

/ 1 + cos 2 0 — sin 2 0 


Pr(0) = 


16 x 


V 


0 0 

- sin 2 0 1 + cos 2 0 0 0 

0 0 2 cos 0 0 

0 0 0 2 cos 0 , 


( 13 ) 


where 0 is the scattering angle. Note that the depolarization factor 14 has been set to zero. 


The VRTE can be solved by the Successive Order Method and the Monte-Carlo Method, in 
a similar way used for solving the scalar radiative transfer equation (SRTE). However, instead of 
solving a single equation in the SRTE approximation, we have to solve four coupled equations. 


Although the upwelling radiance It at the TOA is not the result of just single scattering, we 
gain some understanding by examining the Rayleigh and aerosol single scattering contributions 
to I t . From single scattering theory, with the incident light on top of the atmosphere completely 
unpolarized, the upwelling radiance exiting the top of the atmosphere by direct molecular or aerosol 
scattering (i.e., without considering the reflection from the ocean surface) is given by 14 



F 0 cos 9 qljq 
cos 0 o - cos 9 


[l - exp( 


cos 9 cos $ 


/ Pi u 

_D_ol( p«(0) 

COS 0 O * I -^12(0 

V 0 


Pn(Q) 

Pi2(0) cos 2a 
(0) sin 2a 


( 14 ) 


where w 0 is the single scattering albedo, r! is the optical thickness of the molecules or aerosols, and 
0 is the scattering angle. It can be seen that there is no contribution to Vt from single Rayleigh 
scattering or single aerosol scattering. Observing that there is no P34 term in Rayleigh scattering 
phase matrix, there is no contribution to Vj from pure Rayleigh multiple scattering either. The 
term V t is the result of multiple aerosol scattering (including the aerosol and Rayleigh interaction 
terms) because of the existence of element P iA in the scattering phase matrix of aerosols, or the 
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result of single scattering followed or preceeded by reflection from the sea surface. Our simulation 
results agree with the computations of Kattawar et a/., 15 ’ 16 that V t % (V 0 + V ra ) -» 10 -3 / t . This is 
why we can assume V t ss 0 in Eq. (6). It means circular polarization effects are ignored. 

Prom Eq. (14) we can also see that, in the case of single scattering, Q and U (or the degree of 
polarization) are introduced mainly because of the existence of the scattering phase matrix element 
P 12 for Rayleigh and aerosol scattering. Figure 1 provides P n for Rayleigh scattering and aerosol 
scattering as a function of scattering angle 0. The aerosol models used in the figure are the Shettle 
and Fenn 17 Maritime and Tropospheric models at 80% relative humidity, labeled as M80 and T80 
respectively. It is observed from Figure 1 that P l2 for the aerosol is usually small compared to 
that for Rayleigh scattering, and that this element for Rayleigh scattering is a strong function of 
the scattering angle 0. In the backward directions (90° < 0 < 180°), when 0 is greater than 
~ 160°, the P \2 elements for Rayleigh and aerosol scattering are comparable. In contrast, when 
90° < 0 < 150° Rayleigh scattering has a very strong effect on the polarization state. Thus, when 
single scattering dominates, the degree of polarization will be small for luge scattering angle 0 
(0 > 160°), while when 0 is less than 150°, strong polarization effects will be encountered, caused 
largely by Rayleigh scattering. 

By employing the aerosol models described in Gordon and Wang, 18 the Stokes vector I t can be 
computed by solving the VRTE. Using a two-layer atmosphere model with molecules confined in 
the upper layer, and aerosols confined in the lower layer, bounded by a flat Fresnel- reflecting ocean 
surface, the VRTE was solved using a Monte-Carlo code to provide a pseudo radiance vector It 
received by the sensor at the top of the atmosphere. The simulations we present were carried out 
for the M80 and T80 aerosol models in three wavebands (443,765, and 865 nm), and the aerosol 
optical thickness for A = 865 nm was taken to be 0.2 [r a (865) = 0.2], 

As the error in I m [Eq. (7)] is determined by q t and u t , the ratios of Q t / I t and Ut/It computed 
for M80 at three wavelengths are plotted in Figure 2 as a function of the viewing angle of the 
sensor, 0 V = t — 9, for <fi = 90°. The T80 aerosol model yields similar curves. It should be noted 
that U t changes sign on crossing the principal plane {4> = 0). In the figure, we plot U t for the 
side of the principal plane for which it is negative. In contrast, Q t is symmetric with respect to 
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reflection across the principal plane. These computations show that (1) the degree of polarization, 
P, generally increases with increasing 9 0 and 0 V ; (2) P can be as large as — 0.5-0. 6; (3) Q t and 
U t both contribute significantly to P (Q t for small 0 V and U t for large 0 V ); and the polarization 
properties at 765 and 865 am are similar, but show significantly less polarization than at 443 nm. 


These computations can be used to obtain a coarse estimate of the error in I t induced by the 
polarization sensitivity. Choose l and r so that m 13 = 0. Then, 

A I t _Im~ It 


It 


It 


— m \ 2 ( cos 2a -f u t sin 2a) . 


It is easy to show that the maximum value of (q t cos 2a + u t sin 2a) is P* = y/q \ + u \ . Thus 

AI t 

< m l2 P t . 


As Pi can be as large as 0.5-0. 6, we see that the error in I t is potentially as large as ~ 0.5mi2- For 
mi 2 < 0.02, the error is at most 1%. In contrast if mi 2 = 0.10 the error could be as large as 5%. 
Gordon 19 has shown that errors of 1% can be tolerated in atmospheric correction as long as the 
error has the same sign throughout the spectrum, however, errors — 5% cannot be tolerated. 


3. Effects of sensor polarization sensitivity on atmospheric correction 

In the absence of strongly- absorbing aerosols and instrument polarization sensitivity, the Gor- 
don and Wang atmospheric correction algorithm works well. 19 The error in the retrieved water- 
leaving reflectance at 443nm A/?(443) = £(443) Ap w (443), where the reflectance p w = t/^/Pq cos 0q, 
is about ~ ±0.001 to ±0.002. This meets the requirements of MODIS and SeaWiFS. However, since 
the instrument has a residual polarization sensitivity, and the upwelling radiance I t we intend to 
measure is polarized, we cannot be provided with the true I t . Instead, we will have the biased 
I m . Therefore, it is necessary to assess the influence of polarization sensitivity on the Gordon and 
Wang atmospheric correction algorithm. 

We employed simulated I m pseudo data to study the influence of the polarization sensitivity 
of the instrument on the performance of the Gordon and Wang atmospheric correction algorithm. 
Briefly, l t was computed for the M80 and T80 aerosol models with aerosol optical thickness at 865 
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nm, r a (865) = 0.2. This value of r 0 (865) is somewhat higher than observed over regions with a pure 
maritime atmosphere, i.e., not subjected to anthropogenic aerosol or mineral dust transported over 
the oceans. 20 22 The water-leaving radiance was taken to be zero. The simulations were carried out 
for 0 O = 0, 20°, 40°, and 60°, both at the center, 9 V ss 1°, and the edge, 9 V ^ 45°, of the SeaWiFS 
scan with <j> = 90°. The combination 0 O = 0 and 9 V ^ 1° is omitted because it would be near the 
center of the sun’s glitter pattern. These seven geometries approximately cover the full range of 
sun- view geometries encountered in progressing along a polar orbit from the equator to a latitude 
of 60° at the equinox. 

The Stokes vector I t was then introduced into Eq. (7) and I m computed given the polarization 
sensitivity of the sensor. A rotation of the reference plane used to define m l2 and m 13 does not 
change \/ m* 2 + mj 3 , it only reallocates polarization sensitivity between m 12 and m 13 , i.e., it simply 
changes the definition of the angle a. Thus, for simplicity we set m 13 equal to zero, and present 
the results of our study as a function of a. 

The Gordon and Wang algorithm 18 uses the SeaWiFS bands at 765 and 865 nm (where p w = 
0, except in turbid coastal water) to provide atmospheric correction for the visible. Here, we 
examine the error in the water-leaving reflectance, tAp w = Ap, at 443 nm. To demonstrate the 
effect of the polarization sensitivity, we examined six sets of values of m 12 in the three wavebands 
443, 765 and 865 nm. The six sets of mu [written (m l2 (443), mx 2 (765), mi 2 (865))] were mi 2 = 
(0.10,0.10,0.10), m 12 - (0.05,0.05,0.05), m l2 = (0.05,0.03,0.03), m 12 = (0.02,0.02,0.02), m 12 = 
(0.05,-0.05,-0.05), and m l2 = (0.02,-0.02,-0.02). As mentioned above m 13 = (0,0,0) for 
each set. Note that setting m 13 to zero at all three wavebands, implies that when the plane of 
polarization of radiation (of constant radiance) entering the sensor is rotated through 360°, the 
sinusoidal responses (in 2x the rotation angle) of the detectors for the three bands will be either 
in phase or exactly out of phase with each other. This is not a necessity for the development we 
present, it is used only to simplify the analysis. 

Sample results for the errors in the retrieved water-leaving reflectance at 443 nm are presented 
in Figure 3, for viewing at the scan edge with <f> = 90°, as a function of a. For this fixed viewing 
geometry, a would be constant; however here, varying a is identical to varying the fraction of 
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the polarization sensitivity allocated to m 12 and m 13 . The solid line connecting the solid dots 
on the figure provides the error in tp w in the absence of polarization sensitivity, i.e., the result 
of operating the algorithm with the correct input, I t not 7 m . The simulation results in Figure 3 
suggest that (1) large instrument m 12 causes large errors when polarization state is considered; (2) 
for mi 2 as small as 0.02, the polarization sensitivity effects are not a problem in most sun-viewing 
geometries (error < ±0.002) as long as the m l2 has the same sign in all band; (3) retrieval errors 
for mx 2 = (0.05,0.05,0.05) are about same as that for m l2 = (0.05,0.03,0.03), and retrieval errors 
for 771x2 = (0.05, -0.05, -0.05) sire about same as that for 77i l2 = (0.10,0.10,0.10), suggesting that 
when the tt^’s are wavelength dependent, the retrieval errors are larger than that when they are 
wavelength independent. 

It is clear that the performance of the Gordon and Wang algorithm is degraded in the presence 
of sensor polarization sensitivity. Thus, a method to remove the errors induced by the instrument 
polarization sensitivity is required. 

4. Removal of instrument polarization sensitivity 

To completely remove the instrument polarization sensitivity, the polarization properties of the 
upwelling radiance vector I t are needed. As aerosol scattering is highly variable, measuring only 
the upwelling radiance I t can not provide any information regarding its polarization characteristics. 
Fortunately, the polarization state of the upwelling radiance I t is mainly determined by Rayleigh 
scattering (Figure 1). Figure 4 provides a comparison of Q /I and U/I between Rayleigh scattering 
(without aerosols) and total scattering (Rayleigh plus aerosols) for 0$ — 60°. Associated with 
Q / 7 and U / 7 , the subscript w r” refers to Rayleigh scattering and for total scattering. Figure 4 
shows that, generally in the blue Q r /7 r is close to Qt/It> and U r jl r is close to U t / I t , with \Q r /I r \ > 
\Qt/It\ and \U r /I r \ > \U t /It\- This can be explained as the aerosol-scattering contribution to Q t 
and U t is generally smaller than its contribution to 7 t , which means the existence of the aerosol 
reduces the degree of polarization, due to pure Rayleigh scattering, in the upwelling radiance 
vector It. For the short wavelength (Figure 4a), the Rayleigh scattering dominates, and there is 
little difference between Q r /7 r and Q t j I t , and between U r / I r and Ut/ I t . For the long wavelengths, 
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aerosol scattering significantly reduces the degree of polarization of the upwelling radiance. 

Observing that a significant portion of q t and u t is contributed by Rayleigh scattering, and 
noting that the Rayleigh scattering properties of the air are well known, we would expect that a 
significant amount of the polarization-induced error in ^^(443) could be removed by replacing q t 
and u t in Eq. (7) by their Rayleigh-scattering counterparts, Q r /Ir and U T /I T . When q r close to 
q t and u r close to u t , the corrected radiance should be close to the true radiance I t . Using the 
upwelling radiance at the top of the atmosphere corrected in this m ann er, instead of the biased 7 m , 
we applied Gordon and Wang atmospheric correction algorithm to the same pseudo data used in 
Figure 3. The retrieved errors in water-leaving reflectance at 443 nm are presented in Figure 5. It 
can be observed from the simulations that: (1) generally, the polarization correction works better 
for smaller m 12 ; (2) when mi 2 is not dependent on wavelength, it works very well, even for m !2 
as large as (0.1, 0.1, 0.1); (3) when m 12 is dependent on wavelength, it does not work as well as 
in (1); (4) when values of m i 2 have different signs in different wavebands, even for m l2 as small as 
2%, it cannot perform well (Figures 5e and 5f). Similar results were obtained for the T80 aerosol 
model with r a = 0.2, although the residual error in tp w was larger. For the purpose of removing the 
effects of polarization sensitivity, these simulations show the importance of designing instruments 
for which does not depend significantly on wavelength. 

We have tried to improve on this polarization-sensitivity correction by accounting for the 
presence of aerosols, which cause the difference between the q r -u r pair and the qt~u t pair (Figure 
4). However, as the polarization properties of the aerosol cannot be known prior to atmospheric 
correction, some assumption must be made in this regard. Figure 1 suggests that a reasonable 
assumption would be that aerosol single scattering completely depolarizes incident radiance, i.e., 
that the only nonzero element of P a (0) is at the top left of the matrix. In the single-scattering 
approximation to the TO A radiance, this would correspond to replacing the actual q t by Q r /h 
and u t by U r / I t . Such a replacement would effect a first-order correction for the depolarization by 
aerosols. As I t is unknown, we tried replacing I t by / m , the radiance measured without regard for 
the polarization sensitivity, i.e., in Eq. (7) we replaced q t and U* by Q r / Im and U r /Im , respectively. 
The results were disappointing, there was no improvement in the polarization correction, and in 
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many cases the error after correction was larger than that shown in Figure 5 (but less than that in 
Figure 3). Thus, we rejected this method in favor of replacing q t and u t in Eq. (7) by q r and u r , 
respectively. 

5. Concluding remarks 

We have presented simulations demonstrating the effect of sensor polarization sensitivity on 
the atmospheric correction of ocean color sensors (Figure 3). In addition, we provided a simple 
method — substituting the for polarization properties of I* those of a pure Rayleigh-scattering 
atmosphere — for partially correcting the error induced by polarization sensitivity. This correction 
method was shown to be effective (even for relatively large polarization sensitivities) as long as the 
polarization sensitivity of the instrument does not vary strongly from band to band (Figures 5a-d). 
An attractive feature of this polarization correction is the simplicity of its implementation as part of 
the overall atmospheric correction algorithm. The possible error in the retrieved values of tp w with 
out polarization correction (Figure 3), underscores the importance of a complete characterization 
of the polarization sensitivity of an ocean- color sensor prior to launch so that the polarization- 
sensitivity correction can be applied. 
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Figure Captions 


Figure 1. Scattering phase matrix element P u for aerosols at 443, 765 and 865 nm, and Rayleigh 
(solid line): (a) M80 aerosol model; (b) T80 aerosol model. 

Figure 2. Ratios Q t /I t and U t /I t as a function of 0„ for the M80 aerosol model with <j> = 90°: (a) 
0 O = 20°; (b) 0 O = 60°. 


Figure 3. Retrieved error tAp w at 443 nm as a function of 0o for the M80 aerosol model with 
r 0 (865) = 0.2, in the presence of polarization sensitivity, for viewing at the edge of the scan 
with (a) mi 2 = (0.10,0.10,0.10); (b) m n = (0.05,0.05,0.05); (c) m n = (0.05,0.03,0.03); (d) 
mi 2 = (0.02, 0.02, 0.02); (e) m n = (0.05, -0.05, -0.05); (f) m l2 = (0.02, -0.02, -0.02). 

Figure 4. Ratios Qjl and U/I for Rayleigh scattering and total scattering as a function of viewing 
zenith angle 0„ for the M80 aerosol model with r„(865) = 0.2 and 0 O = 60°: (a) A = 443 nm; (b) 
A = 765 nm; (c) A = 865 nm. 

Figure 5 Retrieved error tAp w at 443 nm as a function of 0o for the M80 aerosol model with 
r 0 (865) = 0.2, after applying the removal algorithm, for viewing at the edge of the scan with 
(a) mij = (0.10,0.10,0.10); (b) m n = (0.05,0.05,0.05); (c) m l2 = (0.05,0.03,0.03); (d) m u = 
(0.02, 0.02, 0.02); (e) m^ = (0.05, -0.05, -0.05); (f) m l2 = (0.02, -0.02, -0.02). 
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Figure 2. Ratios Q t / It and Ut/ It as a function of 0„ for the M80 aerosol model with <f> = 90°: (a) 
9 o = 20°; (b) 0 O = 60°. 
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Figure 3. Retrieved error tAp w at 443 nm as a function of &o for the M80 aerosol model with 
T a (865) = 0.2, in the presence of polarization sensitivity, for viewing at the edge of the scan 
with (a) m 12 = (0.10,0.10,0.10); (b) m l2 = (0.05,0.05,0.05); (c) m l2 = (0.05,0.03,0.03); (d) 
m 12 = (0.02,0.02,0.02); (e) m 12 = (0.05, -0.05, -0.05); (f) m t2 = (0.02, -0.02, -0.02). 
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Figure 4. Ratios Q / 1 and U j I for Rayleigh scattering and total scattering as a function of viewing 
zenith angle 9 V for the M80 aerosol model with r a (865) = 0.2 and 9 0 = 60°: (a) A = 443 nm; (b) 
A = 765 nm; (c) A = 865 nm. 
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Figure 5 Retrieved error tAp w at 443 am as a function of 9q for the M80 aerosol model with 
r tt (865) = 0.2, after applying the removed algorithm, for viewing at the edge of the scan with 
(a) m 12 = (0.10,0.10,0.10); (b) m l2 = (0.05,0.05,0.05); (c) m 12 = (0.05,0.03,0.03); (d) m l2 = 
(0.02,0.02,0.02); (e) m l2 = (0.05, -0.05, -0.05); (f) m l2 = (0.02, -0.02, -0.02). 







